colpitts4w2 .m 



i 



% Program to Calculate Electrical Component Values for Oscillator 
System of Coil #1 

syms t CI C2 C3 LI Rl Rs Rd Rll Vbias Vddl Vdd Vdds IDss Vss Vg II 

12 13 14 Xcl Xc2 Xc3 Xll Gm f pi; 
syms Amplitude; 
syms GmO Vgs Vgoff; 

syms Ml M2 Bl Result 1 M3 M4 B2 Result2 signl sign2; 

syms Fcx; 

pi = 3.14159; 

CI - 0.000000000680*1000; 

C2 - 0.000000000075*1000; 

C3 = 0.000000000047*1000; 

LI = 0.0000013*1000; 

Rl = 1000000*1000; 

Rd = 6800*1000; 

f = 1/ (2*pi* (L1*C1*C2/ (CI + C2)) A (l/2)); 

Rll = 2*pi*f*Ll; 

Vdd = 12*1000; 

Vg = sin (2*pi*f ) ; 

Xcl = 1/ (2*pi*f*Cl) ; 

Xc2 = 1/ (2*pi*f*C2) ; 

Xc3 = 1/ (2*pi*f*C3) ; 

Xll = 2*pi*f*Ll; 

IDss = 0.003*1000; 

Vgoff = 1.5*1000; 

GmO = -2*IDss/Vgof f ; 

Vgs = -2*1000; 

Gm = (Vgoff / (Gm0*Vgs) ) + 1; 

Rs = Rd - (1/Gm) ; 

Vss = IDss*Rs; 

Vddl = Gm*Rd*Vg; 

Amplitude = 0.0039; 

Vbias - 0.0039; 

Ml = [(Xll + Xcl + Xc2 + Rll) (-l*Xcl) (-l*Xc2) (0.0);]; 

Ml = [Ml; (-l*Xcl) (Xcl + Rs) (-l*Rs) (0.0);]; 

Ml = [Ml; (-l*Xc2) (-l*Rs) (Xc2 + Rs + Xc3) (0.0);]; 

Ml = [Ml; (0.0) (0.0) (0.0) (Rl + Rd) ; ] ; 

Bl = [ (Vg) ; (Vss); (-l*Vss + Vdds); (Vdd + Vbias);]; 

M2 = [Bl(l) (-l*Xcl) (-l*Xc2) (0.0);]; 

M2 = [M2; Bl(2) (Xcl + Rs) (-l*Rs) (0.0);]; 

M2 = [M2; Bl(3) (-l*Rs) (Xc2 + Rs + Xc3) (0.0);]; 

M2 = [M2; Bl(4) (0.0) (0.0) (Rl + Rd) ; ] ; 

Resultl = abs (det (M2) /det (Ml) ) ; 

% 

% Now, you want to design a circuit that has Fc times the freq. of 
circuit 1 

% and Ac times the amp. of circuit 1 
syms Fc Ac Acl Xcla Xc2a Rsa; 
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Vg = Acl*sin(2*pi*f*Fc) ; 

% Xc2a is the reactance for the second capacitor adjusted for Fc a 
nd Acl 

Fcx = ( (2*2*pi*pi) / (Fc*Fc*f*f) )* (C1*0. 000000001 + C2*0 . 000000001 ) / 
(C1*0. 00000000 1*C2*0.000000001*L1) ; 
Xcla = l/(2*pi*f*Cl*0.000000001*Fcx) ; 
Xc2a = 1/ (2*pi*f*C2*0.000000001*Fcx) ; 
Rsa = Rd*Acl - (1/Gm) ; 

M3 = [(Xll*Acl + Xcla*Acl + Xc2a + Rll*Acl) (-l*Xcla*Acl) (-l*Xc2a 
) (0.0);]; 

M3 = [M3; (-l*Xcla*Acl) (Xcla*Acl + Rsa) (-l*Rsa) (0.0);]; 
M3 = [M3; (-l*Xc2a) (-l*Rsa) (Xc2a + Rsa + Xc3*Acl) (0.0);]; 
M3 = [M3; (0.0) (0.0) (0.0) (Rl*Acl + Rd*Acl);]; 

B2 = [(Vg*Acl); (IDss*Rsa); (-1* (IDss*Rsa) + Vdds*Acl); (Vdd*Acl + 

Vbias*Acl) ; ] ; 
M4 = [B2(l) (-l*Xcla*Acl) (-l*Xc2a) (0.0);]; 
M4 = [M4; B2(2) (Xcla*Acl + Rsa) (-l*Rsa) (0.0);]; 
M4 = [M4; B2(3) (-l*Rsa) (Xc2a + Rsa + Xc3*Acl) (0.0);]; 
M4 = [M4; B2(4) (0.0) (0.0) (Rl*Acl + Rd*Acl);]; 
Result2 = abs (det (M4) /det (M3) ) ; 
signl = abs (det (M2) /det (Ml) ) / (det (M2) /det (Ml) ) ; 
sign2 = abs (det (M4 ) /det (M3) ) / (det (M4 ) /det (M3) ) ; 
cw2outl ; 

fid = f open ( f f : /drawl/colpitts4w2 . txt f , 1 w 1 ) ; 

syms DesiredAmp DesiredAmpl NewAmp NewAmp2 pAcl Aclx2 Closest Clos 

est Acl NAexp; 

DesiredAmp = 2.564102564; 

Fc = 1000000000/(16983891.6318); 

Closest = 0.0; 

ClosestAcl = 0.0; 

NAexp = abs (eval ( (Result2/sin (2*pi*f*Fc) ) / (Resultl/sin (2*pi*f ) ) ) ) ; 
for Acl = 1:1000, 

NewAmp = abs (eval (NAexp) ) ; 

if abs (NewAmp - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = Acl; 
end 

if abs (NewAmp - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/10; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/10; 
end 
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if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/100; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs(NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/100; 
end 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/1000; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs(NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/1000; 
end 

if abs(NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/10000; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/10000; 
end 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

fprintf (fid, 1 ;\n f ) ; 

fprintf (fid, 1 ; Program colpitts4w2 .m - (c) 2003 Leonard Gojer\n'); 
fprintf (fid, 1 ^-Computation of Coil Component Parameters for Drawl\n 
f ); 

fprintf ( fid, 1 ; \n 1 ) ; 

fprintf (fid, 1 ; Initial Frequency and Amplitude Change To Set Up Dra 
wl\n f ); 

fprintf (fid, 1 ; from Colpitts example in textbookAn 1 ) ; 
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■Fnr*i ri 4~ "F 
JL \J L J_ 1 1 U J_ 


f -Ft H 


^ J J . 1DI / KcSUlLl ] / 






iprinii 


, 1 10/ 


\n ; , 






rpr in l i 


f -F -1 /-J 


/ beconoary Kesuit ^! Ainplituae 




) ; 


■f- y-N y» "1 T"» "f" P 

iprinLi 
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'\n') ; 






rpr m t i 


f -F -i ^ 

, r la, 


-Computed Acl\n'); 






f print f 


' -F A si 
, 1 10/ 


'; (Component Amplification Factor) 




) ; 


rprinti 


( -F -i 

, r iO/ 


•%35.16f ' ,Acl) ; 






rpr in l r 


, r iO/ 


' \n ' ) ; 






"F V\ V* T 4— -p 

rprin r r 


f -F -! 

^ r iO/ 


'; Check of Amplification Factor Desired\n'); 




iprinLi 


, r iO/ 


'; (Result 2/Result 1) 




) ; 


■F y\ v* "i r"» 4- "F 

rpr in l. r 


, r iQ/ 


' %35.16f • ,eval (Result2/Resultl) ) ; 






P *~» v~ "i t*» 4— -F 

iprint r 


, r iO/ 


• \n ' ) ; 






fprintf 


[fid, 


' ; \n • ) ; 






fprintf 


[fid, 


1 (defun crdrawblks ( ) ' ) ; 






fprintf 


[fid, 


'\n' ) ; 






fprintf 


[fid, 


• (progn' ) ; 






fprintf 


[fid, 


•\n' ) ; 







DesiredAmpl = DesiredAmp; 

syms sacm; 

sacm = size(acm); 

for il = l:sacm 

DesiredAmp = acm(il) * DesiredAmpl; 

Fc = il*1000000000/ (16983891,6318) ; 

Closest = 0.0; 

for Acl = 1:1000/ 

NewAmp = abs (eval (NAexp) ) ; 

if abs (NewAmp - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = Acl; 
end 

if abs (NewAmp - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5/ 

Acl = pAcl + Aclx2/10; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 
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ClosestAcl = pAcl + Aclx2/10; 
end 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/100; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) 



ClosestAcl = pAcl + Aclx2/100; 
end 

if abs(NewAmp2 - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 
pAcl = Acl; 
for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/1000; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) 



< abs (Closest 

< abs (Closest 



DesiredAmp) 
DesiredAmp) 



ClosestAcl = pAcl + Aclx2/1000; 
end 

if abs (NewAmp2 - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 
pAcl = Acl; 
for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/10000; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) 



< abs (Closest 

< abs (Closest 



DesiredAmp) 
DesiredAmp) 



ClosestAcl = pAcl + Aclx2/10000; 
end 

if abs(NewAmp2 - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

fprintf (fid, 1 ;\n f ) ; 

f print f (fid, 1 ; Desired Frequency 

fprintf (fid, 1 %35.16f 1 , round (f*Fc) ) ; 

fprintf (fid, ' \n f ) ; 

fprintf (fid, 1 ; Desired Amplification 



< abs (Closest 

< abs (Closest 



DesiredAmp) 
DesiredAmp) 



= f ); 
= ■); 
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fprintf (fid, 1 %35 . 16f 1 , DesiredAmp) ; 
fprintf (fid, 1 \n 1 ) ; 

fprintf (fid, 1 ; Initial Result 1 Amplitude 
fprintf (fid, f %35.16f ' ,Resultl) ; 
fprintf (fid, 1 \n ' ) ; 

fprintf (fid, f ; Secondary Result 2 Amplitude 

fprintf (fid, 1 %35. 16f 1 ,eval (Result2) ) ; 

fprintf (fid, 1 \n' ) ; 

fprintf (fid, 1 ; Computed Acl\n ? ); 

fprintf (fid, 1 ; (Component Amplification Factor) 

fprintf (fid, '%35.16f ' ,Acl) ; 

fprintf (fid, 1 \n 1 ) ; 

fprintf (fid, ' ;Check of Amplification Factor Desired\n' 
fprintf (fid, 1 ; (Result 2/Result 1) 
fprintf (fid, '%35.16f 1 , eval (Result2/Resultl ) ) ; 
fprintf (fid, 1 \n* ) ; 

fprintf (fid, 1 ; Sign 1 (sign of Result 1) 
fprintf (fid, 1 %35.16f 1 ,signl) ; 
fprintf (fid, ' \n' ) ; 

fprintf (fid, 1 ; Sign 2 (sign of Result 2) 

fprintf (fid, 1 %35.16f 1 ,eval(sign2) ) ; 

fprintf (fid, 1 \n 1 ) ; 

fprintf (fid, 1 ;L1 

fprintf (fid, 1 %35 . 16f ' , Ll*Acl ) ; 

fprintf (fid, 1 \n' ) ; 

fprintf (fid, 1 ;C1 

fprintf (fid, 1 %35. 16f f , eval (Cl*Fcx/Acl) ) ; 
fprintf (fid, 1 \n 1 ) ; 
fprintf (fid, 1 ;C2 

fprintf (fid, ' %35. 16f 1 ,eval (C2*Fcx/Acl) ) ; 

fprintf (fid, 'Xn 1 ) ; 

fprintf (fid, 1 ;C3 

fprintf (fid, 1 %35 . 16f ' , C3/Acl ) ; 

fprintf (fid, 1 \n 1 ) ; 

fprintf (fid, ' ;Rd 

fprintf (fid, f %35 . 16f f , Rd*Acl ) ; 

fprintf (fid, 1 \n' ) ; 

fprintf (fid, 1 ;R1 

fprintf (fid, ' %35.16f 1 ,Rl*Acl) ; 

fprintf (fid, ' \n 1 ) ; 

fprintf (fid, 1 ;R11 

fprintf (fid, 1 %35.16f 1 ,Rll*Acl) ; 

fprintf (fid, f \n ! ) ; 

fprintf (f id, 1 ; \n' ) ; 

fprintf (fid, 1 ; \n' ) ; 

fprintf (fid, 1 (setq count f ); 

fprintf (fid, , %22.16f 1 ,il) ; 

fprintf (fid, 1 ) \n' ) ; 
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fprintf 


(fid, 




(setq DesFreq 1 ) ; 


fprintf 


k fid, 


1 %22 . 


lbf ' , round ( f *Fc) ) ; 


fprintf 


( fid, 


' ) \n' 




fprintf 


(fid, 




(setq DesAmp 1 ) ; 


fprintf 


, fid, 


' %22 . 


16f ' , DesiredAmp) ; 


fprintf 


(fid, 


' ) \n' 


) r 


fprintf 


( fid, 




(setq InitlAmp f ); 


fprintf 


, fid, 


1 %22 . 


lbf, Result! ) ; 


fprintf 


(fid, 


f \ V ■ 

1 ) \n' 


) J 


rprintf 


, fid, 




(setq Sec2Amp '); 


rprinti 


, rid, 




lbf , eval (Re suit 2 ) ) ; 


fprintf 


(fid, 


f ) \n ! 




fprintf 


( fid, 




(setq CompAcl 1 ) ; 


fprintf 


, fid, 


1 Q. O O 

1 %22 . 


lbf 1 , Acl) ; 


fprintf 


( fid, 


1 ) \n f 


) r 


fprintf 


(fid, 




(setq CheckAmp 1 ); 


fprintf 


.fid, 


' %22 . 


lbf 1 , eval (Result2/Resultl) ) ; 


fprintf 


!f id, 


1 ) \n 1 


/ r 


fprintf 


( fid, 




(setq Signl 1 ) ; 


fprintf 


( fid, 


' %22 . 


16f 1 , signl) ; 


fprintf 


(fid, 


' ) \n f 




fprintf 


(fid, 




(setq Sign2 1 ) ; 


fprintf 


(fid, 


1 o o o 

1 %22 . 


lbf 1 , eval (sign2) ) ; 


fprintf 


(fid, 


1 ) \n' 


/ / 


fprintf 


( fid, 




(setq LI 1 ) ; 


fprintf 


t fid, 




lbf , Ll*Acl) ; 


fprintf 


( fid, 






fprintf 


( fid, 




(setq CI 1 ) ; 


fprintf 


( fid, 


1 %22 . 


16f 1 ,eval (Cl*Fcx/Acl) ) ; 


fprintf 


( fid, 


1 ) \n' 


) f 


fprintf 


(fid, 




(setq C2 1 ) ; 


fprintf 


( fid, 


1 %22 . 


lbf 1 , eval (C2*Fcx/Acl ) ) ; 


fprintf 


( fid, 


' ) \n f 


) ) 


fprintf 


(fid, 




( setq C3 1 ) ; 


fprintf 


[fid, 


' %22 . 


16f 1 , C3/Acl) ; 


fprintf 


(fid, 


1 ) \n' 




fprintf 


(fid, 




(setq Rd 1 ) ; 


fprintf 


(fid, 


1 %22 . 


16f 1 , Rd*Acl) ; 


fprintf 


(fid, 


1 ) \n ! 




fprintf 


(fid, 




(setq Rl 1 ) ; 


fprintf 


( fid, 


1 o o o 

' %22 . 


lbf 1 , Rl*Acl) ; 


fprintf 


(fid, 


•)\n' 




fprintf 


(fid, 




(setq Rll ' ) ; 


fprintf 


(fid, 


■%22. 


16f ' ,Rll*Acl) ; 


fprintf 


(fid, 


■ )\n' 


\ • 
I r 


fprintf 


(fid, 




(insblkd count DesFreq DesAm 



CompAcl CheckAmp S 
fprintf (fid, 1 LI 



ignl Sign2 1 ) ; 
CI C2 C3 Rd Rl Rll) \n' ) ; 
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fprintf (fid, ' ;\n' ) ; 
end 

fprintf (fid, ' )'); 
fprintf (fid, ' \n' ) ; 
fprintf (fid, • ) ' ) ; 
fprintf (fid, • \n ' ) ; 
f close (fid) ; 
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% Program to Calculate Electrical Component Values for Oscillator 
System of Coil #2 

syms t CI C2 C3 LI Rl Rs Rd Rll Vbias Vddl Vdd Vdds IDss Vss Vg II 

12 13 14 Xcl Xc2 Xc3 Xll Gm f pi; 
syms Amplitude; 
syms GmO Vgs Vgoff; 

syms Ml M2 Bl Resultl M3 M4 B2 Result2 signl sign2; 

syms Fcx; 

pi = 3.14159; 

CI = 0.000000000680*1000; 

C2 = 0.000000000075*1000; 

C3 = 0.000000000047*1000; 

LI - 0.0000013*1000; 

Rl - 1000000*1000; 

Rd = 6800*1000; 

f = 1/ (2*pi* (L1*C1*C2/ (CI + C2)) A (l/2)); 

Rll = 2*pi*f*Ll; 

Vdd = 12*1000; 

Vg = sin (2*pi*f ) ; 

Xcl = 1/ (2*pi*f*Cl) ; 

Xc2 = 1/ (2*pi*f*C2) ; 

Xc3 = 1/ (2*pi*f*C3) ; 

Xll = 2*pi*f*Ll; 

IDss = 0.003*1000; 

Vgoff = 1.5*1000; 

GmO = -2*IDss/Vgof f ; 

Vgs = -2*1000; 

Gm = (Vgoff / (Gm0*Vgs) ) + 1; 

Rs = Rd - (1/Gm) ; 

Vss = IDss*Rs; 

Vddl = Gm*Rd*Vg; 

Amplitude = 0.0039; 

Vbias = 0.0039; 

Ml - [(Xll + Xcl + Xc2 + Rll) (-l*Xcl) (-l*Xc2) (0.0);]; 

Ml = [Ml; (-l*Xcl) (Xcl + Rs) (-l*Rs) (0.0);]; 

Ml = [Ml; (-l*Xc2) (-l*Rs) (Xc2 + Rs + Xc3) (0.0);]; 

Ml = [Ml; (0.0) (0.0) (0.0) (Rl + Rd) ; ] ; 

Bl = [(Vg); (Vss); (-l*Vss + Vdds); (Vdd + Vbias);]; 

M2 = [Bid) (-l*Xcl) (-l*Xc2) (0.0);]; 

M2 = [M2; Bl(2) (Xcl + Rs) (-l*Rs) (0.0);]; 

M2 = [M2; Bl(3) (-l*Rs) (Xc2 + Rs + Xc3) (0.0);]; 

M2 = [M2; Bl(4) (0.0) (0.0) (Rl + Rd) ; ] ; 

Resultl = abs (det (M2) /det (Ml) ) ; 

% 

% Now, you want to design a circuit that has Fc times the freq. of 
circuit 1 

% and Ac times the amp. of circuit 1 
syms Fc Ac Acl Xcla Xc2a Rsa; 
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Vg = Acl*sin (2*pi*f*Fc) ; 

% Xc2a is the reactance for the second capacitor adjusted for Fc a 
nd Acl 

Fcx = ( (2*2*pi*pi) / (Fc*Fc*f*f) )* (C1*0. 000000001 + C2*0 . 000000001 ) / 
(C1*0.000000001*C2*0.000000001*L1) ; 
Xcla = 1/ (2*pi*f*Cl*0.000000001*Fcx) ; 
Xc2a = 1/ (2*pi*f*C2*0.000000001*Fcx) ; 
Rsa = Rd*Acl - (1/Gm) ; 

M3 = [(Xll*Acl + Xcla*Acl + Xc2a + Rll*Acl) (-l*Xcla*Acl) (-l*Xc2a 
) (0.0);]; 

M3 = [M3; (-l*Xcla*Acl) (Xcla*Acl + Rsa) (-l*Rsa) (0.0);]; 
M3 = [M3; (-l*Xc2a) (-l*Rsa) (Xc2a + Rsa + Xc3*Acl) (0.0);]; 
M3 = [M3; (0.0) (0.0) (0.0) (Rl*Acl + Rd*Acl);]; 

B2 = [(Vg*Acl); (IDss*Rsa) ; (-1* (IDss*Rsa) + Vdds*Acl) ; (Vdd*Acl + 

Vbias*Acl) ;] ; 
M4 = [B2(l) (-l*Xcla*Acl) (-l*Xc2a) (0.0);]; 
M4 = [M4; B2(2) (Xcla*Acl + Rsa) (-l*Rsa) (0.0);]; 
M4 = [M4; B2(3) (-l*Rsa) (Xc2a + Rsa + Xc3*Acl) (0.0);]; 
M4 = [M4; B2(4) (0.0) (0.0) (Rl*Acl + Rd*Acl);]; 
Result2 - abs (det (M4) /det (M3) ) ; 
signl = abs (det (M2) /det (Ml) ) / (det (M2) /det (Ml) ) ; 
sign2 = abs (det (M4) /det (M3) ) / (det (M4) /det (M3) ) ; 
cw4outl; 

fid = fopen( ! f :/drawl/colpitts4w4 .txt 1 , 'w f ) ; 

syms DesiredAmp DesiredAmpl NewAmp NewAmp2 pAcl Aclx2 Closest Clos 

estAcl NAexp; 

DesiredAmp = 2.564102564; 

Fc = 1000000000/(16983891.6318); 

Closest = 0.0; 

ClosestAcl = 0.0; 

NAexp = abs (eval ( (Result2/sin (2*pi*f *Fc) ) / (Resultl/sin (2*pi*f ) ) ) ) ; 
for Acl = l:1000 f 

NewAmp = abs (eval (NAexp) ) ; 

if abs (NewAmp - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = Acl; 
end 

if abs (NewAmp - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/10; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/10; 
end 
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if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/100; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs(NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/100; 
end 

if abs(NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5 : 5, 

Acl = pAcl + Aclx2/1000; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/1000; 
end 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/10000; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/10000; 
end 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

fprintf (fid, 1 ;\n' ) ; 

fprintf (fid, 1 ; Program colpitts4w4 .m - (c) 2003 Leonard GojerXn 1 ); 
fprintf (fid, f ; Computation of Coil Component Parameters for Drawl\n 
f ) ; 

fprintf (fid, 'An 1 ) ; 

fprintf (fid, f ; Initial Frequency and Amplitude Change To Set Up Dra 
wl\n f ); 

fprintf (fid, 1 ; from Colpitts example in textbook\n 1 ) ; 
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Fn ri nf f 
ipi 111 Li 


f f i H 








-L JJX XII L 1 


f f i H 

11U; 


' • Haqi v*o^ TT 1 v* £!i 1 1 a r\ f* \ i 

/ ucbiicu i!icL|Ucncy 




) ; 


"Fn ri nf f 
J- pX 111 L X 


i 11U^ 


^ jj . loi / rouna \ x ~r c ) ) r 






f print f 




■ \ n M • 






1|JI _L 1 1 LI 


k X X Ll f 


; uyolicU nlupililCaLlOn 




) ; 


X |J -L _Lil LI 


f f -i H 

, X X Ll / 


^ JJ . 1D1 / XJc b X L eLl/ilTip ) f 






f r\r "i n f f 

ipxiriLX 


{■Fir* 
. X XU r 


■ \ n M - 






f n, r t n t f 

1^1 XXI LL 


11U; 


/ lnillal KcSUll x /iinpxxLuae 




) ; 


f nr 1 nff 
lyl ±11 Ll 


X X Ll , 


f %^S 1 fif 1 Rpcn 1 f 1 \ • 
j . loi ^ rxcolii Ll j r 






-L LJ-L 1 i I L -L 










•F ri nf f 
ipi ±11 L X 


, X Id , 


/ oeconaary Kesuit /. Ampxiuucie 




) ; 


F-pv Y* 1 n f F 
i^/i X11L1 


X X Ll f 


ojji xdx f evdi iKesuiL^ ) ; , 






■Prvyi n f "F 

ipi in L X 


X XU/ 


'\n'); 






"Fttn y n n f F~ 

xpx in L X 


i X XQ / 


1 ; Computed Acl\n' ) ; 






F t"*s ri nf f 
XpX X 11 L X 


k x x a f 


1 ; (Component Amplification Factor) 




) ; 


■Fn v n n f "F 
XpX XIX L X 


f f n H 
IlQf 


•%35.16f ' ,Acl) ; 






Xpx XII L X 


f F-i H 

x xa i 


•\n'); 






xpxxn Li 


f -F -i H 


' /Check of Amplification Factor Desired\n'); 




Xpx XIX LI 


' f i r\ 

{ ria, 


1 /(Result 2/Result 1) 




) ; 


■f r\v«l nf f 
Xpx XXX LI 


f f-i H 
i X 1Q f 


'%35.16f 1 ,eval (Result2/Resultl ) ) / 






"Ft> v *i n f F 
X LvX 111 LI 


f f i H 

k X lLl / 


' \n ' ) ; 






fprintf 


(fid, 


' / \n ' ) / 






fprintf 


(fid, 


1 (defun c:drawblks ( ) 1 ) / 






fprintf 


(fid, 


•\n')/ 






fprintf 


(fid, 


' (progn')/ 






fprintf 


(fid, 


' \n ' ) / 







DesiredAmpl = DesiredAmp; 

syms sacm; 

sacm = size(acm); 

for il = l:sacm 

DesiredAmp = acm(il) * DesiredAmpl; 

Fc = 11*1000000000/(16983891.6318); 

Closest = 0.0; 

for Acl - 1:1000, 

NewAmp = abs (eval (NAexp) ) ; 

if abs (NewAmp - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = Acl; 
end 

if abs (NewAmp - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/10; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs(NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 
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ClosestAcl = pAcl + Aclx2/10; 
end 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/100; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/100; 
end 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/1000; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/1000; 
end 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/10000; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/10000; 
end 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 
fprintf (fid, 1 ;\n f ) ; 

f print f ( fid, 1 ; Desired Frequency = f ); 

fprintf (fid, ' %35.16f 1 , round (f*Fc) ) ; 
fprintf (fid, 1 \n' ) ; 

fprintf (fid, 1 ; Desired Amplification = ? ); 
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fprintf (fid, 1 %35 . 16f 1 , DesiredAmp) ; 
fprintf (f id, 1 \n 1 ) ; 

fprintf (fid, f ; Initial Result 1 Amplitude 
fprintf (fid, 1 %35.16f 1 , Result 1) ; 
fprintf ( fid, 1 \n 1 ) ; 

fprintf (fid, 1 ; Secondary Result 2 Amplitude 

fprintf (fid, 1 %35.16f 1 , eval (Result2) ) ; 

fprintf ( fid, 1 \n 1 ) ; 

fprintf (fid, 1 ; Computed Acl\n ! ); 

fprintf (fid, 1 ; (Component Amplification Factor) 

fprintf (fid, ' %35.16f 1 ,Acl) ; 

fprintf ( fid, 1 \n ' ) ; 

fprintf (fid, ' ;Check of Amplification Factor Desired\n 
fprintf (fid, 1 ; (Result 2/Result 1) 
fprintf (fid, 1 %35 . 16f 1 , eval (Result2/Resultl ) ) ; 
fprintf (fid, 1 \n' ) ; 

fprintf (fid, 1 ; Sign 1 (sign of Result 1) 
fprintf (fid, ' %35.16f ' , signl) ; 
fprintf (fid, 1 \n' ) ; 

fprintf (fid, 1 ;Sign 2 (sign of Result 2) 

fprintf (fid, f %35.16f 1 ,eval (sign2) ) ; 

fprintf ( fid, 1 \n 1 ) ; 

fprintf (fid, f ;L1 

fprintf (fid, 1 %35 . 16f f , Ll*Acl ) ; 

fprintf (fid, '\n'); 

fprintf (fid, 1 ;C1 

fprintf (fid, 1 %35 . 16f 1 , eval (Cl*Fcx/Acl) ) ; 
fprintf ( fid, f \n ' ) ; 
fprintf (fid, 1 ;C2 

fprintf (fid, ' %35 . 16f ' , eval (C2*Fcx/Acl) ) ; 

fprintf (fid, 1 \n' ) ; 

fprintf (fid, 1 ;C3 

fprintf (fid, ■ %35 . 16f ' , C3/Acl) ; 

fprintf (fid, 1 \n f ) ; 

fprintf (fid, ' ;Rd 

fprintf (fid, 1 %35 . 1 6f 1 , Rd*Acl ) ; 

fprintf (f id, 1 \n' ) ; 

fprintf (fid, 1 ;R1 

fprintf (fid, 1 %35 . 16f 1 , Rl*Acl ) ; 

fprintf (fid, ' \n 1 ) ; 

fprintf (fid, f ;R11 

fprintf (fid, 1 %35 . 16f 1 , Rll*Acl ) ; 

fprintf (fid, f \n' ) ; 

fprintf ( f id, 1 ; \n 1 ) ; 

fprintf ( fid, 1 ; \n 1 ) ; 

fprintf (fid, 1 (setq count '); 

fprintf (fid, 1 %22 . 16f 1 , il) ; 

fprintf (fid, 1 ) \n f ) ; 
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fprintf (fid, ' (setq DesFreq '); 
fprintf (fid, ' %22 . 16f ', round (f*Fc) ) ; 
fprintf (fid, 1 ) \n' ) ; 
fprintf (fid, ' (setq DesAmp 1 ) ; 
fprintf (fid, ' %22 . 16f 1 , DesiredAmp) ; 
fprintf (fid, 1 ) \n' ) ; 

fprintf (fid, ' (setq InitlAmp '); 

fprintf (fid, ' %22 . 16f 1 , Resultl ) ; 

fprintf (fid, ' ) \n' ) ; 

fprintf (fid, ' (setq Sec2Amp ' ) ; 

fprintf (fid, ' %22 . 16f ' , eval (Result2) ) ; 

fprintf (fid, ' ) \n' ) ; 

fprintf (fid, ' (setq CompAcl '); 

fprintf (fid, • %22 . 16f • , Acl ) ; 

fprintf (fid, • ) \n' ) ; 

fprintf (fid, ' (setq CheckAmp ' ) ; 

fprintf (fid, 1 %22 . 16f ' , eval (Result2/Resultl) ) ; 

fprintf (fid, 1 ) \n' ) ; 

fprintf (fid, 1 (setq Siqnl ' ) ; 

fprintf (fid, '%22.16f ' , signl) ; 

fprintf (fid, ' ) \n' ) ; 

fprintf (fid, ' (setq Sign2 ' ) ; 

fprintf (fid, ' %22 . 16f ' , eval (sign2) ) ; 

fprintf (fid, ' ) \n ' ) ; 

fprintf (fid, 1 (setq LI '); 

fprintf (fid, ' %22 . 16f ' , Ll*Acl) ; 

fprintf (fid, 1 ) \n' ) ; 

fprintf (fid, • (setq CI ' ) ; 

fprintf (fid, ' %22 . 16f ' , eval (Cl*Fcx/Acl) ) ; 

fprintf (fid, 1 ) \n' ) ; 

fprintf (fid, ' (setq C2 '); 

fprintf (fid, ' %22 . 16f • , eval (C2*Fcx/Acl) ) ; 

fprintf (fid, ' )\n' ) ; 

fprintf (fid, • (setq C3 1 ) ; 

fprintf (fid, ' %22 . 16f ' , C3/Acl ) ; 

fprintf (fid, 1 ) \n' ) ; 

fprintf (fid, ' (setq Rd '); 

fprintf (fid, ' %22 . 16f ' , Rd*Acl) ; 

fprintf (fid, • ) \n' ) ; 

fprintf (fid, ' (setq Rl ' ) ; 

fprintf (fid, • %22 . 16f ' , Rl*Acl ) ; 

fprintf (fid, 1 ) \n' ) ; 

fprintf (fid, ' (setqRll • ); 

fprintf (fid, 1 %22 . 16f • , Rll*Acl ) ; 

fprintf (fid, 1 ) \n' ) ; 

fprintf (fid, • (insblkd count DesFreq DesAmp InitlAmp Sec2Amp 
CompAcl CheckAmp Signl Sign2'); 

fprintf (fid, ' LI CI C2 C3 Rd Rl Rll)\n'); 
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fprintf (fid, ' ;\n' ) ; 
end 

fprintf (fid, ' )'); 
fprintf (fid, ' \n' ) ; 
fprintf (fid, • ) ' ) ; 
fprintf (fid, ' \n' ) ; 
f close (fid) ; 
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% Program to Calculate Electrical Component Values for Oscillator 
System of Coil #3 

syms t CI C2 C3 Ll Rl Rs Rd Rll Vbias Vddl Vdd Vdds IDss Vss Vg II 

12 13 14 Xcl Xc2 Xc3 Xll Gm f pi; 
syms Amplitude; 
syms GmO Vgs Vgoff; 

syms Ml M2 Bl Resultl M3 M4 B2 Result2 signl sign2; 

syms Fcx; 

pi = 3.14159; 

CI - 0.000000000680*1000; 

C2 = 0.000000000075*1000; 

C3 = 0.000000000047*1000; 

Ll = 0.0000013*1000; 

Rl - 1000000*1000; 

Rd = 6800*1000; 

f = 1/ (2*pi* (L1*C1*C2/ (CI + C2)) A (l/2)); 

Rll = 2*pi*f*Ll; 

Vdd = 12*1000; 

Vg = sin (2*pi*f ) ; 

Xcl = 1/ (2*pi*f*Cl) ; 

Xc2 = 1/ (2*pi*f*C2) ; 

Xc3 = 1/ (2*pi*f*C3) ; 

Xll = 2*pi*f*Ll; 

IDss = 0.003*1000; 

Vgoff = 1.5*1000; 

GmO = -2*IDss/Vgof f ; 

Vgs = -2*1000; 

Gm = (Vgoff / (Gm0*Vgs) ) + 1; 

Rs = Rd - (1/Gm) ; 

Vss = IDss*Rs; 

Vddl = Gm*Rd*Vg; 

Amplitude = 0.0039; 

Vbias = 0.0039; 

Ml = [(Xll + Xcl + Xc2 + Rll) (-l*Xcl) (-l*Xc2) (0.0);]; 

Ml = [Ml; (-l*Xcl) (Xcl + Rs) (-l*Rs) (0.0);]; 

Ml = [Ml; (-l*Xc2) (-l*Rs) (Xc2 + Rs + Xc3) (0.0);]; 

Ml = [Ml; (0.0) (0.0) (0.0) (Rl + Rd) ; ] ; 

Bl = [ (Vg) ; (Vss); (-l*Vss + Vdds); (Vdd + Vbias);]; 

M2 = [B1(D (-l*Xcl) (-l*Xc2) (0.0);]; 

M2 = [M2; Bl(2) (Xcl + Rs) (-l*Rs) (0.0);]; 

M2 = [M2; Bl(3) (-l*Rs) (Xc2 + Rs + Xc3) (0.0);]; 

M2 = [M2; Bl(4) (0.0) (0.0) (Rl + Rd);]; 

Resultl = abs (det (M2) /det (Ml) ) ; 

% 

% Now, you want to design a circuit that has Fc times the freq. of 
circuit 1 

% and Ac times the amp. of circuit 1 
syms Fc Ac Acl Xcla Xc2a Rsa; 
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Vg = Acl*sin (2*pi*f *Fc) ; 

% Xc2a is the reactance for the second capacitor adjusted for Fc a 
nd Acl 

Fcx = ( (2*2*pi*pi) / (Fc*Fc*f*f) )* (C1*0. 000000001 + C2*0 . 000000001 ) / 
(C1*0.000000001*C2*0.000000001*L1) ; 
Xcla = 1/ (2*pi*f*Cl*0.000000001*Fcx) ; 
Xc2a = 1/ (2*pi*f*C2*0.000000001*Fcx) ; 
Rsa = Rd*Acl - (1/Gm) ; 

M3 = [(Xll*Acl + Xcla*Acl + Xc2a + Rll*Acl) ( -l*Xcla*Acl ) (-l*Xc2a 
) (0.0);]; 

M3 = [M3; (-l*Xcla*Acl) (Xcla*Acl + Rsa) (-l*Rsa) (0.0);]; 
M3 - [M3; (-l*Xc2a) (-l*Rsa) (Xc2a + Rsa + Xc3*Acl) (0.0);]; 
M3 = [M3; (0.0) (0.0) (0.0) (Rl*Acl + Rd*Acl);]; 

B2 = [(Vg*Acl); (IDss*Rsa) ; (-1* ( IDss*Rsa) + Vdds*Acl) ; (Vdd*Acl + 

Vbias*Acl) ; ] ; 
M4 = [B2(l) (-l*Xcla*Acl) (-l*Xc2a) (0.0);]; 
M4 = [M4; B2(2) (Xcla*Acl + Rsa) (-l*Rsa) (0.0);]; 
M4 = [M4; B2(3) (-l*Rsa) (Xc2a + Rsa + Xc3*Acl) (0.0);]; 
M4 = [M4; B2(4) (0.0) (0.0) (Rl*Acl + Rd*Acl);]; 
Result2 = abs (det (M4) /det (M3) ) ; 

signl = abs (det (M2) /det (Ml) ) / (det (M2) /det (Ml) ) ; 
sign2 = abs (det (M4 ) /det (M3) ) / (det (M4 ) /det (M3) ) ; 
cw6outl; 

fid = f open ( 1 f : /drawl/colpitts4w6 . txt 1 , 1 w 1 ) ; 

syms DesiredAmp DesiredAmpl NewAinp NewAmp2 pAcl Aclx2 Closest Clos 

estAcl NAexp; 

DesiredAmp = 2.564102564; 

Fc = 1000000000/(16983891.6318); 

Closest = 0.0; 

ClosestAcl = 0.0; 

NAexp = abs (eval ( (Result2/sin (2*pi*f *Fc) ) / (Resultl/sin (2*pi*f ) ) ) ) ; 
for Acl = 1:1000, 

NewAmp = abs (eval (NAexp) ) ; 

if abs (NewAmp - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = Acl; 
end 

if abs (NewAmp - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5/ 

Acl = pAcl + Aclx2/10; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/10; 
end 
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if abs(NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/100; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs(NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/100; 
end 

if abs(NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/1000; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs(NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/1000; 
end 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/10000; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/10000; 
end 

if abs(NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

f print f (fid, ' ;\n f ) ; 

fprintf (fid, 1 ; Program colpitts4w6 .m - (c) 2003 Leonard GojerXn 1 ); 
fprintf (fid, 1 ; Computation of Coil Component Parameters for Drawl\n 
1 ) ; 

fprintf (fid, 1 ;\n' ) ; 

fprintf (fid, 1 ; Initial Frequency and Amplitude Change To Set Up Dra 
wl\n'); 

fprintf (fid, 1 ; from Colpitts example in textbook\n 1 ) ; 
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f print f 


(f id, 


' \n 1 ) ; 






f tor int f 

-i- £S X Xll l-» X 


|f id, 


1 .* Dp^irpd Frpfiiipnrv 




1 \ • 
) i 


f nrint f 

1— V_/ 1_ _1_ 11 L> 1. 


If id, 


1 %35 16f 1 . round ( f *Fr ) ) • 

O J J • 1 U 1 f X W ull kX ^ X E \^ J f f 






fnr "i riff 

X X -LAX L. X 


f fid 








f orint f 

K-' -I* «X» A X w -X» 


(fid, 


1 ; Desired Amolif icat i on 




) r 


f orint f 


If id, 


1 %35 . 16f ' . DesiredAmn) : 






f nrint f 

X. fy/ -I— _L 11 


|f id, 


1 \n 1 ) ; 






f or int f 

-t- -1- -I- i- X w J— 


|f id, 


1 ; Initial Rpsul t 1 Amr>l i fnHp 

f XllXk^XUX ij IX X X ill 11 ky X X L~ IX vX v3- 




) / 


f orint f 


If id, 


1 %35 . 16f 1 . Resultl ) ; 






f print f 


If id, 


1 \n 1 ) ; 






fnr -i n t" "F 

J- ky JL JL i 1 L. X 


f f id 


z oct-uiiuax y incouil. ^ /"UiLkyxX L. Livj.tr 




) r 


f print f 


( fid, 


'%35 16f ' pval (Rpsnlr?) } • 






ftjr i nff 

x. ky x. xii x 


f f id 

J — L kX , 








"Ft> t H n t* "F 

j_ ky J Lii u x 


f f i d 

-L X Li. f 


f ^ wlu|y Lx L- C. LJ. r\L> x. \il y / 






fnri nr f 

J- k> X. Xll 1 X. 


f fid 


/ \ l^ ^iLi^y wi id i l, riiUL^x xixua i — l^jii x d ^ u w j_ y 




) / 


"Fdt "i nt f 

x jy x x 1 1 l. x. 


f fid 

, J — L kX f 








f print f 


(fid, 


1 \n 1 ) ; 






"FrNT i n t" "P 

x. ky i xll ui 


f fid 

J L\JL f 


/V^iicv^jv ui nllipill XL-a LlvJll £ au LUI jJcolicU \ il 






for int f 


(fid, 


1 ; (Result ?/Rpsulr 1^ 




) r 


f print f 


If id, 


' %35 16f 1 pval (Re^ul t? /RpshI r1 ) ) • 






for int f 

J_ KS -I— _L. 11 U X. 


! fid. 

, X X vx , 


1 \n 1 ) ; 






fprintf 


(fid, 


' ; \n ' ) ; 






fprintf 


(fid, 


1 (def un c : drawblks ( ) ' ) ; 






fprintf 


(fid, 


•\n'); 






fprintf 


(fid, 


' (progn 1 ); 






fprintf 


(fid, 


•\n'); 







DesiredAmpl = DesiredAmp; 

syms sacm; 

sacm = size(acm); 

for il = lrsacm 

DesiredAmp = acm(il) * DesiredAmpl; 

Fc = 11*1000000000/(16983891.6318); 

Closest = 0.0; 

for Acl = 1:1000, 

NewAmp = abs (eval (NAexp) ) ; 

if abs (NewAmp - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = Acl; 
end 

if abs (NewAmp - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/10; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 
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ClosestAcl = pAcl + Aclx2/10; 
end 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/100; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs(NewAmp2 - DesiredAmp) 



ClosestAcl = pAcl + Aclx2/100; 
end 

if abs (NewAmp2 - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 
pAcl = Acl; 
for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/1000; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs(NewAmp2 - DesiredAmp) 



< abs (Closest 

< abs (Closest 



DesiredAmp) 
DesiredAmp) 



ClosestAcl = pAcl + Aclx2/1000; 
end 

if abs (NewAmp2 - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/10000; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) < abs 
ClosestAcl = pAcl + Aclx2/10000; 

end 

if abs(NewAmp2 - DesiredAmp) < abs 

Closest = NewAmp2; 
end 
end 

Acl - ClosestAcl; 

fprintf (fid, ' ;\n f ) ; 

f printf (fid, 1 ; Desired Frequency 

fprintf (fid, 1 %35 . 16f ', round ( f*Fc) ) ; 

fprintf (fid, 1 \n' ) ; 

fprintf (fid, 1 ; Desired Amplification 



< abs (Closest 

< abs (Closest 



DesiredAmp) 
DesiredAmp) 



(Closest 
(Closest 



DesiredAmp) 
DesiredAmp) 



= '); 



= '); 
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fprintf (fid, ' %35 . 16f ' , DesiredAmp) ; 
fprintf (fid, ' \n' ) ; 

fprintf (fid, ' /Initial Result 1 Amplitude 
fprintf (fid, 1 %35 .-16f ' , Resultl) ; 
fprintf (fid, ' \n' ) ; 

fprintf (fid, ' /Secondary Result 2 Amplitude 

fprintf (fid, ' %35 . 16f ' , eval (Result2) ); 

fprintf (fid, 1 \n' ) / 

fprintf (fid, ' /Computed Acl\n' ) / 

fprintf (fid, '/ (Component Amplification Factor) 

fprintf (fid, ' %35 . 16f ' , Acl) / 

fprintf (fid, ' \n' ) / 

fprintf (fid, ' /Check of Amplification Factor Desired\n 
fprintf (fid, '/ (Result 2/Result 1) 
fprintf (fid, 1 %35 . 16f ' , eval (Result2/Resultl ) ) / 
fprintf (fid, » \n' ) ; 

fprintf (fid, ' /Sign 1 (sign of Result 1) 
fprintf (fid, ' %35. 16f ' , signl) ; 
fprintf (fid, • \n' ) ; 

fprintf (fid, 1 /Sign 2 (sign of Result 2) 

fprintf (fid, ' %35 . 16f ' , eval (sign2) ) / 

fprintf (fid, '\n' ) / 

fprintf (fid, ' /LI 

fprintf (fid, ' %35 . 16f ' , Ll*Acl) / 

fprintf (fid, ' \n u y ; 

fprintf (fid, ' /CI 

fprintf (fid, ' %35 . 16f ' , eval (Cl*Fcx/Acl) ) / 
fprintf (fid, ' \n' ) / 
fprintf (fid, ' /C2 

fprintf (fid, ' %35 . 16f ■ , eval (C2*Fcx/Acl ) ) / 

fprintf (fid, ' \n' ) / 

fprintf (fid, 1 /C3 

fprintf (fid, ' %35 . 16f ' , C3/Acl ) / 

fprintf (fid, ' \n' ) / 

fprintf (fid, ' /Rd 

fprintf (fid, ' %35 . 16f * , Rd*Acl ) / 

fprintf (fid, ' \n 1 ) / 

fprintf (fid, ' /Rl 

fprintf (fid, • %35 . 16f ' , Rl*Acl) / 

fprintf (fid, ' \n' ) / 

fprintf (fid, ' ;R11 

fprintf (fid, ' %35 . 16f ' , Rll*Acl ) / 

fprintf (fid, 1 \n' ) / 

fprintf (fid, ' / \n ' ) / 

fprintf ( fid, ' / \n ' ) / 

fprintf (fid, ' (setq count ')/ 

fprintf (fid, •%22.16f ,il) / 

fprintf (fid, » )\n' ) / 
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fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fpr-in-tf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 
fprintf 



(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 
(fid 



%22 
)\n 

%22 
)\n 

%22 
)\n 

%22 
)\n 

%22 
)\n 

%22 
)\n 

%22 
)\n 

%22 
)\n 

%22 
)\n 

%22 
)\n 

%22 
)\n 

%22 
)\n 

%22 
)\n 

%22 
)\n 

%22 
)\n 



CompAcl CheckAmp S 
fprintf (fid, ' LI 



(setq DesFreq '); 
.16f round (f*Fc) ) ; 
'); 

(setq DesAmp ' ) ; 
. 16f ' , DesiredAmp) ; 
') ; 

(setq InitlAmp '); 
.16f 1 ,Resultl) ; 
') ; 

(setq Sec2Amp ' ) ; 
.16f • ,eval (Result2) ) ; 
'); 

(setq CompAcl 1 ) ; 
.16f ' ,Acl) ; 
') ; 

(setq CheckAmp '); 
.16f ' ,eval (Result2/Resultl) ) ; 
' ) ; 

(setq Signl 1 ) ; 
. 16f 1 , signl ) ; 
') ; 

(setq Sign2 ' ) ; 
. 16f ' , eval (sign2) ) ; 
*) ; 

(setq LI » ) ; 
.16f ' ,Ll*Acl) ; 

■) ; 

-(setq CI 1 ) ; 

.16f ' f eval (Cl*Fcx/Acl) ) ; 

•) ; 

~(setq C2 ' ) ; 

. 16f ' ,eval (C2*Fcx/Acl) ) ; 

') ; 

(setq C3 ' ) ; 
.16f ' ,C3/Acl) ; 
') ; 

(setq Rd ' ) ; 
.16f ' ,Rd*Acl) ; 
'); 

(setq Rl ' ) ; 
.16f * ,Rl*Acl) ; 
' ) ; 

(setq Rll ' ) ; 
.16f • ,Rll*Acl) ; 
' ) ; 

(insblkd count DesFreq DesAmp InitlAmp Sec2Amp 
ignl Sign2 ' ) ; 
CI C2 C3 Rd Rl Rll) \n' ) ; 
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fprintf (fid, ' ;\n' ) ; 
end 

f print f (fid, ' )'); 
fprintf (fid, • \n' ) ; 
f print f (fid, • ) ' ) ; 
fprintf (fid, '\n' ) ; 
f close (fid) ; 
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% Program to Calculate Electrical Component Values for Oscillator 
System of Coil #4 

syms t CI C2 C3 LI Rl Rs Rd Rll Vbias Vddl Vdd Vdds IDss Vss Vg II 

12 13 14 Xcl Xc2 Xc3 Xll Gm f pi; 
syms Amplitude; 
syms GmO Vgs Vgoff; 

syms Ml M2 Bl Resultl M3 M4 B2 Result2 signl sign2; 

syms Fcx; 

pi = 3.14159; 

CI = 0.000000000680*1000; 

C2 = 0.000000000075*1000; 

C3 = 0.000000000047*1000; 

LI = 0.0000013*1000; 

Rl = 1000000*1000; 

Rd = 6800*1000; 

f = 1/ (2*pi* (L1*C1*C2/ (CI + C2)) A (l/2)); 

Rll = 2*pi*f*Ll; 

Vdd = 12*1000; 

Vg = sin(2*pi*f ) ; 

Xcl = 1/ (2*pi*f*Cl) ; 

Xc2 = 1/ (2*pi*f*C2) ; 

Xc3 = 1/ (2*pi*f*C3) ; 

Xll = 2*pi*f*Ll; 

IDss = 0.003*1000; 

Vgoff = 1.5*1000; 

GmO = -2*IDss/Vgof f ; 

Vgs = -2*1000; 

Gm = (Vgof f / (Gm0*Vgs) ) + 1; 

Rs = Rd - (1/Gm) ; 

Vss = IDss*Rs; 

Vddl = Gm*Rd*Vg; 

Amplitude = 0.0039; 

Vbias = 0.0039; 

Ml = [(Xll + Xcl + Xc2 + Rll) (-l*Xcl) (-l*Xc2) (0.0);]; 

Ml = [Ml; (-l*Xcl) (Xcl + Rs) (-l*Rs) (0.0);]; 

Ml = [Ml; (-l*Xc2) (-l*Rs) (Xc2 + Rs + Xc3) (0.0);]; 

Ml = [Ml; (0.0) (0.0) (0.0) (Rl + Rd) ; ] ; 

Bl = [ (Vg) ; (Vss); (-l*Vss + Vdds); (Vdd + Vbias);]; 

M2 = [B1(D (-l*Xcl) (-l*Xc2) (0.0);]; 

M2 = [M2; Bl(2) (Xcl + Rs) (-l*Rs) (0.0);]; 

M2 = [M2; Bl(3) (-l*Rs) (Xc2 + Rs + Xc3) (0.0);]; 

M2 = [M2; Bl(4) (0.0) (0.0) (Rl + Rd);]; 

Resultl = abs (det (M2) /det (Ml) ) ; 

% 

% Now, you want to design a circuit that has Fc times the freq. of 
circuit 1 

% and Ac times the amp. of circuit 1 
syms Fc Ac Acl Xcla Xc2a Rsa; 
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Vg = Acl*sin(2*pi*f*Fc) ; 

% Xc2a is the reactance for the second capacitor adjusted for Fc a 
nd Acl 

Fcx = ( (2*2*pi*pi) / (Fc*Fc*f*f) )* (C1*0. 000000001 + C2*0 . 000000001 ) / 
(C1*0.000000001*C2*0.000000001*L1) ; 
Xcla = 1/ (2*pi*f*Cl*0.000000001*Fcx) ; 
Xc2a = l/(2*pi*f*C2*0.000000001*Fcx); 
Rsa = Rd*Acl - (1/Gm) ; 

M3 = [(Xll*Acl + Xcla*Acl + Xc2a + Rll*Acl) ( -l*Xcla*Acl ) (-l*Xc2a 
) (0.0);]; 

M3 = [M3; (-l*Xcla*Acl) (Xcla*Acl + Rsa) (-l*Rsa) (0.0);]; 
M3 = [M3; (-l*Xc2a) (-l*Rsa) (Xc2a + Rsa + Xc3*Acl) (0.0);]; 
M3 = [M3; (0.0) (0.0) (0.0) (Rl*Acl + Rd*Acl);]; 

B2 = [(Vg*Acl); (IDss*Rsa) ; (-1* (IDss*Rsa) + Vdds*Acl) ; (Vdd*Acl + 

Vbias*Acl) ; ] ; 
M4 = [B2(l) (-l*Xcla*Acl) (-l*Xc2a) (0.0);]; 
M4 = [M4; B2(2) (Xcla*Acl + Rsa) (-l*Rsa) (0.0);]; 
M4 = [M4; B2(3) (-l*Rsa) (Xc2a + Rsa + Xc3*Acl) (0.0);]; 
M4 = [M4; B2(4) (0.0) (0.0) (Rl*Acl + Rd*Acl);]; 
Result2 = abs (det (M4) /det (M3) ) ; 

signl = abs (det (M2) /det (Ml) ) / (det (M2) /det (Ml) ) ; 
sign2 - abs (det (M4 ) /det (M3) ) / (det (M4 ) /det (M3) ) ; 
cw8outl; 

fid = fopen ( ' f : /drawl/colpitts4w8 . txt 1 / 1 w 1 ) ; 

syms DesiredAmp DesiredAmpl NewAmp NewAmp2 pAcl Aclx2 Closest Clos 

estAcl NAexp; 

DesiredAmp = 2.564102564; 

Fc = 1000000000/(16983891.6318); 

Closest = 0.0; 

ClosestAcl = 0.0; 

NAexp = abs (eval ( (Result2/sin (2*pi*f *Fc) ) / (Resultl/sin (2*pi*f ) ) ) ) ; 
for Acl - 1:1000, 

NewAmp = abs (eval (NAexp) ) ; 

if abs (NewAmp - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = Acl; 
end 

if abs (NewAmp - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/10; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp 2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/10; 
end 
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if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/100; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/100; 
end 

if abs(NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/1000; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs(NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/1000; 
end 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/10000; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs(NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = pAcl + Aclx2/10000; 
end 

if abs (NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

fprintf (fid, ? ;\n' ) ; 

fprintf (fid, 1 ; Program colpitts4w8 .m - (c) 2003 Leonard GojerXn 1 ); 
fprintf (fid, 1 ; Computation of Coil Component Parameters for Drawl\n 
1 ) ; 

fprintf (fid, 1 ;\n f ) ; 

fprintf (fid, 1 ; Initial Frequency and Amplitude Change To Set Up Dra 
wl\n' ) ; 

fprintf (fid, 1 ; from Colpitts example in textbook\n 1 ) ; 



Page 3 



colpitts4w8 .m 



■p ir\ >»* A y~\ 4— 4^ 

ipririL x 


i X 1U , 








J_ y>L XII U X 


k X 1U. , 


' • O "J >~ <J TP >~ /~» ■« f /-\ ym. ~ , * 

f uesirea r requency 




) ; 


p v~\ k i n 1" f 

iprinti 


i X XU/ 


1 §r "3 R 1^471 mnnrl / f * \ \ . 

^jO.ioi / rounu (i T i:c) ) ; 






iprinur 


f -Ft H 


\n ) , 






p t*n y "i ri I" f 
J_ px 




, uesirea MipiiriCaLion 




) ; 


"Ft^ v i n 1" ~F 
ijjiinLi 


k X Xvu / 


^ jo. 101 / uesxrea/imp ) , 






fprintf 


f -F -J H 


\n ) ; 






iprinur 


f -F -1 H 
i -LIU/ 


, iniuidi Kesuit 1 Amplitude 




) / 


ipi III LI 


, X lU , 


tijj. lor , KesuiLi ) , 






f r\ y! riff 

ipx XII U X 


( -Fi H 


\n ) , 






ipi inn 




f beconaary Kesuit z Amplitude 




) ; 


ipi XII L X 


f f -i H 
,11(1; 


"5 0 o • 1 ox r evai [ KeSUltZ ; ) f 






P y-\ y -1 n +• f 

iprinLi 


r -Fn H 

, x la, 


' \n • ) ; 






iprinur 


r -f -i h 
, x la , 


' ; Computed Acl\n' ) ; 






p t~\ v* n 4~ "F 

iprinui 


f -F -1 H 

, iia, 


' ; (Component Amplification Factor) 




) ; 


4~ *i *-» 4— 4- 

iprinir 


t ria, 


•%35.16f ' ,Acl) ; 






rpriiiL r 


f -F A H 

k iia, 


•\n'); 






f print f 


, iia, 


'; Check of Amplification Factor Desired\n'); 




x p X X 11 U X 


x ia , 


1 ; (Result 2/Result 1) 




) ; 


xpx XII X, X 


k x xa , 


1 %35. 16f ' , eval (Result2/Resultl) ) ; 






Xpx HIT- X 


f -Ft H 


' \n') ; 






f prin-tf 


(fid, 


■ ; \n • ) 7 






fprintf 


(fid, 


' (def un c : drawblks ( ) ' ) ; 






fprintf 


(fid, 


•\n'); 






fprintf 


(fid, 


' (progn'); 






fprintf 


(fid, 


' \n ' ) ; 







DesiredAmpl = DesiredAmp; 
syms sacm; ^ 
sacm = size(acm); 
for il = l:sacm 

DesiredAmp = acm(il) * DesiredAmpl; 

Fc = il*1000000000/ (16983891.6318) ; 

Closest = 0.0; 

for Acl = 1:1000, 

NewAmp = abs (eval (NAexp) ) ; 

if abs (NewAmp - DesiredAmp) < abs (Closest - DesiredAmp) 

ClosestAcl = Acl; 
end 

if abs (NewAmp - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/10; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs(NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 
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ClosestAcl = pAcl + Aclx2/10; 
end 

if abs(NewAmp2 - DesiredAmp) < abs (Closest - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 

pAcl = Acl; 

for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/100; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs(NewAmp2 - DesiredAmp) < abs 
ClosestAcl = pAcl + Aclx2/100; 

end 

if abs (NewAmp2 - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 
pAcl = Acl; 
for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/1000; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs(NewAmp2 - DesiredAmp) 



(Closest 
< abs (Closest 



DesiredAmp) 
DesiredAmp) 



ClosestAcl = pAcl + Aclx2/1000; 
end 

if abs (NewAmp2 - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 
pAcl = Acl; 
for Aclx2 = -5:5, 

Acl = pAcl + Aclx2/10000; 

NewAmp2 = abs (eval (NAexp) ) ; 

if abs (NewAmp2 - DesiredAmp) 



< abs (Closest 

< abs (Closest 



DesiredAmp) 
DesiredAmp) 



ClosestAcl = pAcl + Aclx2/10000; 
end 

if abs (NewAmp2 - DesiredAmp) 

Closest = NewAmp2; 
end 
end 

Acl = ClosestAcl; 



< abs (Closest 

< abs (Closest 



DesiredAmp) 
DesiredAmp) 



fprintf (fid, 
fprintf (fid, 
fprintf ( fid, 
fprintf (fid, 
fprintf (fid, 



;\n' ) ; 

; Desired Frequency 
%35.16f 1 ,round(f*Fc) ) ; 
\n') ; 

; Desired Amplification 



= f ); 
= Mr- 
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fprintf (fid, 1 %35 » 16f 1 , DesiredAmp) ; 
fprintf (fid, ? \n f ) ; 

fprintf (fid, 1 ; Initial Result 1 Amplitude 
fprintf (fid, 1 %35 . 16f 1 , Result 1) ; 
fprintf (fid, 1 \n f ) ; 

fprintf (fid, 1 /Secondary Result 2 Amplitude 

fprintf (fid, 1 %35 . 16f 1 , eval (Result2) ) ; 

fprintf (fid, 1 \n' ) ; 

fprintf (fid, 1 ; Computed Acl\n f ); 

fprintf (fid, 1 ; (Component Amplification Factor) 

fprintf (fid, ' %35.16f ■ ,Acl) ; 

fprintf (fid, 1 \n' ) ; 

fprintf (fid, 1 ;Check of Amplification Factor Desired\ 
fprintf (fid, ' ; (Result 2/Result 1) 
fprintf (fid, '%35.16f 1 , eval (Result2/Resultl ) ) ; 
fprintf (f id, 1 \n 1 ) ; 

fprintf (fid, 1 ; Sign 1 (sign of Result 1) 
fprintf (fid, ' %35.16f 1 ,signl) ; 
fprintf (fid, 1 \n 1 ) ; 

fprintf (fid, 1 ; Sign 2 (sign of Result 2) 

fprintf (fid, ? %35.16f 1 f eval (sign2) ) ; 

fprintf ( fid, 1 \n 1 Y-; 

fprintf (fid, ';L1 

fprintf (fid, '%35.16f 1 ,Ll*Acl) ; 

fprintf (fid, f \n^Yr 

fprintf (fid, 1 ;C1 

fprintf (fid, 1 %35 . 16f 1 , eval (Cl*Fcx/Acl) ) ; 
fprintf ( fid, ? \n 1 } ; 
fprintf (fid, 1 ;C2 

fprintf (fid, 1 %35 . 16f f , eval (C2*Fcx/Acl) ) ; 

fprintf (fid, 1 \n f Y; 

fprintf (fid, 1 ;C3 

fprintf (fid, 1 %35 . 16f ' , C3/Acl ) ; 

fprintf (f id, 1 \n f ) ; 

fprintf ( fid, 1 ; Rd 

fprintf (fid, 1 %35 . 16f 1 , Rd*Acl ) ; 

fprintf (fid, 1 \n f ) ; 

fprintf (fid, ' ;R1 

fprintf (fid, 1 %35 . 16f f , Rl*Acl ) ; 

fprintf (f id, 1 \n ! ) ; 

fprintf (fid, ' ;R11 

fprintf (fid, 1 %35 . 16f 1 , Rll*Acl ) ; 

fprintf (fid, 1 \n' ) ; 

fprintf (f id, 1 ; \n' ) ; 

fprintf (fid, 1 ; \n 1 ) ; 

fprintf ( fid, 1 (setq count 

fprintf (fid, ! %22.16f 1 ,il) ; 

fprintf (fid, 1 ) \n' ) ; 
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fprintf (f id, ' (setq DesFreq '); 
fprintf (fid, ' %22 . 16f ', round (f*Fc) ) ; 
fprintf (fid, 1 ) \n' ) ; 
fprintf (fid, ' (setq DesAmp '); 
fprintf (fid, ' %22 . 16f • , DesiredAmp) ; 
fprintf (fid, ' ) \n' ) ; 

fprintf (fid, ' (setq InitlAmp ' ) ; 

fprintf (fid, 1 %22 . 16f ', Result 1 ) ; 

fprintf (fid, ' ) \n' ) ; 

fprintf (fid, ' (setq Sec2Amp ' ) ; 

fprintf (fid, ' %22 . 16f ' , eval (Result2 ) ) ; 

fprintf (fid, ' ) \n' ) ; 

fprintf (fid, ' (setq CompAcl '); 

fprintf (fid, • %22 . 16f • , Acl) ; 

fprintf (fid, • ) \n' ) ; 

fprintf ( fid, * (setq CheckAmp ' ) ; 

fprintf (fid, 1 %22 . 16f 1 , eval (Result2/Resultl) ) ; 

fprintf (fid, ' ) \n ! ) ; 

fprintf (fid, ' (setq Signl '); 

fprintf (fid, 1 %22 . 16f ' , signl) ; 

fprintf (fid, ' ) \n' ) ; 

fprintf (fid, ' (setq Sign2 » ) ; 

fprintf (fid, '%22.16f ' ,eval(sign2) ) ; 

fprintf (fid, ') \n' ) ; 

fprintf (fid, ' (setq LI '); 

fprintf (fid, ' %22 . 16f • , Ll*Acl) ; 

fprintf (fid, • )\n' ) ; 

fprintf (fid, • -(setq CI '); 

fprintf (fid, ' %22 . 16f ' , eval (Cl*Fcx/Acl ) ) ; 

fprintf (fid, ' ) \n' ) ; 

fprintf (fid, • (setq C2 ' ) ; 

fprintf (fid, • %22 . 16f * , eval (C2*Fcx/Acl ) ) ; 

fprintf (fid, 1 ) \n' ) ; 

fprintf (fid, ' (setq C3 ' ); 

fprintf (fid, '%22.16f ' ,C3/Acl) ; 

fprintf (fid, » ) \n' ) ; 

fprintf (fid, ' (setq Rd 1 ) ; 

fprintf (fid, ' %22 . 16f • , Rd*Acl) ; 

fprintf (fid, • ) \n' ) ; 

fprintf (fid, ' (setq Rl ' ) ; 

fprintf (fid, ' %22 . 16f ' , Rl*Acl) ; 

fprintf (fid, 1 ) \n' ) ; 

fprintf (fid, 1 (setq Rll • ); 

fprintf (fid, ' %22 . 16f 1 , Rll*Acl ) ; 

fprintf (fid, 1 ) \n' ) ; 

fprintf (fid, ' (insblkd count DesFreq DesAmp InitlAmp Sec2Amp 

CompAcl CheckAmp Signl Sign2 ' ) ; 

fprintf (fid, 1 LI CI C2 C3 Rd Rl Rll)\n'); 
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fprintf (fid, • ; \n'J ; 
end 

fprintf (fid, ' ) » ) ; 
fprintf (fid, ' \n' ) ; 
fprintf -(f id, 1 ) ' ) ; 
fprintf (fid, ' \n' ) ; 
f close (fid) ; 
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%coilwave4.m - Program to Compute Amplitude Values for Oscillators of Coil #4 of 
Magnetic Bottle 

% " - waveform is Specified by interpolation of 200 term fourier series in 

Matlab 

% 



syms ax ay bx by cx cy dx dy ex ey fx fy; 
syms gx gy hx hy ix iy jx jy kx ky lx ly; 

syms mx my nx ny ox oy xof xofa yof yof2 ml m2 dxl dxldyl xvar yvar; 

syms al bl cl fl a2 b2 c2 a3 b3 c3 sumll suml2 terml term2; 

syms a4 b4 c4 ala bla cla fla alb bib clb exprl expr2 pi ; 

syms ale blc clc fib flc f2a f2b alq blq clq fq asubn bsubn; 

syms L k f3a f3b f4a f4b f4c x il 12 kl loops iiil foura fourb; 

pi = 3.141592654; 

L = 10; 

k = 2; 

kl = 6.25; 

ml = 2.5; 

m2 = -0.5; 

ax = 0; 

bx s= k - kl + ax; 

cx = k + kl + ax; 

dx = L - k - kl + ax; 

ex = L - k + kl + ax; 

fx = L + k - kl + ax; 

gx = L + k + kl + ax; 

hx = 3*L - gx; 

ix = 3*L - fx; 

jx = 2*L; 

kx = k + ax; 

lx = 2*L - k + ax; 

mx = L - k + ax; 

nx = L + k + ax; 

ox = 2*L - k; 

ay = 0; 

by = ml*(k - kl) + ay; 

cy = m2*(kl) + ml*k + ay; 

dy = m2*(L - 2*k - kl) + ml*k + ay; 

ey = m2*(L - k + kl) + ml*k + ay; 

fy = m2*(L + k - kl) + ml*k + ay; 

gy = m2*(k + kl) + ay; 

hy = gy; 

iy = fy; 

jy = 0; 

ky = ml*k + ay; 

ly = -ml*k + ay; 

my = m2*(L - 2*k) + ml*k + ay; 

ny = -l*(dy+ey)/2; 

oy = -my; 

fid = f open ( , g:/drawl/coilwave4. out 1 , 'w') ; 
fprintf (fid, 1 Fourier Series of Coil #4 WaveForm 1 ); 
fprintf(fid, An 1 ); 

fprintf (fid, 1 '); 

fprintf (fid, 'XrT) ; 
fprintf (fid, 'VT); 
dxl = 2*L/1000; 
for iil = 0:100 

sumll =0.0; 

sum!2 = 0.0; 

for iiil = 1:500 
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xofa = i-iil*2*L/1000; 

xof = (L/2) + xofa; 
% _ 

%fl 

%2al*bx + bl = ml 
%2al*cx + bl = m2 
%cl = -al*bxA2 - bl*bx 

% 

al = (ml - m2)/(2*bx - 2*cx); 

bl = (2*bx*m2 - 2*cx*ml)/(2*bx - 2*cx) ; 

cl = (-l*al)*bx*bx - bl*bx; 

fl = al*xof*xof + bl*xof + cl; 



%fla 

%is fl flipped 

% 

fla = al*(L - xof)*(L - xof) + bl*(L - xof) + cl; 

% 

%f2a 
%dy+fla 

% 

f2a = fla + dy; 

% 

%f2b 
%my+fla 

% 

f2b = fla + my; 

% 

%flb 

%-fl of (xof + L) 



alb = (-ml + m2)/(2*fx - 2*gx) ; 

bib = (2*fx*(-m2) - 2*gx*(-ml))/(2*fx - 2*gx) ; 

clb = (-l*alb)*fx*fx - blb*fx; 

fib = alb*xof*xof + blb*xof + clb; 

ale = (ml - m2)/(2*hx - 2*ix); 

blc = (2*hx*(m2) - 2*i x*(ml))/(2*hx - 2*ix); 
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clc = (-l*alc)*hx*hx - blc*hx; 

% 

%f3a 
%fy-flb 

% 

f3a = f y - fib; 

% 

%f3b 
%ny-flb 

% 

f3b = ny - fib; 

% 

%f4c 

%-fla at 2nd point 

% 

f4c = -alc*(xof)*(xof) - blc*(xof) - clc; 

% 

%f4a 

%-cy + f4c 

% 

f4a = -f4c - cy; 

% 

%f4b 

%f4c -cy 

% 

f4b = -f4c - cy; 
% _ 

%returns 

%x,y for x, cycle # 

% 

if (xof < bx) 

yof = ay + ml*xof ; 
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end 

if (xof >= bx) 

if (xof < cx) 

yof = by + al*xof*xof + bl*xof + cl; 

end 
end 

if (xof >= cx) 

if (xof < dx) 

yof = cy + m2*(xof - cx) ; 

end 
end 

if (xof >= dx) 

if (xof < (L - k)) 
yof = f2a; 

end 
end 

if (xof >= (L - k)) 

if (xof < ex) 
yof = f2b; 

end 
end 

if (xof >= ex) 

if (xof < fx) 

yof = ey + m2*(xof - ex); 

end 
end 

if (xof >= fx) 

if (xof < (L + k)) 
yof = f3a; 

end 
end 

if (xof >= (L + k)) 

if (xof < gx) 
yof = f3b; 

end 
end 

if (xof >= gx) 

if (xof < hx) 
yof = gy + m2*(xof - gx) ; 

end 
end 

if (xof >= hx) 

if (xof < ox) 
yof = f4a; 

end 
end 

if (xof >= ox) 

if (xof < ix) 
yof = f4b; 

end 
end 

if (xof >= ix) 

if (xof < jx) 

yof = -by + ml* (xof - ix); 

end 
end 

if (xof == jx) 

yof =0.0; 
end 

terml = cos(pi*iil*xof/L) ; 
term2 = sin(pi*iil*xof/l_) ; 
yof 2 = -l*yof; 
dxldyl = dxl*yof2; 
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exprl = terml*dxldyl; 
expr2 = term2*dxldyl; 
sumll = sumll + exprl; 
suml2 = suml2 + expr2; 
end 

for iiil = 501:1000 

xofa = iiil*2*L/1000; 

xof = xofa + (L/2) ; 
% 

%fl 

%2al*bx + bl = ml 
%2al*cx + bl = m2 
%cl = -al*bxA2 - bl*bx 

% 

al = (ml - m2)/(2*bx - 2*cx) ; 

bl = (2*bx*m2 - 2*cx*ml)/(2*bx - 2*cx) ; 

cl = (-l*al)*bx*bx - bl*bx; 

fl = al*xof*xof + bl*xof + cl; 

% 

%fla 

%is fl flipped 

% 

fla = al*(L - xof)*(L - xof) + bl*(L - xof) + cl; 

% 

%f2a 
%dy+fla 

% 

f2a = fla + dy; 

% 

%f2b 
%my+fla 

% 

f2b = fla + my; 

% 

%flb 

%-fl of (xof + L) 

% 
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alb = (-ml + m2)/(2*fx - 2*gx) ; 

bib = (2*fx*(-m2) - 2*gx*(-ml))/(2*fx - 2*gx) ; 

clb = (-l*alb)*fx*fx - blb*fx; 

fib = alb*xof*xof + blb*xof + clb; 

ale = (ml - m2)/(2*hx - 2*ix); 

blc = (2*hx*(m2) - 2*i x*(ml))/(2*hx - 2*ix); 

clc = (-l*alc)*hx*hx - blc*hx; 

% 

%f3a 
%fy-flb 

% 

f3a = f y - fib; 

% 

%f3b 
%ny-flb 

% 

f3b = ny - fib; 

% 

%f4c 

%-fla at 2nd point 

% 

f4c = -alc*(xof)*(xof) - blc*(xof) - clc; 

% 

Xf4a 

%-cy + f4c 

% 

f4a = -f4c - cy; 

% 

%f4b 

%f4c -cy 

% 

f4b = -f4c - cy; 

% 

%returns 
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%x,y for x, cycle # 

% 



if (xof < bx) 

yof = ay + ml*xof ; 
end 

if (xof >= bx) 

if (xof < cx) 

yof = by + al*xof*xof + bl*xof + cl; 

end 
end 

if (xof >= cx) 

if (xof < dx) 

yof = cy + m2*(xof - cx) ; 

end 
end 

if (xof >= dx) 

if (xof < (L - k)) 
yof = f2a; 

end 
end 

if (xof >= (L - k)) 

if (xof < ex) 
yof = f2b; 

end 
end 

if (xof >= ex) 

if (xof < fx) 

yof = ey + m2*(xof - ex); 

end 
end 

if (xof >= fx) 

if (xof < (L + k)) 
yof = f3a; 

end 
end 

if (xof >= (L + k)) 

if (xof < gx) 
yof = f3b; 

end 
end 

if (xof >= gx) 

if (xof < hx) 

yof = gy + m2*(xof - gx) ; 

end 
end 

if (xof >= hx) 

if (xof < ox) 
yof = f4a; 

end 
end 

if (xof >= ox) 

if (xof < ix) 
yof = f4b; 

end 
end 

if (xof >= ix) 

if (xof < jx) 

yof = -by + ml*(xof - ix); 

end 
end 

if (xof == jx) 
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yof =0.0; 
end 

terml = cos(pi*iil*xof/L) ; 
term2 = sin(pi*iil*xof/L) ; 
yof 2 = -l*yof; 
dxldyl = dxl*yof2; 
exprl = terml*dxldyl; 
expr2 = term2*dxldyl; 
sumll = sumll + exprl; 
suml2 = suml2 + expr2; 
end 

asubn = (l/l_)*sumll; 
bsubn = (1/L)*suml2; 
if (iil == 0) 

asubn = asubn/2 ; 

bsubn = 0.0; 
end 

foura = asubn*cos(pi*iil*xvar/L) ; 
fourb = bsubn*cos(pi*iil*yvar/L) ; 
f pri ntf (fi d , ' N = '); 
f pri ntf (fid, '%4i ' ,iil); 
f pri ntf (fid, '\n'); 
f pri ntf (fid, 'a Term = '); 
f pri ntf (fid, '%24.20f ' , asubn); 
fprintf(fid, '*'); 
fprintf(fid, 'cos('); 
f pri ntf (fid, '%8.7f ' ,pi) ; 
fprintf(fid, '*'); 
f pri ntf (fid, '%li ' ,iil); 
fprintf(fid, '*x/') ; 
fprintf(fid, '%9.8f ' ,L); 
f pri ntf (fid, 'An') ; 
f pri ntf (fid, 'B Term = '); 
fprintf(fid, '%24.20f ' , bsubn); 
fprintf(fid, ■*'); 
fprintf(fid, 'sin('); 
fprintf(fid, '%8. 7f ' , pi); 
fprintf(fid, '*'); 
f pri ntf (fid, "%li ' ,iil); 
fprintf(fid, ^x/ 1 ); 
fprintf(fid, '%9.8f ' ,L); 
fprintf(fid, ')\n'); 
f pri ntf (fid, '\n'); 
end 

fclose(fid); 



coilwave3.ni 

% 



%coilwave3.m - Program to Compute Amplitude Values for Oscillators of Coil #3 of 
Magnetic Bottle 

% - waveform is Specified by interpolation of 200 term fourier series in 

Matlab 

% 



syms ax ay bx by cx cy dx dy ex ey fx fy; 
syms gx gy hx hy ix iy jx jy kx ky lx ly; 

syms mx my nx ny ox oy xof xofa yof yof2 ml m2 dxl dxldyl xvar yvar; 

syms al bl cl fl a2 b2 c2 a3 b3 c3 sumll sum! 2 terml term2; 

syms a4 b4 c4 ala bla cla fla alb bib clb exprl expr2 pi ; 

syms ale blc clc fib flc f2a f2b alq blq clq fq asubn bsubn; 

syms L k f3a f3b f4a f4b f4c x il i2 kl loops iiil foura fourb; 

pi = 3.141592654; 

L = 10; 

k = 2; 

kl = 0.25; 

ml = 2.5; 

m2 = -0.5; 

ax = 0; 

bx = k - kl + ax; 

cx = k + kl + ax; 

dx = L - k - kl + ax; 

ex = L - k + kl + ax; 

fx = L + k - kl + ax; 

gx = l_ + k + kl + ax; 

hx = 3*L - gx; 

ix = 3*L - fx; 

jx = 2*L; 

kx = k + ax; 

lx = 2*L - k + ax; 

mx = L - k + ax; 

nx = L + k + ax; 

ox = 2*L - k; 

ay = 0; 

by = ml*(k - kl) + ay; 

cy = m2*(kl) + ml*k + ay; 

dy = m2*(L - 2*k - kl) + ml*k + ay; 

ey = m2*(L - k + kl) + ml*k + ay; 

fy = m2*(i_ + k - kl) + ml*k + ay; 

gy = m2*(k + kl) + ay; 

hy = gy; 

iy = fy; 

jy = 0; 

ky = ml*k + ay; 

ly = -ml*k + ay; 

my = m2*(L - 2*k) + ml*k + ay; 

ny = -l*(dy+ey)/2 ; 

oy = -my; 

fid = f open ('g:/drawl/coilwave3. out 1 , 'w 1 ) ; 
fprintf (fid, 'Fourier Series of Coil #3 WaveForm 1 ); 
fprintf(fid, 'Xn 1 ); 

fprintf(fid, 1 '); 

fprintfCfid, 'Xn 1 ); 
fprintf(fid, 'Xn'); 
dxl = 2*L/1000; 
for iil = 0:100 

sumll = 0.0; 

suml2 = 0.0; 

for iiil = 1:500 
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xofa = iiil*2*L/1000; 

xof = 2*L - xofa; 

% 

%fl 

%2al*bx + bl = ml 
%2al*cx + bl = m2 
%cl = -al*bxA2 - bl*bx 

% 

al = (ml - m2)/(2*bx - 2*cx) ; 

bl = (2*bx*m2 - 2*cx*ml)/(2*bx - 2*cx) ; 

cl = (-l*al)*bx*bx - bl*bx; 

fl = al*xof*xof + bl*xof + cl; 

% 

%fla 

%is fl flipped 

% 

fla = al*(L - xof)*(l_ - xof) + bl*(l_ - xof) + cl; 

% 

%f2a 
%dy+fla 

% 

f2a = fla + dy; 

% 

%f2b 
%my+fla 

% 

f2b = fla + my; 

% 

%flb 

%-fl of (xof + L) 



alb = (-ml + m2)/(2*fx - 2*gx) ; 

bib = (2*fx*(-m2) - 2*gx*(-ml))/(2*fx - 2*gx) ; 

clb = (-l*alb)*fx*fx - blb*fx; 

fib = alb*xof*xof + blb*xof + clb; 

ale = (ml - m2)/(2*hx - 2*ix); 

blc = (2*hx*(m2) - 2*ix*(ml))/(2*hx - 2*ix); 
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clc = (-l*alc)*hx*hx - blc*hx; 

% 

%f3a 
%fy-flb 

% 

f3a = f y - fib; 

% 

%f3b 
Xny-flb 

% 

f3b = ny - fib; 

% 

%f4c 

%-fla at 2nd point 

% 

f4c = -alc*(xof)*(xof) - blc*(xof) - clc; 

% 

%f4a 

%-cy + f4c 

% _ 

f4a = -f4c - cy; 

% 

%f4b 

%f4c -cy 

% 

f4b = -f4c - cy; 

% 

%returns 

%x,y for x, cycle # 

% 

if (xof < bx) 
yof = ay + ml*xof ; 
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end 

if (xof >= bx) 

if (xof < cx) 

yof = by + al*xof*xof + bl*xof + cl; 

end 
end 

if (xof >= cx) 

if (xof < dx) 

yof = cy + m2*(xof - cx) ; 

end 
end 

if (xof >= dx) 

if (xof < (L - k)) 
yof = f2a; 

end 
end 

if (xof >= (L - k)) 

if (xof < ex) 
yof = f2b; 

end 
end 

if (xof >= ex) 

if (xof < fx) 

yof = ey + m2*(xof - ex); 

end 
end 

if (xof >= fx) 

if (xof < (L + k)) 
yof = f3a; 

end 
end 

if (xof >= (L + k)) 

if (xof < gx) 
yof = f3b; 

end 
end 

if (xof >= gx) 

if (xof < hx) 

yof = gy + m2*(xof - gx) ; 

end 
end 

if (xof >= hx) 

if (xof < ox) 
yof = f4a; 

end 
end 

if (xof >= ox) 

if (xof < ix) 
yof = f4b; 

end 
end 

if (xof >= ix) 

if (xof < jx) 

yof = -by + ml* (xof - ix); 

end 
end 

if (xof == jx) 

yof =0,0; 
end 

terml = cos(pi*iil*xof/L) ; 
term2 = sin(pi*iil*xof/L) ; 
dxldyl = dxl-yof ; 
exprl = terml*dxldyl; 

Page 4 



coilwave3 .m 

expr2 = term2*dxldyl; 
sumll = sumll + exprl; 
sum12 = suml2 + expr2; 
end 

for -nil = 501:1000 

xofa = iiil*2*L/1000; 

xof = xofa - 2*L; 
% 

%fl 

%2al*bx + bl = ml 
%2al*cx + bl = ml 
%cl = -al*bxA2 - bl*bx 

% 

al = (ml - m2)/(2*bx - 2*cx) ; 

bl = (2*bx*m2 - 2*cx*ml)/(2*bx - 2*cx); 

cl = (-l*al)*bx*bx - bl*bx; 

fl = al*xof*xof + bl*xof + cl; 



%fla 

%is fl flipped 

% 

fla = al*(L - xof)*(L - xof) + bl*(L - xof) + cl; 

% 

%f2a 
%dy+fla 

% 

f2a = fla + dy; 

% 

%f2b 
%my+fla 

% 

f2b = fla + my; 

% 

%flb 

%-fl of (xof + L) 

% 

alb = (-ml + m2)/(2*fx - 2*gx) ; 
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bib = (2*fx*(-m2) - 2*gx*(-ml))/(2*fx - 2*gx); 

clb = (-l*alb)*fx*fx - blb*fx; 

fib = alb*xof*xof + blb*xof + clb; 

ale = (ml - m2)/(2*hx - 2*ix); 

blc = (2*hx*(m2) - 2*ix*(ml))/(2*hx - 2*ix); 

clc = (-l*alc)*hx*hx - blc*hx; 

% 

%f3a 
%fy-flb 

% 

f3a = f y - fib; 

% 

%f3b 
%ny-flb 

% 

f3b = ny - fib; 

% 

%f4c 

%-fla at 2nd point 

% 

f4c = -alc*(xof)*(xof) - blc*(xof) - clc; 

% 

%f4a 

%-cy + f4c 

% 

f4a = -f4c - cy; 

% 

%f4b 

%f4c -cy 

% 

f4b = -f4c - cy; 

% 

%returns 

%x,y for x, cycle # 
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coilwave3.m 

% 



if (xof < bx) 

yof = ay + ml*xof ; 
end 

if (xof >= bx) 

if (xof < cx) 

yof = by + al*xof*xof + bl*xof + cl; 

end 
end 

if (xof >= cx) 

if (xof < dx) 

yof = cy + m2*(xof - cx) ; 

end 
end 

if (xof >= dx) 

if (xof < (L - k)) 
yof = f2a; 

end 
end 

if (xof >= (L - k)) 

if (xof < ex) 
yof = f2b; 

end 
end 

if (xof >= ex) 

if (xof < fx) , 
yof = ey + m2*(xof - ex); 

end 
end 

if (xof >= fx) 

if (xof < (L + k)) 
yof = f3a; 

end 
end 

if (xof >= (L + k)) 

if (xof < gx) 
yof = f3b; 

end 
end 

if (xof >= gx) 

if (xof < hx) 

yof = gy + m2*(xof - gx) ; 

end 
end 

if (xof >= hx) 

if (xof < ox) 
yof = f4a; 

end 
end 

if (xof >= ox) 

if (xof < ix) 
yof = f4b; 

end 
end 

if (xof >= ix) 

if (xof < jx) 

yof = -by + ml*(xof - ix); 

end 
end 

if (xof == jx) 
yof = 0.0; 
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end 

terml = cos(pi*iil*xof/L) ; 
term2 = sin(pi*iil*xof/L) ; 
yof2 = l*yof; 
dxldyl = dxl*yof2; 
exprl = terml*dxldyl; 
expr2 = term2*dxldyl; 
sumll = sumll + exprl; 
sum12 = suml2 + expr2; 
end 

asubn = (1/L)*sum"ll; 
bsubn = (1/L)*suml2; 
if (iii == 0) 

asubn = asubn/2; 

bsubn = 0.0; 
end 

foura = asubn*cos(pi *iil*xvar/L) ; 

fourb = bsubn*cos(pi *iil*yvar/L) ; 

fpri ntf (fid , 1 N = '); 

fprintf(fid, '%4i ' ,iil); 

fpri ntf (fid, '\n'); 

fpri ntf (fid, 'A Term = '); 

fprintf(fid, '%24.20f ' , asubn); 

fpri ntf (fid, '*'); 

fpri ntf (fid, 'cos('); 

fprintf(fid,'%8.7f , ,pi); 

fpri ntf (fid, '*'); 

fprintf(fid,'%li , ,iil); 

fprintf(fid,'*x/'); 

fpri ntf (fid, '%9.8f ,L); 

fprintf(fid, *)\n'); 

fpri ntf (fid , 1 B Term = '); 

fprintf(fid, '%24.20f , bsubn) ; 

fprintf(fid, '*'); 

fpri ntf(fid , 'sinC ); 

fprintf(fid, '%8.7f ' ,pi); 

fprintf(fid, '*'); 

fprintf(fid,'%li , ,iil); 

fprintf(fid, , *x/'); 
fprintf(fid, '%9.8f ' ,L); 
fprintf(fid, ')\n'); 
fpri ntf (fid, '\n'); 
end 

fclose(fid) ; 



coilwave2.m 

% 



%coilwave2.m - Program to Compute Amplitude Values for Oscillators of Coil #2 of 
Magnetic Bottle 

% - waveform is Specified by interpolation of 200 term fourier series in 

Matlab 

% 



syms ax ay bx by cx cy dx dy ex ey fx fy; 
syms gx gy hx hy ix iy jx jy kx ky lx ly; 

syms mx my nx ny ox oy xof xofa yof yof2 ml m2 dxl dxldyl xvar yvar; 

syms al bl cl fl a2 b2 c2 a3 b3 c3 sumll suml2 terml term2; 

syms a4 b4 c4 ala bla cla fla alb bib clb exprl expr2 pi ; 

syms ale blc clc fib flc f2a f2b alq blq clq fq asubn bsubn; 

syms L k f3a f3b f4a f4b f4c x il i2 kl loops iiil foura fourb; 

pi = 3.141592654; 

L = 10; 

k = 2; 

kl = 0.25; 

ml = 2.5; 

m2 = -0.5; 

ax = 0; 

bx = k - kl + ax; 

cx = k + kl + ax; 

dx = L - k - kl + ax; 

ex = L - k + kl + ax; 

fx = L + k - kl + ax; 

gx = L + k + kl + ax; 

hx = 3*L - gx; 

ix = 3*L - fx; 

ix = 2*l; 

kx = k + ax; 

lx = 2*L - k + ax; 

mx = L - k + ax; 

nx = L + k + ax; 

ox = 2*L - k; 

ay = 0; 

by = ml' v (k - kl) + ay; 

cy = m2*(kl) + ml*k + ay; 

dy = m2*(L - 2*k - kl) + ml*k + ay; 

ey = m2*(L - k + kl) + ml*k + ay; 

fy = m2*(L + k - kl) + ml*k + ay; 

gy = m2*(k + kl) + ay; 

hy = gy; 

iy = fy; 

jy = 0; 

ky = ml*k + ay; 

ly = -ml*k + ay; 

my = m2*(L - 2*k) + ml*k + ay; 

ny = -l*(dy+ey)/2; 

oy = -my; 

fid = f open ( , g:/drawl/coilwave2. out 1 , 'w 1 ) ; 
fprintf(fid, 'Fourier Series of Coil #2 waveForm 1 ); 
fprintf(fid, An 1 ); 

fprintf(fid, 1 '); 

fprintfCfid/Xn 1 ); 
fprintf(fid, An 1 ); 
dxl = 2*L/1000; 
for iil = 0:100 

sumll = 0.0; 

sum!2 = 0.0; 

for iiil = 1:500 

Page 1 



coilwave2 .m 

xofa = iiil*2*L/1000; 
xof = L - xofa; 



%fl 

%2al*bx + bl = ml 
%2al*cx + bl = m2 
%cl = -al*bxA2 - bl*bx 
% 



al = (ml - m2)/(2*bx - 2*cx) ; 

bl = (2*bx*m2 - 2*cx*ml)/(2*bx - 2*cx); 

cl = (-l*al)*bx*bx - bl*bx; 

fl = al*xof*xof + bl*xof + cl; 

% 

%fla 

%is fl flipped 



fla = al*(L - xof)*(l_ - xof) + bl*(L - xof) + cl; 

% 

%f2a 
%dy+fla 

% 

f2a = fla + dy; 

% 

%f2b 
%my+fla 

% 

f2b = fla + my; 

%— _ 

%flb 

%-fl of (xof + L) 

% 

alb = (-ml + m2)/(2*fx - 2*gx); 

bib = (2*fx*(-m2) - 2*gx*(-ml))/(2*fx - 2*gx); 

clb = (-l*alb)*fx*fx - blb*fx; 

fib = alb*xof*xof + blb*xof + clb; 

ale = (ml - m2)/(2*hx - 2*ix); 

blc = (2*hx*(m2) - 2*ix*(ml))/(2*hx - 2*ix); 
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clc = (-l*alc)*hx*hx - blc*hx; 



% 

%f3a 
%fy-flb 

% 

f3a = f y - fib; 

% 

%f3b 
%ny-flb 

%___ 

f3b = ny - fib; 

% 

%f4c 

%-fla at 2nd point 

% 

f4c = -alc*(xof)*Cxof) - blc*(xof) - clc; 

% 

%f4a 

%-cy + f4c 

% 

f4a = -f4c - cy; 

% 

%f4b 
%f4c -cy 

% 

f4b = -f4c - cy; 

% 

^returns 

%x,y for x, cycle # 

% 

if (xof < bx) 

yof = ay + ml*xof ; 
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end 

if (xof >= bx) 

if (xof < cx) 

yof = by + al*xof*xof + bl*xof + cl; 

end 
end 

if (xof >= cx) 

if (xof < dx) 

yof = cy + m2*(xof - cx) ; 

end 
end 

if (xof >= dx) 

if (xof < (L - k)) 
yof = f2a; 

end 
end 

if (xof >= (L - k)) 

if (xof < ex) 
yof = f2b; 

end 
end 

if (xof >= ex) 

if (xof < fx) 

yof = ey + m2*(xof - ex); 

end 
end 

if (xof >= fx) 

if (xof < (L + k)) 
yof = f3a; 

end 
end 

if (xof >= (L + k)) 

if (xof < gx) 
yof = f3b; 

end 
end 

if (xof >= gx) 

if (xof < hx) 

yof = gy + m2*(xof - gx) ; 

end 
end 

if (xof >= hx) 

if (xof < ox) 
yof = f4a; 

end 
end 

if (xof >= ox) 

if (xof < ix) 
yof = f4b; 

end 
end 

if (xof >= ix) 

if (xof < jx) 

yof = -by + ml*(xof - ix); 

end 
end 

if (xof == jx) 

yof =0.0; 
end 

terml = cos(pi*iil*xof/l_) ; 
term2 = sin (pi *iil*xof/L) ; 
yof 2 = -l*yof; 
dxldyl = dxl*yof2; 
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exprl = terml*dxldyl; 
expr2 = term2*dxldyl; 
sumll = sumll + exprl; 
suml2 = sum! 2 + expr2; 
end 

for iiil = 501:1000 
xofa = i-iil*2*L/1000; 

xof = xofa - L ; 

% 

%fl 

%2al*bx + bl = ml 
%2al*cx + bl = m2 
%cl = -al*bxA2 - bl*bx 

% 

al = (ml - m2)/(2*bx - 2*cx) ; 

bl = (2*bx*m2 - 2*cx*ml)/(2*bx - 2*cx) ; 

cl = (-l*al)*bx*bx - bl*bx; 

fl = al*xof*xof + bl*xof + cl; 

o/ 0 

%fla 

%is fl flipped 

% 

fla = al*(L - xof)*(L - xof) + bl*(L - xof) + cl; 

% 

%f2a 
%dy+fla 

% 

f2a = fla + dy; 

% 

%f2b 
%my+fla 

% 

f2b = fla + my; 

% 

%flb 

%-fl of (xof + L) 

% 
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alb = (-ml + m2)/(2*fx - 2*gx); 

bib = (2*fx*(-m2) - 2*gx*(-ml))/(2*f x - 2*gx); 

clb = (-l*alb)*fx*fx - blb*fx; 

fib = alb*xof*xof + blb*xof + clb; 

ale = (ml - m2)/(2*hx - 2*ix); 

blc = (2*hx*(m2) - 2*ix*(ml))/(2*hx - 2*ix); 

clc = (-l*alc)*hx*hx - blc*hx; 

% 

%f3a 
%fy-flb 

% 

f3a = fy - fib; 

% 

%f3b 
%ny-flb 

% 

f3b = ny - fib; 

% 

%f4c 

%-fla at 2nd point 

% 

f4c = -alc*(xof)*(xof) - blc*(xof) - clc; 
% 

%f4a 

%-cy + f4c 

% 

f4a = -f4c - cy; 

% 

%f4b 

%f4c -cy 

% 

f4b = -f4c - cy; 

% 

%returns 
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%x,y for x, cycle # 



if (xof < bx) 

yof = ay + ml*xof ; 
end 

if (xof >= bx) 

if (xof < cx) 

yof = by + al*xof*xof + bl*xof + cl; 

end 
end 

if (xof >= cx) 

if (xof < dx) 

yof = cy + m2*(xof - cx) ; 

end 
end 

if (xof >= dx) 

if (xof < (L - k)) 
yof = f2a; 

end 
end 

if (xof >= (L - k)) 

if (xof < ex) 
yof = f2b; 

end 
end 

if (xof >= ex) 

if (xof < fx) 

yof = ey + m2*(xof - ex); 

end 
end 

if (xof >= fx) 

if (xof < (L + 10) 
yof = f3a; 

end 
end 

if (xof >= (L + k)) 

if (xof < gx) 
yof = f3b; 

end 
end 

if (xof >= gx) 

if (xof < hx) 

yof = gy + m2*(xof - gx) ; 

end 
end 

if (xof >= hx) 

if (xof < ox) 
yof = f4a; 

end 
end 

if (xof >= ox) 

if (xof < ix) 
yof = f4b; 

end 
end 

if (xof >= ix) 

if (xof < jx) 

yof = -by + ml* (xof - ix); 

end 
end 

if (xof == jx) 
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%x,y for x, cycle # 

% 



if (xof < bx) 

yof = ay + ml*xof ; 
end 

if (xof >= bx) 

if (xof < cx) 

yof = by + al*xof*xof + bl*xof + cl; 

end 
end 

if (xof >= cx) 

if (xof < dx) 

yof = cy + m2*(xof - cx) ; 

end 
end 

if (xof >= dx) 

if (xof < (L - k)) 
yof = f2a; 

end 
end 

if (xof >- (L - k)) 

if (xof < ex) 
yof = f2b; 

end 
end 

if (xof >= ex) 

if (xof < fx) 

yof = ey + m2*(xof - ex); 

end 
end 

if (xof >= fx) 

if (xof < (L + k)) 
yof = f3a; 

end 
end 

if (xof >= (L + k)) 

if (xof < gx) 
yof = f3b; 

end 
end 

if (xof >= gx) 

if (xof < hx) 

yof = gy + m2*(xof - gx) ; 

end 
end 

if (xof >= hx) 

if (xof < ox) 
yof = f4a; 

end 
end 

if (xof >= ox) 

if (xof < ix) 
yof = f4b; 

end 
end 

if (xof >= ix) 

if (xof < jx) 

yof = -by + ml*(xof - ix); 

end 
end 

if (xof == jx) 
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yof = 0.0; 
end 

terml = cos (pi *i il*xof/L) ; 
term2 = sin (pi *i il*xof/L) ; 
yof 2 = -l*yof ; 
dxldyl = dxl*yof2; 
exprl = terml*dxldyl; 
expr2 = term2*dxldyl; 
sumll = sumll + exprl; 
suml2 = suml2 + expr2; 
end 

asubn = (1/L)*sumll; 
bsubn = (1/L)*suml2; 
if (iii == 0) 

asubn = asubn/2; 

bsubn = 0.0; 
end 

foura = asubn*cos(pi*iil*xvar/L) ; 
fourb = bsubn*cos(pi*iil*yvar/L) ; 
fprintf (fid, 'N = '); 
fprintf(fid, '%4i ' ,iil); 
fprintf(fid, '\n'); 
fprintf (fid, 'A Term = '); 
fprintf(fid, , %24.20f , asubn); 
fprintf (fid, '*') ; 
fprintf (fid, 'cos(') ; 
fprintf (fid, '%8.7f ' ,pi); 
fprintf(fid, , *') ; 
fprintf(fid, '%li ' ,iil); 
fprintf (fid, '*x/'); 
fprintf(fid, '%9.8f ' ,L); 
fprintf (fid, ')\n') ; 
fpri ntf (f id , 1 B Term = '); 
fprintf(fid, '%24.20f ' , bsubn); 
fprintf(fid, '*') ; 
fprintf(fid, , sin('); 
fprintf(fid,'%8.7f,pi); 
fpri ntf (fid, '*') ; 
fpri ntf (fid, '%li ' ,iil) ; 
fpri ntf (fid, '*x/') ; 
fpri ntf (fid, '%9.8f ' ,L); 
fprintf(fid, ')\n'); 
fpri ntf (fid, '\n'); 
end 

fclose(fid); 
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% 



%coilwavel.m - Program to Compute Amplitude values for Oscillators of Coil #1 of 
Magnetic Bottle 

% " - Waveform is Specified by Interpolation of 200 term fourier series in 

Matlab 

% 



syms ax ay bx by cx cy dx dy ex ey fx fy; 

syms gx gy hx hy ix iy jx jy kx ky lx ly; 

syms mx my nx ny ox oy xof yof ml m2 dxl dxldyl xvar yvar; 

syms al bl cl fl a2 b2 c2 a3 b3 c3 sumll suml2 terml term2; 

syms a4 b4 c4 ala bla cla fla alb bib clb exprl expr2 pi ; 

syms ale blc clc fib flc f2a f2b alq blq clq fq asubn bsubn; 

syms L k f3a f3b f4a f4b f4c x il i2 kl loops iiil foura fourb; 

pi = 3.141592654; 

L = 10; 

k = 2; 

kl = 0.25; 

ml = 2.5; 

m2 = -0.5; 

ax = 0; 

bx = k - kl + ax; 

cx = k + kl + ax; 

dx = L - k - kl + ax; 

ex = L - k + kl + ax; 

fx = L + k - kl + ax; 

gx = L + k + kl + ax; 

hx = 3*L - gx; 

ix = 3*L - fx; 

ix = 2*L; 

kx = k + ax; 

lx = 2*L - k + ax; 

mx = L - k + ax; 

nx = L + k + ax; 

ox = 2*L - k; 

ay = 0; 

by = ml*(k - kl) + ay; 

cy = m2*(kl) + ml*k + ay; 

dy = m2*(L - 2*k - kl) + ml*k + ay; 

ey = m2*(L - k + kl) + ml*k + ay; 

fy = m2*(L + k - kl) + ml*k + ay; 

gy = m2*(k + kl) + ay; 

hy = gy; 

iy = fy; 

jy = 0; 

ky = ml*k + ay; 

ly = -ml*k + ay; 

my = m2*(L - 2*k) + ml*k + ay; 

ny = -l*(dy+ey)/2; 

oy = -my; 

fid = fopenC'gr/drawl/coilwavel.out 1 , 'w 1 ) ; 
fprintf (fid, 'Fourier Series of Coil #1 waveForm 1 ); 
fprintf(fid, An 1 ); 

fprintf(fid, 1 '); 

fprintfCfid/Xn 1 ); 
fprintf(fid, An 1 ); 
dxl = 2*L/1000; 
for iil = 0:100 

sumll = 0.0; 

suml2 = 0.0; 

for iiil = 1:1000 
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xof = iiil*2*L/1000; 

% 



coilwavel.m 



%fl 

%2al*bx + bl = ml 
%2al*cx + bl = m2 
%cl = -al*bxA2 - bl*bx 

% 

al = (ml - m2)/(2*bx - 2*cx) ; 

bl = (2*bx*m2 - 2*cx*ml)/(2*bx - 2*cx) ; 

cl = (-l*al)*bx*bx - bl*bx; 

fl = al*xof*xof + bl*xof + cl; 

% 

%fla 

%is fl flipped 

% 

fla = al*(L - xof)*(L - xof) + bl*(L - xof) + cl; 

% 

%f2a 
%dy+fla 

% 

f2a = fla + dy; 

% 

%f2b 
%my+fla 

% _ 

f2b = fla + my; 

% 

%flb 

%-fl of (xof + L) 



alb = (-ml + m2)/(2*fx - 2*gx); 

bib = (2*fx*(-m2) - 2*gx*(-ml))/(2*fx - 2*gx) ; 

clb = (-l*alb)*fx*fx - blb*fx; 

fib = alb*xof*xof + blb*xof + clb; 

ale = (ml - m2)/(2*hx - 2*ix); 

blc = (2*hx*(m2) - 2*ix*(ml))/(2*hx - 2*ix); 

clc = (-l*alc)*hx*hx - blc*hx; 
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% 

%f3a 
%fy-flb 

% 

f3a = f y - fib; 

% 

%f3b 
%ny-flb 

% 

f3b = ny - fib; 

% 

%f4c 

%-fla at 2nd point 

% _ 

f4c = -alc*(xof)*(xof) - blc*(xof) - clc; 

% 

%f4a 

%-cy + f4c 

% 

f4a = -f4c - cy; 

% 

%f4b 
%f4c -cy 

% 

f4b = -f4c - cy; 

% 

%returns 

%x,y for x, cycle # 

% 

if (xof < bx) 

yof = ay + ml*xof ; 
end 
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if (xof >= bx) 

if (xof < cx) 

yof = by + al*xof*xof + bl*xof + cl; 

end 
end 

if (xof >= cx) 

if (xof < dx) 

yof = cy + m2*(xof - cx) ; 

end 
end 

if (xof >= dx) 

if (xof < (L - k)) 
yof = f2a; 

end 
end 

if (xof >= (L - k)) 

if (xof < ex) 
yof = f2b; 

end 
end 

if (xof >= ex) 

if (xof < fx) 

yof = ey + m2*(xof - ex); 

end 
end 

if (xof >= fx) 

if (xof < (L + k)) 
yof = f3a; 

end 
end 

if (xof >= (L + k)) 

if (xof < gx) 
yof = f3b; 

end 
end 

if (xof >= gx) 

if (xof < hx) 

yof = gy + m2*(xof - gx) ; 

end 
end 

if (xof >= hx) 

if (xof < ox) 
yof = f4a; 

end 
end 

if (xof >= ox) 

if (xof < ix) 
yof = f4b; 

end 
end 

if (xof >= ix) 

if (xof < jx) 

yof = -by + ml*(xof - ix); 

end 
end 

if (xof == jx) 

yof =0.0; 
end 

terml = cos(pi*iil*xof/L) ; 
term2 = sin(pi*iil*xof/L) ; 
dxldyl = dxl*yof ; 
exprl = terml*dxldyl; 
expr2 = term2*dxldyl; 
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sumll = sumll + exprl; 
suml2 = sum12 + expr2; 
end 

asubn = (1/L)*suml 1; 
bsubn = (1/L)*suml2; 
if (iil == 0) 

asubn = asubn/2; 

bsubn = 0.0; 
end 

foura = asubn*cos(pi*iil*xvar/L); 
fourb = bsubn*cos(pi*iil*yvar/L); 
fprintfCfid.'N = '); 
fprintf(fid, '%4i ' ,iil); 
fprintfCfid.'Xn'); 
fprintf (fid, 'A Term = '); 
fprintfCfid, '%24.20f , asubn); 
fprintf (fid, '*'); 
fprintf (fid.'cosC); 
fprintf(fid,'%8.7f ,pi); 
fprintf(fid,'*'); 
fprintf (fid, '%li ' ,iil); 
fprintf(fid,'*x/'); 
fprintf(fid, '%9.8f ' ,L); 
fprintf (fid, 'An'); 
fprintf (fid, 'B Term = '); 
fprintf (fid, '%24.20f ' , bsubn) ; 
fprintf(fid,'*'); 
fprintf (fid, 'sinO ; 
fprintf (fid, '%8.7f ' ,pi); 
fprintf (fid, '*'); 
fprintf (fid, '%li ' ,iil); 

fprintf(fid, , *x/'); 
fprintf (fid, '%9.8f ' ,L); 
fprintf(fid,')\n'); 
fprintf (fid, '\n'); 
end 

fclose(fid) ; 



Fourier Series of Coil #1 waveForm 



coilwavel.out 



N = 0 
A Term = 
B Term = 

N = 1 
A Term = 
B Term = 

N = 2 
A Term = 
B Term = 

N = 3 
A Term = 
B Term = 

N = 4 
A Term = 
B Term = 

N = 5 
A Term = 
B Term = 

N = 6 
A Term = 
B Term = 

N = 7 
A Term = 
B Term = 

N = 8 
A Term = 
B Term = 

N = 9 
A Term = 
B Term = 

N = 10 
A Term = 
B Term = 

N = 11 
A Term = 
B Term = 

N = 12 
A Term = 
B Term = 

N = 13 
A Term = 
B Term = 

N = 14 

A Term = 
B Term = 



0 . 32151000000000080000*cos 
0.00000000000000000000*sin 



-0.05486015000563282400*cos 
3 . 59960771976998780000*si n 



-0 . 28606410623826406000*cos 
1 . 31790776415870600000*si n 



0.02127969436458058600*cos 
0.68364708872142377000*sin 



-0 . 12370743263194371000*cos 
0 . 32864714547868346000*si n 



0.06624931097526866300*cos 
-0.00202007584984204580*sin 



0.02232680536410826600*cos 
-0.00262224676588311510*sin 



0.01892140443709789500*cos 
-0 . 15843 526909556116000*si n 



0.07502481994225042300*cos 
-0.08235904454089781900*sin 



-0.05130372995787575700*cos 
-0.06472085885657588600*sin 



0.03963932407652325000*cos 
-0 . 04893408424081575200*si n 



-0.04741530329485993400*cos 
0.05525853461247107900*sin 



-0.01948282817234723400*cos 
0.00185072898389278590*sin 



0.01923342680901290300*cos 
0 .06647372279285850000*si n 



-0.04283176627266559900*cos 
0.02297951841575042500*sin 



3 . 1415927*0*x/10 . 00000000) 
3 . 1415927*0*x/10 . 00000000) 



3 . 1415927*l*x/10 . 00000000) 
3 . 1415927*l*x/10 . 00000000) 



3 . 1415927*2*x/10 .00000000) 
3 . 1415927*2*x/10 . 00000000) 



3 . 1415927*3*x/10 . 00000000) 
3 . 1415927*3*x/10 . 00000000) 



3 . 1415927*4*x/10 . 00000000) 
3 . 1415927*4*x/10 . 00000000) 



3 . 1415927*5*x/10 .00000000) 
3 . 1415927*5*x/10 .00000000) 



3 . 1415927*6*x/10 . 00000000) 
3 . 1415927*6*x/10 . 00000000) 



3 . 1415927*7*x/10 .00000000) 
3 . 1415927*7*x/10 . 00000000) 



3 . 1415927*8*x/10 .00000000) 
3 . 141592 7*8*x/10 .00000000) 



3 . 1415927*9*x/10 . 00000000) 
3. 1415927*9*x/10. 00000000) 



3 . 141592 7 -10*x/10 . 00000000) 
3 . 141592 7*10*x/10 . 00000000) 



3 . 141592 7*ll*x/10 . 00000000) 
3 . 1415927*ll*x/10 . 00000000) 



3 . 1415927*12*x/10 . 00000000) 
3. 1415927*12*x/10. 00000000) 



3. 1415927*13*x/10. 00000000) 
3. 1415927*13*x/10. 00000000) 



3 . 1415927*14*x/10 . 00000000) 
3 . 1415927*14*x/10 . 00000000) 
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I 



N = 15 
A Term = 
B Term = 

N = 16 
A Term = 
B Term = 

N = 17 
A Term = 
B Term = 

N = 18 
A Term = 
B Term = 

N = 19 
A Term = 
B Term = 

N = 20 
A Term = 
B Term = 

N = 21 
A Term = 
B Term = 

N = 22 
A Term = 
B Term = 

N = 23 
A Term = 
B Term = 

N = 24 
A Term = 
B Term = 

N = 25 
A Term = 
B Term = 

N = 26 
A Term = 
B Term = 

N = 27 
A Term = 
B Term = 

N = 28 
A Term = 
B Term = 

N = 29 
A Term = 
B Term = 

N = 30 
A Term = 
B Term = 



coil wavel. out 



0.05318462413540643600*cos 
-0.00555508117389191380*sin 



-0.01966295136883588000*cos 
0.01622586292254936700*sin 



0 . 0132 5992646096467800*cos 
-0.05034501867654114800*sin 



0.01739919222397608500*cos 
O.O0231256317432152740*sin 



-0.03784310796783909400*cos 
-0.02026438592921991900*sin 



0.03054541665948651300*cos 
-0.00685818717707770060*sin 



-0.03175351239314501400*cos 
0 . 02760073001586024400*si n 



0 . 01097494260098336200*cos 
-0.00998983404396259350*sin 



0.01326094124579632100*cos 
0 . 02764880007382740700*si n 



-0.01736471933434124600*cos 
-0.00639988620723064890*sin 



O.03163208011535725100*cos 
-0.00775663763005859670*sin 



-0.02432662006972984000*cos 
O.003276131O661935O36O*sin 



0.00647780727335501610*cos 
-0.02273594482710212700*sin 



-0.00448567468792871430*cos 
0. 01081935869544950500*si n 



-O.01901986543122055300*cos 
-0.00321493589944499970*sin 



0.01847807299492187200*cos 
0.00680464415618746090*sin 



3 . 1415927*15*x/10 . 00000000) 
3 . 1415927*15*x/10 . 00000000) 



3 . 1415927*16*x/10 . 00000000) 
3 . 1415927*16*x/10 . 00000000) 



3 . 1415927*17*x/10 . 00000000) 
3 . 1415927*17*x/10 . 00000000) 



3 . 1415927*18*x/10 . 00000000) 
3 . 1415927*18*x/10 . 00000000) 



3 . 1415927*19*x/10 . 00000000) 
3 . 1415927*19*x/10 . 00000000) 



3. 1415927*20*x/10. 00000000) 
3 . 1415927*20*x/10 . 00000000) 



3 . 1415927*21*x/10 . 00000000) 
3 . 1415927*21*x/10 . 00000000) 



3 . 1415927*22*x/10 .00000000) 
3 . 1415927*22*x/10 . 00000000) 



3. 1415927*23*x/10. 00000000) 
3 . 1415927*23*x/10 .00000000) 



3 . 1415927*24*x/10 . 00000000) 
3 . 1415927*24*x/10 . 00000000) 



3. 1415927*25*x/10. 00000000) 
3. 1415927*25*x/10. 00000000) 



3. 1415927*26*x/10. 00000000) 
3 . 1415927*26*x/10 . 00000000) 



3 . 1415927*27*x/10 . 00000000) 
3 . 1415927*27*x/10 . 00000000) 



3 . 1415927*28*x/10 . 00000000) 
3 . 1415927*28*x/10 . 00000000) 



3 . 1415927*29*x/10 . 00000000) 
3 . 1415927*29*x/10 . 00000000) 



3. 1415927*30*x/10. 00000000) 
3 . 1415927*30*x/10 . 00000000) 
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( 



coil wavel. out 



N = 31 
A Term = 
B Term = 

N = 32 
A Term = 
B Term = 

N = 33 
A Term = 
B Term = 

N = 34 
A Term = 
B Term = 

N = 35 
A Term = 
B Term = 

N = 36 
A Term = 
B Term = 

N = 37 
A Term = 
B Term = 

N = 38 
A Term = 
B Term = 

N = 39 
A Term = 
B Term = 

N = 40 
A Term = 
B Term = 

N = 41 
A Term = 
B Term = 

N = 42 
A Term = 
B Term = 

N = 43 
A Term = 
B Term = 

N = 44 
A Term = 
B Term = 

N = 45 
A Term = 
B Term = 

N = 46 
A Term = 



-0.01335542965674995700*cos 
0 . 01397175417416619100*si n 



0 . 01895225373098448400*cos 
-0.00540025280609754060*sin 



O.00525239301261612970*cos 
0.00539726191049372510*sin 



-0.00181500168679699890*cos 
-0.01111043822350311700*sin 



0.00884024988157692920*cos 
-0.00826646205777484790*sin 



-0.01853781900409428100*cos 
-0 . 00326373777604168450*si n 



0.00090858332093911483*cos 
-0.00369006018012366060*si n 



-0.01260582756548626700*cos 
0.00766565536988162570*sin 



-0.00095468235656087977*cos 
0.00759183857054968050*sin 



0.00682897771037084020*cos 
0.00802591596683270920*sin 



0.00186672892665128720*cos 
0 . 00361371471043719660*si n 



0.01612278857333271400*cos 
-0.00087980358413492659*sin 



-0.00225531366458747200*cos 
-O.0O90188O27572160531O*sin 



0.00618230530849278250*cos 
-0.00653418431482273060*sin 



-0.00821891555902279710*cos 
-0.0O72906752868269545O*sin 



3 . 1415927*31*x/10 . 00000000) 
3 . 1415927*31*x/10 . 00000000) 



3 . 1415927*32*x/10 . 00000000) 
3 . 1415927*32*x/10 . 00000000) 



3 . 1415927*33*x/10 . 00000000) 
3 . 1415927*33*x/10 . 00000000) 



3 . 1415927*34*x/10 . 00000000) 
3 . 1415927*34*x/10 .00000000) 



3 . 1415927*3 5*x/10 . 00000000) 
3. 1415927*35*x/10. 00000000) 



3 . 1415927*36*x/10 . 00000000) 
3 . 1415927*36*x/10 . 00000000) 



3 . 1415927*37*x/10 . 00000000) 
3 . 1415927*37*x/10 . 00000000) 



3 . 1415927*38*x/10 . 00000000) 
3 . 1415927*38*x/10 . 00000000) 



3. 1415927*39*x/10. 00000000) 
3 . 1415927*39*x/10 . 00000000) 



3 . 1415927*40*x/10 . 00000000) 
3 . 1415927*40*x/10 . 00000000) 



3. 1415927*41*x/10. 00000000) 
3 . 1415927*41*x/10 . 00000000) 



3 . 1415927*42*x/10 . 00000000) 
3 . 1415927*42*x/10 . 00000000) 



3 . 1415927*43*x/10 . 00000000) 
3 . 1415927*43*x/10 . 00000000) 



3 . 1415927*44*x/10 . 00000000) 
3 . 1415927*44*x/10 . 00000000) 



3 . 1415927*45*x/10 . 00000000) 
3. 1415927*45*x/10. 00000000) 



-0 . 00900751434422850970*cos (3 . 1415927*46*x/10 . 00000000) 
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B Term = 

N = 47 
A Term = 
B Term = 

N = 48 
A Term = 
B Term = 

N = 49 
A Term = 
B Term = 

N = 50 
A Term = 
B Term = 

N = 51 
A Term = 
B Term = 

N = 52 
A Term = 
B Term = 

N = 53 
A Term = 
B Term = 

N = 54 
A Term = 
B Term = 

N = 55 
A Term = 
B Term = 

N = 56 
A Term = 
B Term = 

N = 57 
A Term = 
B Term = 

N = 58 
A Term = 
B Term = 

N = 59 
A Term = 
B Term = 

N = 60 
A Term = 
B Term = 

N = 61 
A Term = 
B Term = 

N = 62 



coilwavel.out 

-0.00337286820234028320*sin(3.1415927*46*x/10. 00000000) 



-0 . 00190034301289706460*cos 
0 . 00803873711547283850*si n 



-0.01187033516838867900*cos 
0 . 00186524040013030970*si n 



0 . 01110848179612070700*cos 
0 . 01241246934037173100*si n 



-0.00159117423189477580*cos 
0 . 00260416208199160420*si n 



0.00989529039500905950*cos 
-0.0025433136145555793O*sin 



0.00834195745433642490*cos 
0 . 00100061304297829000*si n 



-0 . 00709848487851465120*cos 
-0 . 01463862004224305100*si n 



0.00780453689952040340*cos 
0.00093723523788297267*sin 



-0.01552742091581945200*cos 
-0 . 00546954959031795660*si n 



-0.00031068263743843705*cos 
0.00045535291016017012*sin 



-0 . 00171980053120679780*cos 
0 . 01134967159167115000*si n 



-0.00650118847783573860*cos 
-0 . 00282284986098885150*si n 



0.01476205547024808100*cos 
0 .01152390686665279900*si n 



-0 . 00546928763965044780*cos 
-0 . 00449115806730165590*si n 



0 . 00998143469133775000*cos 
-0 . 00384198825224952960*si n 



3 . 1415927*47*x/10 . 00000000) 
3. 1415927*47*x/10. 00000000) 



3. 1415927*48*x/10. 00000000) 
3 . 1415927*48* x/10 . 00000000) 



3 . 1415927*49*x/10 . 00000000) 
3 . 1415927*49*x/10 . 00000000) 



3. 1415927*50*x/10. 00000000) 
3 . 1415927*50*x/10 . 00000000) 



3. 1415927*51*x/10. 00000000) 
3. 1415927*51*x/10. 00000000) 



3 . 1415927*52*x/10 . 00000000) 
3 . 1415927*52*x/10 . 00000000) 



3 . 1415927*5 3*x/10 . 00000000) 
3 . 1415927*5 3*x/10 . 00000000) 



3 . 1415927*54*x/10 . 00000000) 
3 . 1415927*54*x/10 . 00000000) 



3 . 1415927*5 5*x/10 . 00000000) 
3. 1415927*55*x/10. 00000000) 



3 . 1415927*56*x/10 . 00000000) 
3 . 1415927*56*x/10 . 00000000) 



3 . 1415927*57*x/10 . 00000000) 
3 . 1415927*57*x/10 . 00000000) 



3 . 1415927*58*x/10 . 00000000) 
3 . 1415927*58*x/10 . 00000000) 



3 . 1415927*59*x/10 . 00000000) 
3 . 1415927*59*x/10 . 00000000) 



3. 1415927*60*x/10. 00000000) 
3. 1415927*60*x/10. 00000000) 



3 . 1415927*61*x/10 . 00000000) 
3 . 1415927*61*x/10 . 00000000) 
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A Term = 
B Term = 

N = 63 
A Term = 
B Term = 

N = 64 
A Term = 
B Term = 

N = 65 
A Term = 
B Term = 

N = 66 
A Term = 
B Term = 

N = 67 
A Term = 
B Term = 

N = 68 
A Term = 
B Term = 

N = 69 
A Term = 
B Term = 

N = 70 
A Term = 
B Term = 

N = 71 
A Term = 
B Term = 

N = 72 
A Term = 
B Term = 

N = 73 
A Term = 
B Term = 

N = 74 
A Term = 
B Term = 

N = 75 
A Term = 
B Term = 

N = 76 
A Term = 
B Term = 

N = 77 
A Term = 
B Term = 



0.00047679155335931677*cos 
0.00041315228612778701*sin 



-0.00830793700369917950-cos 
-0 .01216080976658390900*si n 



0.00531478250010750540*cos 
0 .00691916775064072180*si n 



-0.01327873704499691600*cos 
-0.00358982664265015830*sin 



0.00470244489436939810*cos 
0.00499158341120479330*sin 



0.00044186078444939566*cos 
0.00754977211548748040*sin 



-0.00061916606725582958-cos 
-0.00474470156635318360*sin 



0.01107194058225292700*cos 
0.00682202118963723140*sin 



-0.00524047959257894010*cos 
-0.00931792666515206080*sin 



0.00452833197511323060*cos 
-0.00129668626795288990*sin 



-0.00414237595537471530*cos 
-0.00136299709646369310*sin 



-0.00638004680540780500*cos 
-0.00483632724329400870*sin 



0.00155576684012201500*cos 
0.00903414933907483940*sin 



-0.00525803316638073350*cos 
-0.00227114132641837710*sin 



O.00530363606356375540*cos 
0.00760031213570345680*sin 



0.00294788000475155200*cos 
0.00025259568186932000*sin 



coil wavel. out 



3 . 1415927*62*x/10 . 00000000) 
3 . 1415927*62*x/10 . 00000000) 



3 . 1415927*63*x/10 . 00000000) 
3 . 1415927*63*x/10 . 00000000) 



3 . 1415927*64*x/10 . 00000000) 
3 . 1415927*64*x/10 . 00000000) 



3 . 1415927*65*x/10 . 00000000) 
3 . 1415927*65*x/10 . 00000000) 



3 . 1415927*66*x/10 . 00000000) 
3 . 1415927*66*x/10 . 00000000) 



3 . 1415927*67*x/10 . 00000000) 
3 . 1415927*67*x/10 . 00000000) 



3 . 1415927*68*x/10 . 00000000) 
3 . 1415927*68*x/10 . 00000000) 



3 . 1415927*69*x/10 . 00000000) 
3 . 1415927*69*x/10 . 00000000) 



3 . 1415927*70*x/10 . 00000000) 
3 . 1415927*70*x/10 . 00000000) 



3 . 1415927*71*x/10 . 00000000) 
3 . 1415927*71*x/10 . 00000000) 



3 . 1415927*72*x/10 . 00000000) 
3 . 1415927*72*x/10 . 00000000) 



3. 1415927*73*x/10. 00000000) 
3 . 1415927*73*x/10 . 00000000) 



3 . 1415927*74*x/10 . 00000000) 
3 . 1415927*74*x/10 . 00000000) 



3 . 1415927*75*x/10 . 00000000) 
3 . 1415927*75*x/10 . 00000000) 



3. 1415927*76*x/10. 00000000) 
3. 1415927*76*x/10. 00000000) 



3 . 1415927*77*x/10 . 00000000) 
3 . 1415927*77*x/10 . 00000000) 
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N = 78 
A Term = 
B Term = 

N = 79 
A Term = 
B Term = 

N = 80 
A Term = 
B Term = 

N = 81 
A Term = 
B Term = 

N = 82 
A Term = 
B Term = 

N = 83 
A Term = 
B Term = 

N = 84 
A Term = 
B Term = 

N = 85 
A Term = 
B Term = 

N = 86 
A Term = 
B Term = 

N = 87 
A Term = 
B Term = 

N = 88 
A Term = 
B Term = 

N = 89 
A Term = 
B Term = 

N = 90 
A Term = 
B Term = 

N = 91 
A Term = 
B Term = 

N = 92 
A Term = 
B Term = 

N = 93 
A Term = 
B Term = 



coil wavel. out 



0.00280764231574142300*cos 
-0 .00390452186095862710*si n 



0.00369946527672188700*cos 
0 . 00120461488497710730*51 n 



-0.00262932464685880370*cos 
-0.O0995793363575586160*sin 



-0.00235225704843293250*cos 
0 . 00280121879824070950*si n 



-0.00443582836621307310*cos 
-0.00288959821728245350*sin 



-0.00290887390425370550*cos 
0 .00287886478308926360*si n 



-0.00109739835456087440*cos 
0.00717063047961887660*sin 



0.00323820353726577990*cos 
-0.00176020640344507190*sin 



0.00268203541840628200*cos 
0.00718653169857251960*sin 



0.00434887683760564960*cos 
-0.00626093030519325660*sin 



0.00276227379572644180*cos 
-0.00146935962970417940*sin 



-0.00298820460291454070*cos 
-0.00266707707982515550*sin 



0.00025423899620244172*cos 
-0.00684243892835677620*sin 



-0.00685319685956094820*cos 
0.00598253495912291770*sin 



-0.00148730332746402350*cos 
-0.00326519497112327070*sin 



0.00021463015728170366*cos 
0.00728853904814733870*sin 



3 . 1415927*78*x/10 . 00000000) 
3 . 1415927*78*x/10 . 00000000) 



3. 1415927*79*x/10. 00000000) 
3. 1415927*79*x/10. 00000000) 



3 . 1415927*80*x/10 . 00000000) 
3 . 1415927*80*x/10 . 00000000) 



3 . 1415927*81*x/10 . 00000000) 
3 . 1415927*81*x/10 . 00000000) 



3 . 1415927*82*x/10 . 00000000) 
3 . 1415927*82*x/10 . 00000000) 



3 . 1415927*83*x/10 . 00000000) 
3 . 1415927*83*x/10 . 00000000) 



3 . 1415927*84*x/10 . 00000000) 
3 . 1415927*84*x/10 . 00000000) 

3. 1415927*85*x/10. 00000000) 
3. 1415927*85*x/10. 00000000) 



3 . 1415927*86*x/10 .00000000) 
3. 1415927*86*x/10. 00000000) 



3 . 1415927*87*x/10 . 00000000) 
3 . 1415927*87*x/10 . 00000000) 



3 . 1415927*88*x/10 . 00000000) 
3. 1415927*88*x/10. 00000000) 



3. 1415927*89*x/10. 00000000) 
3. 1415927*89*x/10. 00000000) 



3 . 1415927*90*x/10 . 00000000) 
3 . 1415927*90*x/10 . 00000000) 



3 . 1415927*91*x/10 . 00000000) 
3. 1415927*91*x/10. 00000000) 



3 . 1415927*92*x/10 . 00000000) 
3 . 1415927*92*x/10 . 00000000) 



3. 1415927*93*x/10. 00000000) 
3. 1415927*93*x/10. 00000000) 
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coil wavel. out 



N = 94 

A Term = -0.00153328816747250450*cos(3.1415927*94*x/10. 00000000) 

B Term = 0.00302744005210805630*sin(3.1415927*94*x/10. 00000000) 

N = 95 

A Term = 0.00789054521888109960*cos(3.1415927*95*x/10. 00000000) 

B Term = -0.00190828057130283600*sin(3 .1415927*95*x/10. 00000000) 

N = 96 

A Term = -0.00093151629604490036*cos(3.1415927*96*x/10. 00000000) 

B Term = 0.00437374052172606520*sin(3.1415927*96*x/10. 00000000) 

N = 97 

A Term = 0.00401213500758983420*cos(3.1415927*97*x/10. 00000000) 

B Term = -0.00893521242518463540*sin(3.1415927*97*x/10. 00000000) 

N = 98 

A Term = 0.00014458920017085305*cos(3.1415927*98*x/10. 00000000) 

B Term = 0.00098812893815039180*sin(3.1415927*98*x/10. 00000000) 

N = 99 

A Term = -0.00588939443032514100*cos(3.1415927*99*x/10. 00000000) 

B Term = -0.00343719833602320820*sin(3.1415927*99*x/10. 00000000) 

N = 100 

A Term = 0.00183194952704563850*cos(3.1415927*100*x/10. 00000000) 

B Term = -0.00210907083750795560*sin(3.1415927*100*x/10. 00000000) 
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coilwave2.out 

Fourier Series of Coil #2 waveForm 



N = 0 
A Term = 
B Term = 

N = 1 
A Term = 
B Term = 

N = 2 
A Term = 
B Term = 

N = 3 
A Term = 
B Term = 

N = 4 
A Term = 
B Term = 

N = 5 
A Term = 
B Term = 

N = 6 
A Term = 
B Term = 

N = 7 
A Term = 
B Term = 

N = 8 
A Term = 
B Term = 

N = 9 
A Term = 
B Term = 

N = 10 
A Term = 
B Term = 

N = 11 
A Term = 
B Term = 

N = 12 
A Term = 
B Term = 

N = 13 
A Term = 
B Term = 

N = 14 
A Term = 
B Term = 



-2 . 74361999999999990000*cos 
0 . 00000000000000000000*si n 



-1. 15811931552203710000*cos 
-4 . 13845525093240150000*si n 



1. 32342418551513720000*cos 
-1.25340345275166750000*sin 



0 . 43062339326826637000*cos 
-0 . 66035027064241847000*si n 



0 . 84711439711922676000*cos 
-0 . 28985906421062146000*si n 



0 . 30784754656924740000*cos 
0 . 12045508076403209000*si n 



0. 31598393736477626000*cos 
-0.03663733665335258700*sin 



0 . 1690802 5881034028000*cos 
0 . 222O5448731130051OOO*si n 



-0.02500733130668391900*cos 
0 . 02227928441091642200* si n 



0 . 09859030544601543600*cos 
0.04376589727745757600*sin 



-0.07858550079168165300*cos 
0 . 05042126348264409900*si n 



0.03923496228642797000*cos 
-0 . 10986526678221016000*si n 



0 . 04169530091430402600*cos 
0 . 05410387764194173700*si n 



-0 . 02184716823760496400*cos 
-0.09358682604005536OOO*sin 



0 . 12973417978276905000*cos 
0 . 00780322487170587770*si n 



3 . 1415927*0*x/10 . 00000000) 
3 . 1415927*0*x/10 . 00000000) 



3 . 1415927*l*x/10 . 00000000) 
3 . 1415927*l*x/10 . 00000000) 



3 . 1415927*2*x/10 . 00000000) 
3 . 1415927*2*x/10 . 00000000) 



3 . 1415927*3*x/10 . 00000000) 
3 . 1415927*3*x/10 .00000000) 



3 . 1415927*4*x/10 . 00000000) 
3 . 1415927*4*x/10 . 00000000) 



3. 1415927*5*x/10. 00000000) 
3. 1415927*5*x/10. 00000000) 



3. 1415927*6*x/10. 00000000) 
3 . 1415927*6*x/10 . 00000000) 



3 . 1415927*7*x/10 . 00000000) 
3 . 1415927*7*x/10 . 00000000) 



3 . 1415927*8*x/10 . 00000000) 
3 . 1415927*8*x/10 . 00000000) 



3. 1415927*9*x/10. 00000000) 
3 . 1415927*9* x/10 . 00000000) 



3 . 1415927*10*x/10 . 00000000) 
3 . 1415927*10*x/10 . 00000000) 



3 . 1415927*ll*x/10 . 00000000) 
3 . 1415927*ll*x/10 . 00000000) 



3 . 1415927*12*x/10 . 00000000) 
3 . 1415927*12*x/10 . 00000000) 



3 . 1415927*13*x/10 . 00000000) 
3 . 1415927*13*x/10 .00000000) 



3 . 1415927*14*x/10 . 00000000) 
3 . 1415927*14*x/10 . 00000000) 
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coil wave2 .out 

N = 15 

A Term = -0.03512903837810938500*cos(3 . 1415927*15*x/10 .00000000) 

B Term = 0.02377366768572087100*sin(3 . 1415927*15*x/10 .00000000) 

N = 16 

A Term = 0.08461697643085192000*cos (3. 1415927*16*x/10. 00000000) 

B Term = -0.04802885001383761000*sin(3.1415927*16*x/10. 00000000) 

N = 17 

A Term = 0.01346508119979164900*cos (3. 1415927*17*x/10 . 00000000) 

B Term = 0.08577663069815681700*sin(3.1415927*17*x/10. 00000000) 

N = 18 

A Term = -0.01924792143832004800*cos(3.1415927*18*x/10. 00000000) 

B Term = -0.04664597881991877600*sin(3.1415927*18*x/10. 00000000) 

N = 19 

A Term = 0.06062152266883882900*cos(3.1415927*19*x/10. 00000000) 

B Term = 0.03498715537731213900*sin(3.1415927*19*x/10. 00000000) 

N = 20 

A Term = -0.05564215108891779100*cos(3 . 1415927*20*x/10 .00000000) 

B Term = 0.00903635704816874680*sin(3.1415927*20*x/10. 00000000) 

N = 21 

A Term = O.04559523293367111900*cos(3.1415927*21*x/10. 00000000) 

B Term = -O.04473185345555226500*sin(3.1415927*21*x/10. 00000000) 

N = 22 

A Term = -0. 00207830510075500020*cos (3 . 1415927*22*x/10 .00000000) 

B Term = 0.04735349259883375000*sin(3.1415927*22*x/10. 00000000) 

N = 23 

A Term = -0.00958332265630864630*cos (3. 1415927*23*x/10 . 00000000) 

B Term = -0.05474109917508092300*sin (3. 1415927*23*x/10 . 00000000) 

N = 24 

A Term = 0.05571730548432476900*cos(3 . 1415927*24*x/10 .00000000) 

B Term = 0.02465410998371071700*sin(3 . 1415927*24*x/10 .00000000) 

N = 25 

A Term = -O.03524028832314136900*cos(3 . 1415927*25*x/10 .00000000) 

B Term = 0.00007207656450507609*sin(3 . 1415927*25*x/10 .00000000) 

N = 26 

A Term = 0.04789741788620281700*cos(3.1415927*26*x/10. 00000000) 

B Term = -O.02211360821528025900*sin(3.1415927*26*x/10. 00000000) 

N = 27 

A Term = -0.00808618491691309900*cos (3. 1415927*27*x/10. 00000000) 

B Term = 0.04064425795521726300*sin(3.1415927*27*x/10. 00000000) 

N = 28 

A Term = -0. 00129488993250056680*cos (3 . 1415927*28*x/10 .00000000) 

B Term = -0.03361868854387191400*sin(3 . 1415927*28*x/10 .00000000) 

N = 29 

A Term = 0.02783762211641897100*cos(3 . 1415927*29*x/10 .00000000) 

B Term = 0.02493656147269140300*sin(3 . 1415927*29*x/10 .00000000) 

N = 30 

A Term = -0.02621844534886263000*cos (3. 1415927*30*x/10. 00000000) 

B Term = -0.00515772548244608500*sin(3.1415927*30*x/10. 00000000) 
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f 



coi lwave2 .out 



N = 31 
A Term = 
B Term = 

N = 32 
A Term = 
B Term = 

N = 33 
A Term = 
B Term = 

N = 34 
A Term = 
B Term = 

N = 35 
A Term = 
B Term = 

N = 36 
A Term = 
B Term = 

N = 37 
A Term = 
B Term = 

N = 38 
A Term = 
B Term = 

N = 39 
A Term = 
B Term = 

N = 40 
A Term = 
B Term = 

N = 41 
A Term = 
B Term = 

N = 42 
A Term = 
B Term = 

N = 43 
A Term = 
B Term = 

N = 44 
A Term = 
B Term = 

N = 45 
A Term = 
B Term = 

N = 46 
A Term = 



0.02891337156763972000*cos 
-0.01268678848227504200*sin 



-0 . 00830920443493120800*cos 
0.02043594662800166100*sin 



0.00385996174731023050*cos 
-0.02416761808845794800*sin 



0.01720029205073521700*cos 
0.01665958564798683200*sin 



-0.01233775033100043000*cos 
-0.00762781601253642930*sin 



0.02075354033223669700*cos 
-0.00134898158507033630*sin 



-0 . 00790484178793027410*cos 
0.00828933254464486740*sin 



0.00923848689191959090*cos 
-0 . 01025909079018400500*si n 



0 . 00264053096314467890*cos 
0 . 01005807152201909600*si n 



0 . 00028011886422642002*cos 
-0.00799445845787278600*sin 



0 . 00814800814730738 500*cos 
0.00566753357585613770*sin 



-0.00351369993873072890*cos 
-0 . 00305828484174879100*si n 



0 . 01011493619404944300*cos 
O.00069723301150333844*sin 



-0 . 00518921397483862110*cos 
0 . 00310564167662667630*si n 



0 . 00832740996511603490* cos 
-0.00682123892233490750*sin 



3 . 1415927*31*x/10 . 00000000) 
3 . 1415927*31*x/10 . 00000000) 



3 . 1415927*32*x/10 . 00000000) 
3 . 1415927*32*x/10 . 00000000) 



3 . 1415927*33*x/10 . 00000000) 
3 . 1415927*33*x/10 . 00000000) 



3 . 1415927*34*x/10 . 00000000) 
3 . 1415927*34*x/10 . 00000000) 



3 . 1415927*35*x/10 . 00000000) 
3. 1415927*35*x/10. 00000000) 



3 . 1415927*36*x/10 . 00000000) 
3 . 1415927*36*x/10 . 00000000) 



3 . 1415927*37*x/10 . 00000000) 
3 . 1415927*37*x/10 . 00000000) 



3 . 1415927*38*x/10 .00000000) 
3 . 1415927*38*x/10 .00000000) 



3 . 1415927*39*x/10 . 00000000) 
3 . 1415927*39*x/10 .00000000) 



3 . 1415927*40*x/10 . 00000000) 
3 . 1415927*40*x/10 . 00000000) 



3 . 1415927*41*x/10 . 00000000) 
3 . 1415927*41*x/10 . 00000000) 



3 . 1415927*42*x/10 . 00000000) 
3 . 1415927*42*x/10 . 00000000) 



3 . 1415927*43*x/10 . 00000000) 
3 . 1415927*43*x/10 . 00000000) 



3 . 1415927*44*x/10 . 00000000) 
3 . 1415927*44*x/10 . 00000000) 



3 . 1415927*45*x/10 .00000000) 
3 . 1415927*45*x/10 .00000000) 



-0 . 00032768809059656960*cos(3 . 1415927*46*x/10 .00000000) 
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B Term = 



coilwave2 .out 

0.00988727722111274430*sin(3.1415927*46*x/10. 00000000) 



M 


— A 7 


A 


Term = 


B 


Term = 


N 


A Q 


A 


Term = 


B 


Term = 


N 


= 4y 


A 


Term = 


B 


Term = 


N 


50 


A 


Term = 


B 


Term = 


N 


51 


A 


Term = 


B 


Term = 


N 


52 


A 


Term = 


B 


Term = 


N 


53 


A 


Term = 


B 


Term = 


N 


54 


A 


Term = 


B 


Term = 


N 


55 


A 


Term = 


B 


Term = 


N 


56 


A 


Term = 


B 


Term = 


N 


57 


A 


Term = 


B 


Term = 


N 


r o 

= 58 


A 


Term = 


B 


Term = 


N 


= 59 


A 


Term = 


B 


Term = 


N 


= 60 


A 


Term = 


B 


Term = 


N 


- 61 


A 


Term = 


B 


Term = 


N 


= 62 



-0.00097582794472149360*cos 
-O.01159966053475196300*sin 



O.01158983616738575700*cos 
0.00857826534050223390*sin 



-O.01110373118138918000*cos 
-0.00346001676823489870*sin 



0 . 0165 5 569436061189000*cos 
-0.00464252736236808300*sin 



-0.00751168404419990700*cos 
0.01222777075120149400*sin 



0.00378205799511278920*cos 
-0.01547588193192155400*sin 



0.00926775425763013210*cos 
0.01496482252288258200*sin 



-0.01340001816532341400*cos 
-0.00763775351273100500*sin 



O.01830221973188598700*cos 
-0.00196562842521091690*sin 



-0.01189718600162438300*cos 
0.01131937258747469000*sin 



0.00586830870420980330*cos 
-0.01773761656012484000*sin 



0.00783255050350688440*cos 
0 . 01672097O828546483O0*si n 



-0.01300548059275832100-cos 
-O.OllO234670O15O3215O0*sin 



0.02001989560323314000*cos 
0.00075102584276270283*sin 



-0.01370761922602530000*cos 
0.00959766986324264890*sin 



3 . 1415927*47*x/10 . 00000000) 
3 . 1415927*47*x/10 . 00000000) 



3 . 1415927*48*x/10 . 00000000) 
3 . 1415927*48*x/10 . 00000000) 



3 . 1415927*49*x/10 . 00000000) 
3 . 1415927*49*x/10 . 00000000) 



3 . 1415927* 50*x/10 . 00000000) 
3 . 1415927*50*x/10 . 00000000) 



3 . 1415927* 51*x/10 . 00000000) 
3 . 1415927*51*x/10 . 00000000) 



3 . 1415927* 52*x/10 . 00000000) 
3 . 1415927* 52*x/10 . 00000000) 



3 . 141592 7* 53*x/10 . 00000000) 
3. 1415927* 53*x/10. 00000000) 



3 . 1415927* 54*x/10 . 00000000) 
3 . 1415927* 54*x/10 . 00000000) 



3. 1415927*55*x/10. 00000000) 
3 . 1415927*5 5*x/10 . 00000000) 



3 . 1415927* 56*x/10 . 00000000) 
3. 1415927*56*x/10. 00000000) 



3 . 1415927*57*x/10 . 00000000) 
3 . 1415927* 57*x/10 . 00000000) 



3 . 1415927* 58*x/10 . 00000000) 
3 . 1415927*58*x/10 . 00000000) 



3 . 1415927*59*x/10 . 00000000) 
3 . 1415927* 59*x/10 . 00000000) 



3. 1415927*60*x/10. 00000000) 
3 . 1415927*60*x/10 . 00000000) 



3 . 1415927*61*x/10 . 00000000) 
3 . 1415927*61*x/10 . 00000000) 
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coilwave2 .out 

A Term = 0.00875217772440555450*cos(3.1415927*62*x/10. 00000000) 

B Term = -O.015700738456653O44OO*sin(3.1415927*62*x/lO. 00000000) 

N = 63 

A Term = 0.00406007566968070830*cos(3 .1415927*63*x/10. 00000000) 

B Term = 0.01712474329421512200*sin(3 .1415927*63*x/10. 00000000) 

N = 64 

A Term = -0.01045890113075900700*cos(3 .1415927*64*x/10. 00000000) 

B Term = -0.01169095657029163600*sin(3 .1415927*64*x/10. 00000000) 

N = 65 

A Term = 0.01657614104308649200*cos(3 .1415927*65*x/10. 00000000) 

B Term = 0.00312610880966039280*sin(3 .1415927*65*x/10. 00000000) 

N = 66 

A Term = -0.01314054353460840700*cos(3 .1415927*66*x/10. 00000000) 

B Term = 0.00592605293524422690*sin(3.1415927*66*x/10. 00000000) 

N = 67 

A Term = 0.00903150642165032910*cos(3 .1415927*67*x/10. 00000000) 

B Term = -0.01272730332340424100*sin(3 .1415927*67*x/10. 00000000) 

N = 68 

A Term = 0.00154723980059945430*cos(3 .1415927*68*x/10. 00000000) 

B Term = 0.01408033400703412100*sin(3.1415927*68*x/10. 00000000) 

N = 69 

A Term = -0.00660192133034991080*cos(3 .1415927*69*x/10. 00000000) 

B Term = -0.01130733450086722300*sin(3 .1415927*69*x/10. 00000000) 

N = 70 

A Term = 0.01323494767243209100*cos(3 .1415927*70*x/10. 00000000) 

B Term = 0.00488025411003456620*sin(3 .1415927*70*x/10. 00000000) 

N = 71 

A Term = -0.01086117802871475900*cos(3 . 1415927*71*x/10 . 00000000) 

B Term = 0.00229653403162496030*sin(3 .1415927*71*x/10. 00000000) 

N = 72 

A Term = 0 . 00937328935657385200*cos(3 . 1415927*72*x/10 .00000000) 

B Term = -0.00757446122364574540-sin(3.1415927*72*x/10. 00000000) 

N = 73 

A Term = -0.00183264048480829940*cos(3.1415927*73*x/10. 00000000) 

B Term = 0.01017542459615867500*sin(3.1415927*73*x/10. 00000000) 

N = 74 

A Term = -0.00182193996849948920*cos(3.1415927*74*x/10. 00000000) 

B Term = -0.00911792243181569480*sin(3.1415927*74*x/10. 00000000) 

N = 75 

A Term = 0.00730073558510019660*cos(3.1415927*75*x/10. 00000000) 

B Term = 0.00585174229159036390*sin(3 .1415927*75*x/10. 00000000) 

N = 76 

A Term = -0.00695045763453506700*cos(3.1415927*76*x/10. 00000000) 

B Term = -0.00158513372314571870*sin(3.1415927*76*x/10. 00000000) 

N = 77 

A Term = 0.00777216245233860490*cos(3.1415927*77*x/10. 00000000) 

B Term = -0.00240737547212706680*sin(3.1415927*77-x/10. 00000000) 
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N = 78 
A Term = 
B Term = 

N = 79 
A Term = 
B Term = 

N = 80 
A Term = 
B Term = 

N = 81 
A Term = 
B Term = 

N = 82 
A Term = 
B Term = 

N = 83 
A Term = 
B Term = 

N = 84 
A Term = 
B Term = 

N = 85 
A Term = 
B Term = 

N = 86 
A Term = 
B Term = 

N = 87 
A Term = 
B Term = 

N = 88 
A Term = 
B Term = 

N = 89 
A Term = 
B Term = 

N = 90 
A Term = 
B Term = 

N = 91 
A Term = 
B Term = 

N = 92 
A Term = 
B Term = 

N = 93 
A Term = 
B Term = 



coi lwave2 .out 



-0 . 00367980693096718080*cos 
0.00502551597939557630*sin 



0.00239031393832929300*cos 
-0.00621821150778771590*sin 



0.00215980368637591970*cos 
0 . 00600002604147404410*si n 



-O.00297933191507433360*cos 
-0.00472677784687316940*sin 



0.00619979051524979180*cos 
0.00262558276476868810*sin 



-0.00540773386751254670*cos 
-0 . 00004720366443422604*si n 



0.00647911453000702030*cos 
-0.00286642983112907460*sin 



-0.00312095836221843370*cos 
0.00539640498238052900*sin 



0.00156110297820006500*cos 
-0.00688292614293857690*sin 



O.00362582333822096750*cos 
0.00675761097804106200*sin 



-0.00563368694505811050*cos 
-0.00435817347312019830*sin 



0.00867727448905274840*cos 
0.00047053983625271748*sin 



-0.00684310309148810930*cos 
0 . 00420312461844598360*si n 



0.00471027591097354910*cos 
-0.00798498583283479070*sin 



0.00147686294893744340*cos 
0.00925537900037681250*sin 



-O.00543605435333257250*cos 
-O.00770577504671658350*sin 



3. 1415927*78*x/10. 00000000) 
3 . 1415927*78' v x/10 .00000000) 



3. 1415927*79*x/10. 00000000) 
3. 1415927*79*x/10. 00000000) 



3. 1415927*80*x/10. 00000000) 
3 . 1415927*80*x/10 . 00000000) 



3 . 1415927*81*x/10 . 00000000) 
3 . 1415927*81*x/10 . 00000000) 



3. 1415927*82*x/10. 00000000) 
3 . 1415927*82*x/10 . 00000000) 



3. 1415927*83*x/10. 00000000) 
3. 1415927*83*x/10. 00000000) 



3. 1415927*84*x/10. 00000000) 
3 . 1415927*84*x/10 . 00000000) 



3 . 1415927*85*x/10 . 00000000) 
3 . 1415927*85*x/10 . 00000000) 



3 . 1415927*86*x/10 . 00000000) 
3 . 1415927*86*x/10 . 00000000) 



3 . 1415927*87*x/10 .00000000) 
3 . 1415927*87*x/10 . 00000000) 



3 . 1415927*88*x/10 . 00000000) 
3 . 1415927*88*x/10. 00000000) 



3. 1415927*89*x/10. 00000000) 
3. 1415927*89*x/10. 00000000) 



3. 1415927*90*x/10. 00000000) 
3 . 1415927*90*x/10 . 00000000) 



3 . 1415927*91*x/10 . 00000000) 
3 . 1415927*91*x/10 . 00000000) 



3. 1415927*92*x/10. 00000000) 
3 . 1415927*92*x/10 . 00000000) 



3 . 1415927*93*x/10 .00000000) 
3 . 1415927*93*x/10 .00000000) 
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coilwave2 .out 

A Term = 0.00875217772440555450*cos(3 .1415927*62*x/10. 00000000) 

B Term = -0.015700738456653044OO*sin(3 .1415927*62*x/10. 00000000) 

N = 63 

A Term = 0.00406007566968070830*cos(3 .1415927*63*x/10. 00000000) 

B Term = 0.01712474329421512200*sin(3 .1415927*63*x/10. 00000000) 

N = 64 

A Term = -0. 01045890113075900700*cos (3 .1415927*64*x/10. 00000000) 

B Term = -O.O1169095657O291636OO*sin(3 .1415927-64*x/10. 00000000) 

N = 65 

A Term = 0.01657614104308649200*cos(3 .1415927*65*x/10. 00000000) 

B Term = 0.00312610880966039280*sin(3 .1415927*65*x/10. 00000000) 

N = 66 

A Term = -0.01314054353460840700*cos(3 .1415927*66*x/10. 00000000) 

B Term = 0.00592605293524422690*si n(3 .1415927*66*x/10 .00000000) 

N = 67 

A Term = 0.00903150642165032910*cos(3 .1415927*67*x/10. 00000000) 

B Term = -0.01272730332340424100*si n(3 .1415927*67*x/10 .00000000) 

N = 68 

A Term = 0.00154723980059945430*cos(3 .1415927*68*x/10. 00000000) 

B Term = 0.01408033400703412100*si n(3 . 1415927*68*x/10 .00000000) 

N = 69 

A Term = -0. 00660192133034991080*cos (3 .1415927*69*x/10. 00000000) 

B Term = -0.01130733450086722300*sin(3 .1415927*69*x/10. 00000000) 

N = 70 

A Term = 0.01323494767243209100*cos(3 .1415927*70*x/10. 00000000) 

B Term = 0.00488025411003456620*sin(3.1415927*70*x/10. 00000000) 

N = 71 

A Term = -0.01086117802871475900*cos(3 .1415927*71*x/10. 00000000) 

B Term = O.0O229653403162496O3O*sin(3.1415927*71*x/lO. 00000000) 

N = 72 

A Term = 0.00937328935657385200*cos(3 .1415927*72*x/10. 00000000) 

B Term = -0.00757446122364574540*si n(3 . 1415927*72*x/10 . 00000000) 

N = 73 

A Term = -0. 00183264048480829940*cos (3 .1415927*73*x/10. 00000000) 

B Term = 0.01017542459615867500*sin(3.1415927*73*x/10. 00000000) 

N = 74 

A Term = -0.00182193996849948920*cos(3 .1415927*74*x/10. 00000000) 

B Term = -0.00911792243181569480*sin(3.1415927*74*x/10. 00000000) 

N = 75 

A Term = 0.00730073558510019660*cos(3 . 1415927*75*x/10 . 00000000) 

B Term = 0.00585174229159036390*sin(3 .1415927*75*x/10. 00000000) 

N = 76 

A Term = -0.00695045763453506700*cos (3. 1415927*76*x/10. 00000000) 

B Term = -0.00158513372314571870*sin (3. 1415927*76*x/10 . 00000000) 

N = 77 

A Term = 0.00777216245233860490*cos(3 .1415927*77*x/10. 00000000) 

B Term = -0.00240737547212706680*sin(3.1415927*77*x/10. 00000000) 
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N = 78 
A Term = 
B Term = 

N = 79 
A Term = 
B Term = 

N = 80 
A Term = 
B Term = 

N = 81 
A Term = 
B Term = 

N = 82 
A Term = 
B Term = 

N = 83 
A Term = 
B Term = 

N = 84 
A Term = 
B Term = 

N = 85 
A Term = 
B Term = 

N = 86 
A Term = 
B Term = 

N = 87 
A Term = 
B Term = 

N = 88 
A Term = 
B Term = 

N = 89 
A Term = 
B Term = 

N = 90 
A Term = 
B Term = 

N = 91 
A Term = 
B Term = 

N = 92 
A Term = 
B Term = 

N = 93 
A Term = 
B Term = 
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-0.00367980693096718080*cos 
0.0O50255159793955763O*sin 



O.00239031393832929300*cos 
-0.00621821150778771590*sin 



0.00215980368637591970*cos 
0 . 00600002604147404410*si n 



-0.00297933191507433360*cos 
-0.00472677784687316940*sin 



0.00619979051524979180*cos 
0.00262558276476868810*sin 



-0.00540773386751254670*cos 
-0 . 00004720366443422604*si n 



0.00647911453000702030*cos 
-0.00286642983112907460*sin 



-0.00312095836221843370*cos 
0.00539640498238052900*sin 



0 . 00156110297820006500*cos 
-0 . 00688292614293857690*si n 



0.00362582333822096750*cos 
0 .00675761097804106200*si n 



-0.00563368694505811050*cos 
-O.0043581734731201983O*sin 



0.00867727448905274840*cos 
0.00047053983625271748*sin 



-0 . 00684310309148810930*cos 
0.00420312461844598360*sin 



0.00471027591097354910*cos 
-0.00798498583283479070*sin 



0.00147686294893744340*cos 
0. 009255379000376812 50*sin 



-0.00543605435333257250*cos 
-0.00770577504671658350*sin 



3 . 1415927*78*x/10 . 00000000) 
3 . 141592 7*78*x/10 . 00000000) 



3 . 1415927*79*x/10 . 00000000) 
3 . 1415927*79*x/10 . 00000000) 



3 . 1415927*80*x/10 . 00000000) 
3 . 1415927*80*x/10 . 00000000) 



3 . 1415927*81*x/10 . 00000000) 
3 . 1415927*81*x/10 . 00000000) 



3 . 1415927*82*x/10 . 00000000) 
3 . 1415927*82*x/10 . 00000000) 



3 . 1415927*83*x/10 . 00000000) 
3 . 1415927*83*x/10 . 00000000) 



3. 1415927*84*x/10. 00000000) 
3 . 1415927*84*x/10 . 00000000) 



3. 1415927*85*x/10. 00000000) 
3 . 1415927*85*x/10 . 00000000) 



3 . 1415927*86*x/10 . 00000000) 
3. 1415927*86*x/10. 00000000) 



3 . 1415927*87*x/10 . 00000000) 
3 . 1415927*87*x/10 . 00000000) 



3 . 1415927*88*x/10 . 00000000) 
3 . 141592 7*88*x/10 . 00000000) 



3. 1415927*89*x/10. 00000000) 
3 . 1415927*89*x/10 . 00000000) 



3 . 1415927*90*x/10 . 00000000) 
3 . 1415927*90*x/10 . 00000000) 



3 . 1415927*91*x/10 . 00000000) 
3 . 1415927*91*x/10 . 00000000) 



3 . 1415927*92*x/10 . 00000000) 
3 . 1415927*92*x/10 . 00000000) 



3 . 1415927*93*x/10 . 00000000) 
3 . 1415927*93*x/10 . 00000000) 
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N = 94 

A Term = 0. 01018870208156160900*cos (3 .1415927*94*x/10. 00000000) 

B Term = 0.00311540145693248000*sin(3 .1415927*94*x/10. 00000000) 

N = 95 

A Term = -0.00924508141264390390*cos(3.1415927*95*x/10. 00000000) 

B Term = 0.00264373533982534100*sin(3.1415927*95*x/10. 00000000) 

N = 96 

A Term = 0.00751081839151258610*cos(3 .1415927*96*x/10. 00000000) 

B Term = -0.00772018440263961150*sin(3.1415927*96*x/10. 00000000) 

N = 97 

A Term = -0.00081043351747648196*cos(3 .1415927*97*x/10. 00000000) 

B Term = 0.01040496000703253400*sin(3 .1415927*97*x/10. 00000000) 

N = 98 

A Term = -0.00411922464695734220*cos(3 .1415927*98-x/10. 00000000) 

B Term = -0.00923823644255666360*sin(3 .1415927*98*x/10. 00000000) 

N = 99 

A Term = 0.00951601802805959910*cos(3 .1415927*99"'x/10. 00000000) 

B Term = 0.00508904349955530130*sin(3 .1415927*99*x/10. 00000000) 

N = 100 

A Term = -0.00979379124793574620*cos(3 .-1415927*100*x/10. 00000000) 

B Term = 0.00096861240696119878*sin(3.1415927*100*x/10. 00000000) 
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Fourier Series of Coil #3 waveForm 



coilwave3 .out 
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-7.28779999999999360000*cos 
0 . 00000000000000000000*si n 



4 . 40715604299113120000*cos 
9 . 48810106687465190000*si n 



0 . 4006479832 5 502660000*cos 
-3 . 28761517848058430000*si n 



0.77450352575835657000*cos 
3.00597579732173470000*sin 



0. 32484976266351878000*cos 
-1.80561445393920030000*sin 



0. 39783212224486947000*cos 
1.64962599323958160000*sin 



0.20531877078256025000*cos 
-1. 34707502388659780000*sin 



0.18186716825823435000*cos 
1.08922973093976600000*sin 



O.08752115102527345500*cos 
-1.06572834815167330000*sin 



O.03555203586400892300*cos 
0.84112058592427630000*sin 



0.02534657041743482200*cos 
-0.81923635123851324000*sin 



-0.01091289054795027900*cos 
0.72346945839718757000*sin 



0 . 02636481902202903600*cos 
-0.63392873932575311000*sin 



O.01330319972389277600*cos 
0 . 63147432919845603000*si n 



0.04703532035650099400*cos 
-0. 54116281376694564000*sin 



3 . 1415927*0*x/10 . 00000000) 
3 . 1415927*0*x/10 . 00000000) 



3 . 1415927*l*x/10 . 00000000) 
3 . 1415927*l*x/10 . 00000000) 



3 . 1415927*2*x/10 . 00000000) 
3 . 1415927*2*x/10 . 00000000) 



3 . 1415927*3*x/10 . 00000000) 
3 . 1415927*3*x/10 . 00000000) 



3 . 1415927*4*x/10 . 00000000) 
3 . 1415927*4*x/10 . 00000000) 



3 . 1415927*5*x/10 . 00000000) 
3 . 1415927*5*x/10 . 00000000) 



3 . 1415927*6*x/10 . 00000000) 
3 . 1415927*6*x/10 . 00000000) 



3 . 1415927*7*x/10 . 00000000) 
3 . 1415927*7*x/10 . 00000000) 



3 . 1415927*8*x/10 . 00000000) 
3 . 1415927*8*x/10 . 00000000) 



3 . 1415927*9*x/10 . 00000000) 
3 . 1415927*9*x/10 . 00000000) 



3 . 1415927*10*x/10 . 00000000) 
3 . 1415927*10*x/10 . 00000000) 



3 . 1415927*ll*x/10 . 00000000) 
3 . 1415927*ll*x/10 . 00000000) 



3 . 1415927*12*x/10 . 00000000) 
3 . 1415927*12*x/10 . 00000000) 



3 . 1415927*13*x/10 . 00000000) 
3 . 1415927*13*x/10 . 00000000) 



3. 1415927*14*x/10. 00000000) 
3. 1415927*14*x/10. 00000000) 
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Fourier Series of Coil #3 waveForm 



N = 0 
A Term = 
B Term = 
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A Term = 
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B Term = 

N = 8 
A Term = 
B Term = 
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B Term = 

N = 11 

A Term = 
B Term = 

N = 12 
A Term = 
B Term = 

N = 13 

A Term = 
B Term = 

N = 14 
A Term = 
B Term = 



-7 . 28779999999999360000*cos 
0 . 00000000000000000000*si n 



4.40715604299113120000*cos 
9 . 48810106687465190000*si n 



0.40064798325502660000*cos 
-3.28761517848058430000*sin 



O.77450352575835657000*cos 
3.O05975797321734700O0*sin 



0 . 32484976266351878000*cos 
-1.80561445393920030000*sin 



0 . 39783212224486947000*cos 
1.64962599323958160000*sin 



0 . 20531877078256025000*cos 
-1.34707502388659780000*sin 



0 . 18186716825823435000*cos 
1.08922973093976600000*sin 



0.08752115102527345500*cos 
-1.06572834815167330000*sin 



0.03555203586400892300*cos 
0.84112058592427630000*sin 



0.02534657041743482200*cos 
-0.819236351238513240OO*sin 



-0.01091289054795027900*cos 
0.72346945839718757000*sin 



0.02636481902202903600*cos 
-0.63392873932575311000*sin 



0.01330319972389277600*cos 
0 . 63147432919845603000*si n 



0.04703532035650099400*cos 
-0. 54116281376694564000*sin 



3 . 1415927*0*x/10 . 00000000) 
3 . 1415927*0*x/10 . 00000000) 



3 . 1415927*l*x/10 .00000000) 
3 . 1415927*l*x/10 .00000000) 



3 . 1415927*2*x/10 . 00000000) 
3 . 1415927*2*x/10 . 00000000) 



3 . 1415927*3*x/10 . 00000000) 
3. 1415927*3*x/10. 00000000) 



3 . 1415927*4*x/10 .00000000) 
3 . 1415927*4*x/10 .00000000) 



3 . 1415927*5*x/10 .00000000) 
3. 1415927*5*x/10. 00000000) 



3. 1415927*6*x/10. 00000000) 
3. 1415927*6*x/10. 00000000) 



3. 1415927*7*x/10. 00000000) 
3. 1415927*7*x/10. 00000000) 



3. 1415927*8*x/10. 00000000) 
3 . 1415927*8*x/10 . 00000000) 



3 . 1415927*9*x/10 . 00000000) 
3 . 1415927*9*x/10 . 00000000) 



3 . 1415927*10*x/10 . 00000000) 
3 . 1415927*10*x/10 . 00000000) 



3 . 1415927*ll*x/10 . 00000000) 
3 . 1415927*ll*x/10 . 00000000) 



3. 1415927*12*x/10. 00000000) 
3. 1415927*12*x/10. 00000000) 



3 . 1415927*13*x/10 . 00000000) 
3 . 1415927*13*x/10 . 00000000) 



3 . 1415927*14*x/10 . 00000000) 
3 . 1415927*14*x/10 . 00000000) 
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N = 15 
A Term = 
B Term = 

N = 16 
A Term = 
B Term = 

N = 17 
A Term = 
B Term = 

N = 18 
A Term = 
B Term = 

N = 19 

A Term = 
B Term = 

N = 20 
A Term = 
B Term = 
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A Term = 
B Term = 

N = 22 
A Term = 
B Term = 

N = 23 
A Term = 
B Term = 

N = 24 
A Term = 
B Term = 

N = 25 
A Term = 
B Term = 

N = 26 
A Term = 
B Term = 

N = 27 
A Term = 
B Term = 

N = 28 
A Term = 
B Term = 

N = 29 
A Term = 
B Term = 

N = 30 
A Term = 
B Term = 



coilwave3.out 



0.03315260086104406700*cos 
0. 53645547242637026000*sin 



0.04764553358501690000*cos 
-0 . 50472880968666611000*si n 



O.01253876304330971100*cos 
0 . 46020092691204306000*si n 



0.03277522824397354600*cos 
-0.46263615105982886000*sin 



-0 . 01848225989021693900*cos 
0 . 415 56045242696493000*si n 



0.02772433785500234800*COS 
-0.39970362985068530000*sin 



-0.02245154894023093300*cos 
0.38362520579149861000*sin 



0.03493578679148168900*cos 
-0.347452725473860520OO*sin 



-0.00693736131994262400*cos 
0 . 34566490907174585000*si n 



0 . 03549393014968986500*cos 
-0.32501707269493757000*sin 



-0.00286567885084546870*cos 
0 . 30993451922964416000*si n 



0.02462208561631050100*cos 
-0 . 31316689012176907000*si n 



-0.01559926452876011900*cos 
0 . 29161036281315478000*si n 



0 . 01986687708992120000*cos 
-0 . 28946182460697789000*si n 



-0.02406050516438000600*cos 
0 . 28289878068399082000*si n 



0.03036884706583470000*cos 
-0.26024659365835329000*sin 



3 . 1415927*15*x/10 . 00000000) 
3 . 1415927*15*x/10 . 00000000) 



3 . 1415927*16*x/10 . 00000000) 
3 . 1415927*16*x/10 . 00000000) 



3 . 1415927*17*x/10 . 00000000) 
3. 1415927*17*x/10. 00000000) 



3 . 1415927*18*x/10 . 00000000) 
3 . 1415927*18*x/10 . 00000000) 



3 . 1415927*19*x/10 . 00000000) 
3 . 1415927*19*x/10 . 00000000) 



3 . 1415927*20*x/10 . 00000000) 
3 . 1415927*20*x/10 . 00000000) 



3 . 1415927*21*x/10 . 00000000) 
3 . 1415927*21*x/10 . 00000000) 



3 . 1415927*22*x/10 . 00000000) 
3 . 1415927*22*x/10 . 00000000) 



3 . 1415927*23*x/10 . 00000000) 
3 . 1415927*23*x/10 . 00000000) 



3 . 1415927*24*x/10 . 00000000) 
3 . 1415927*24*x/10 . 00000000) 



3 . 1415927*2 5*x/10 . 00000000) 
3. 1415927*2 5*x/10. 00000000) 



3 . 1415927*26*x/10 . 00000000) 
3. 1415927*26*x/10. 00000000) 



3 . 1415927*27*x/10 . 00000000) 
3 . 1415927*27*x/10 . 00000000) 



3 . 1415927*28*x/10 . 00000000) 
3 . 1415927*28*x/10 . 00000000) 



3 . 1415927*29*x/10 . 00000000) 
3 . 1415927*29*x/10 . 00000000) 



3 . 1415927*30*x/10 . 00000000) 
3 . 1415927*30*x/10 . 00000000) 
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N = 31 

A Term = -0.01861038861489003000*cos(3 . 1415927*31*x/10 .00000000) 

B Term = 0.26351778956237620000*sin(3.1415927*31*x/10. 00000000) 

N = 32 

A Term = 0. 03979764826020104600*cos (3 .1415927*32*x/10. 00000000) 

B Term = -0.24302355529934427000*sin(3.1415927*32*x/10. 00000000) 

N = 33 

A Term = -0.01314889167406609700*cos(3 .1415927*33*x/10. 00000000) 

B Term = 0. 23359274865663179000*si n(3 . 1415927*33*x/10. 00000000) 

N = 34 

A Term = 0.03178514108667273400*cos(3 .1415927*34*x/10. 00000000) 

B Term = -0.23594123689237462000*sin(3.1415927*34*x/10. 00000000) 

N = 35 

A Term = -0.01817635647281173900*cos(3 .1415927*35*x/10. 00000000) 

B Term = 0. 21436679739124920000*sin(3 . 1415927*35*x/10. 00000000) 

N = 36 

A Term = 0.01683894791163258200*cos(3.1415927*36*x/10. 00000000) 

B Term = -0. 22404347849280248000*sin(3 . 1415927*36*x/10 .00000000) 

N = 37 

A Term = -0.02432656797033769600*cos(3 . 1415927*37*x/10 .00000000) 

B Term = 0.21455932244368531000*sin(3.1415927*37*x/10. 00000000) 

N = 38 

A Term = 0.01701341263167277200*cos(3.1415927*38*x/10. 00000000) 

B Term = -0.20588250499755475000*sin(3.1415927*38*x/10. 00000000) 

N = 39 

A Term = -0.02128695440828249800*cos(3.1415927*39*x/10. 00000000) 

B Term = 0.21564363262638464000*sin(3.1415927*39*x/10. 00000000) 

N = 40 

A Term = 0.03196903389535456200*cos(3.1415927*40*x/10. 00000000) 

B Term = -0.19386669082772631000*sin(3.1415927*40*x/10. 00000000) 

N = 41 

A Term = -0.01602882410045956500*cos(3.1415927*41*x/10. 00000000) 

B Term = 0.19946430930028430000*sin(3.1415927*41*x/10. 00000000) 

N = 42 

A Term = O.03936593535861420600*cos(3 .1415927*42*x/10. 00000000) 

B Term = -0.19077828096657160000*sin(3 . 1415927*42*x/10 . 00000000) 

N = 43 

A Term = -0.01944122325976594900*cos(3 .1415927*43*x/10. 00000000) 

B Term = 0.17526659214039758000*sin(3 .1415927*43*x/10. 00000000) 

N = 44 

A Term = 0.02858769507756694600*cos(3 . 1415927*44*x/10 . 00000000) 

B Term = -0.18468586642086063000*sin(3.1415927*44*x/10. 00000000) 

N = 45 

A Term = -0.02653671621534021300*cos(3 . 1415927*45*x/10 . 00000000) 

B Term = 0.16495786139304283000*sin(3 . 1415927*45*x/10 . 00000000) 



N = 46 
A Term = 



0.01582863836890429500*cos(3.1415927*46*x/10. 00000000) 
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= 51 


A 


Te rm = 


B 


Term = 


M 

IN 


52 


A 
/A 


Tp r m — 
1 CI III — 


B 


Term = 


N 


— j j 




Te rm = 


B 


Term = 


IN 


54 


A 
rA 


Tp rm — 
1 C 1 1 II — 


B 


Term = 


N 


= 55 


A, 


Tp rm — 

1 C 1 1 II — 


B 


Term = 


N 


56 


A 


Term = 


B 


Term = 


N 


57 


A 


Te rm — 


B 


Term = 


N 


58 


A 


Term = 


B 


Term = 


N 


59 


A 


Term = 


B 


Term = 


N 


60 


A 


Term = 


B 


Term = 


N 


= 61 


A 


Term = 


B 


Term = 


N 


= 62 



coil wave3 .out 

-0 . 17021777617144418000*si n (3 . 1415927*46*x/10 . 00000000) 



-0.02507814195665177000*cos 
0.17032042726322175000*sin 



0.01892457967576626600*cos 
-0 . 15837348021016978000*si n 



-0.01731665576742365300*cos 
0 . 17180061639874977000*si n 



0.03168666971098311300*cos 
-0.15756088993692205000*sin 



-0.01689606063328836700*cos 
0 . 15826736268999339000*si n 



0.03523298321666368200*cos 
-0 . 15840714248111584000*si n 



-0.02564433517396170300*cos 
0.14159985639404982000*sin 



O.02610452458390985200*cos 
-O.14883089350222131000*sin 



-0.02968481085162888100*cos 
0 . 13679117789105730000*si n 



0.01874072113117799000*cos 
-0 . 13451847846681267000*si n 



-0.02220960163398961000*cos 
0.14059500669783626000*sin 



O.02241508354580401500*cos 
-0 . 13014309409883673000*si n 



-O.01526855385520465400*cos 
0.13934107854961181000*sin 



0.02954065693639712200*cos 
-0.13517023525401087000*sin 



-0 . 02049410402033369200*cos 
0.12981081875725270000*sin 



3 . 1415927*47*x/10 . 00000000) 
3 . 1415927*47*x/10 . 00000000) 



3 . 1415927*48*x/10 . 00000000) 
3 . 1415927*48*x/10 . 00000000) 



3 . 1415927*49*x/10 . 00000000) 
3 . 1415927*49*x/10 . 00000000) 



3 . 1415927*50*x/10 . 00000000) 
3 . 1415927*50*x/10 . 00000000) 



3 . 1415927*51*x/10 . 00000000) 
3 . 1415927*51*x/10 . 00000000) 



3 . 1415927*52*x/10 . 00000000) 
3 . 1415927*52*x/10 . 00000000) 



3 . 1415927*53*x/10 . 00000000) 
3 . 1415927*53*x/10 . 00000000) 



3 . 1415927*54*x/10 . 00000000) 
3 . 1415927*54*x/10 . 00000000) 



3 . 1415927*55*x/10 • 00000000) 
3 . 1415927*55*x/10 . 00000000) 



3 . 1415927*56*x/10 ■ 00000000) 
3 . 1415927*56*x/10 . 00000000) 



3 . 1415927*57*x/10 . 00000000) 
3 . 1415927*57*x/10 . 00000000) 



3 . 1415927*58*x/10 . 00000000) 
3. 1415927*58*x/10. 00000000) 



3 . 1415927*59*x/10 . 00000000) 
3 . 1415927*59*x/10 . 00000000) 



3 . 1415927*60*x/10 . 00000000) 
3 . 1415927*60*x/10 . 00000000) 



3 . 1415927*61*x/10 . 00000000) 
3 . 1415927*61*x/10 . 00000000) 



Page 4 



coi 1wave3 .out 



A 


Te rm ~ 


B 


Term = 


N 
in 


63 


A 
M 


1 CI III — 


B 


Term = 


in 


64 


A 
M 


1 C 1 ill — 


B 


Term = 


M 
IN 


— O J 


A 
A 


1 CI III — 


B 


Term = 


M 
IN 


— DO 


A 

A 


Te rm = 


B 


Term = 


M 

IN 


= 0/ 


A 

A 


Te rm = 


B 


Term = 


M 


= DO 


A 

A 


Te rm = 


B 


Term = 


M 

N 




A 

A 


Te rm = 


B 


Term = 


M 

N 




A 

A 


Te rm = 


B 


Term = 


M 


— 71 

— / ± 


A 

A 


Te rm = 


B 


Term = 


M 
IN 


— I L 


A 
A 


*T* r\ f* m 

1 tr 1 III — 


B 


Term = 


M 
IN 


— 73 

— / J 


A 


1 c I III — 


B 


Term = 


M 
IN 


- 74 

— / 1 


A 
M 


To r>m — 
1 Cllll — 


B 


Term = 


N 

IN 


- 7^ 


A 


Tp r»m — 
1 C 1 III — 


R 


To rrn 

I c I III — 


N 


= 76 


A 


Term = 


B 


Term = 


N 


77 


A 


Term = 


B 


Term = 



0.02985287719262436800*cos 
-0 . 13416070163878444000*si n 



-0.02998469094031682900*cos 
0 . 12106134987919698000*si n 



0.02508532871450993100*cos 
-0 . 12158606129414624000*si n 



-0.02877478945390350000*cos 
0 . 11869366220116620000*si n 



0.02313216990965811700*cos 
-0 . 11088446870413478000*si n 



-0.01889701778896906200*cos 
0 . 11819920798825204000*si n 



0.02515445061853360500*cos 
-0 . 1129445 52878517460OO*si n 



-0 . 01614811779447110000*cos 
0.11468592350343210000*sin 



0.02637699103210623000*cos 
-0 . 11872137102811642000*si n 



-0.02488044892808013700*cos 
0 . 110067 5 5700208100000*si n 



0.02554426551445442000*cos 
-0.11378307790927356000*sin 



-0.03132886187627011800*cos 
0.10734375040054522000*sin 



0.02564479365057211300*cos 
-0 . 10111760445300262000*si n 



-0.02569017606163896000*cos 
0 . 10478722481254854000*si n 



0.02682840404423353400*cos 
-0.09590217657033799100*sin 



-0 . 01729471778039086800*cos 
0.10037251087629728000*sin 



3 . 1415927*62*x/10 . 00000000) 
3 . 1415927*62*x/10 . 00000000) 



3 . 1415927*63*x/10 . 00000000) 
3 . 1415927*63*x/10 . 00000000) 



3 . 1415927*64*x/10 . 00000000) 
3 . 1415927*64*x/10 . 00000000) 



3 . 1415927*65*x/10 . 00000000) 
3 . 1415927*65*x/10 . 00000000) 



3 . 1415927*66*x/10 . 00000000) 
3 . 1415927*66*x/10 . 00000000) 



3 . 1415927*67*x/10 . 00000000) 
3 . 1415927*67*x/10 . 00000000) 



3 . 1415927*68*x/10 . 00000000) 
3 . 1415927*68*x/10 . 00000000) 



3 . 1415927*69*x/10 . 00000000) 
3 . 1415927*69*x/10 . 00000000) 



3 . 1415927*70*x/10 . 00000000) 
3. 1415927*70*x/10. 00000000) 



3 . 1415927*71*x/10 . 00000000) 
3 . 1415927*71*x/10 . 00000000) 



3 . 1415927*72*x/10 . 00000000) 
3. 1415927*72*x/10. 00000000) 



3 . 1415927*73*x/10 . 00000000) 
3 . 1415927*73*x/10 . 00000000) 



3 . 1415927*74*x/10 . 00000000) 
3 . 141592 7*74*x/10 . 00000000) 



3 . 1415927*75*x/10 . 00000000) 
3 . 1415927*75*x/10 . 00000000) 



3 . 1415927*76*x/10 . 00000000) 
3 . 1415927*76*x/10 . 00000000) 



3 . 1415927*77*x/10 . 00000000) 
3 . 1415927*77*x/10. 00000000) 
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N = 78 
A Term = 
B Term = 

N = 79 
A Term = 
B Term = 

N = 80 
A Term = 
B Term = 

N = 81 

A Term = 
B Term = 

N = 82 
A Term = 
B Term = 

N = 83 
A Term = 
B Term = 

N = 84 
A Term = 
B Term = 

N = 85 
A Term = 
B Term = 

N = 86 
A Term = 
B Term = 

N = 87 
A Term = 
B Term = 

N = 88 
A Term = 
B Term = 

N = 89 
A Term = 
B Term = 

N = 90 
A Term = 
B Term = 

N = 91 
A Term = 
B Term = 

N = 92 
A Term = 
B Term = 

N = 93 
A Term = 
B Term = 



coil wave 3. out 



0.02596773565153813500*cos 
-0 . 10136390731208866000*si n 



-0 . 01927676217266644600*cos 
0.09674972831371569900*sin 



0.02345057400101445100*cos 
-O.10432649278214973O00*sin 



-0.02805288764499059400*cos 
0.09655183639764985900*sin 



0.02366406369962355300*cos 
-0.09646619154288227400*sin 



-0.02986045804636907100*cos 
0.09654888694749914900*sin 



O.02714215572224735300*cos 
-0 . 08678827765174096400- si n 



-0 . 02260418071865823200*cos 
0.09232260580069852300*sin 



0.02846258372300060300*cos 
-0.08652436790794175400*sin 



-0.01815197280754313200*cos 
0 .08629715498591289400" si n 



0.02494542714246842200*cos 
-0.09176173262592352000*sin 



-0.02299306047023928400*cos 
O.O846388338352068030O*sin 



0.02183268427359664300*cos 
-0 . 09079144193510409200*si n 



-0.02886934404226929700*cos 
0.08704235678518203200*sin 



0.02425112497408456200*cos 
-0.08271268263635517300*sin 



-0.02690069944200416300*cos 
0.08655412262207529100*sin 



3. 1415927*78*x/10. 00000000) 
3. 1415927*78*x/10. 00000000) 



3 . 1415927*79*x/10 . 00000000) 
3 . 1415927*79*x/10 . 00000000) 



3. 1415927*80*x/10. 00000000) 
3. 1415927*80*x/10. 00000000) 



3 . 1415927*81*x/10 . 00000000) 
3 . 1415927*81*x/10 . 00000000) 



3. 1415927*82*x/10. 00000000) 
3. 1415927*82*x/10. 00000000) 



3. 1415927*83*x/10. 00000000) 
3. 1415927*83*x/10. 00000000) 



3 . 1415927*84*x/10 . 00000000) 
3 . 1415927*84*x/10 . 00000000) 



3. 1415927*85*x/10. 00000000) 
3 . 1415927*85*x/10 . 00000000) 



3 . 1415927*86*x/10 . 00000000) 
3 . 1415927*86*x/10 . 00000000) 



3 . 1415927*87*x/10 . 00000000) 
3 . 1415927*87- x/10 . 00000000) 



3. 1415927*88*x/10. 00000000) 
3. 1415927*88*x/10. 00000000) 



3. 1415927*89*x/10. 00000000) 
3. 1415927*89*x/10. 00000000) 



3. 1415927*90*x/10. 00000000) 
3. 1415927*90*x/10. 00000000) 



3. 1415927*91*x/10. 00000000) 
3. 1415927*91*x/10. 00000000) 



3. 1415927*92*x/10. 00000000) 
3. 1415927*92*x/10. 00000000) 



3 . 1415927*93*x/10 . 00000000) 
3 . 1415927*93*x/10 . 00000000) 
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coil wave 3 .out 

N = 94 

A Term = 0.02856105970443219300*cos(3.1415927*94*x/10. 00000000) 

B Term = -0.07759640326448347100*sin(3.1415927*94*x/10. 00000000) 

N = 95 

A Term = -0.02115368829517351400*cos(3.1415927*95*x/10. 00000000) 
B Term = 0.08067735808281140100*sin(3.1415927*95*x/10. 00000000) 

N = 96 

A Term = 0.02782388973497394000*cos(3.1415927*96*x/10. 00000000) 

B Term = -0.07985090460762139800*sin(3.1415927*96*x/10. 00000000) 

N = 97 

A Term = -0.02083767080448749300*cos(3.1415927*97*x/10. 00000000) 
B Term = 0.07575058316382737900*sin(3.1415927*97*x/10. 00000000) 

N = 98 

A Term = 0.02308497371616000100*cos(3 .1415927*98*x/10. 00000000) 

B Term = -0.08225049738435669100*sin(3.1415927*98*x/10. 00000000) 

N = 99 

A Term = -0.02559749606225375200*cos(3.1415927*99*x/10. 00000000) 
B Term = 0.07687992528559334200*sin(3.1415927*99*x/10. 00000000) 

N = 100 

A Term = 0.02193505074687203800*cos(3 .1415927*100*x/10. 00000000) 

B Term = -0.07856685406744359800*sin(3 .1415927*100*x/10. 00000000) 
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Fourier Series of Coil #4 wave Form 



coilwave4.out 



N = 0 
A Term = 
B Term = 

N = 1 
A Term = 
B Term = 

N = 2 
A Term = 
B Term = 

N - 3 
A Term = 
B Term = 

N = 4 
A Term = 
B Term = 

N = 5 
A Term = 
B Term = 

N = 6 
A Term = 
B Term = 

N = 7 

A Term = 
B Term = 

N = 8 
A Term = 
B Term = 

N = 9 
A Term = 
B Term = 

N = 10 
A Term = 
B Term = 

N = 11 
A Term = 
B Term = 

N = 12 
A Term = 
B Term = 

N = 13 
A Term = 
B Term = 

N = 14 
A Term = 
B Term = 



0. 56617249999999986000*cos 
0 . 00000000000000000000*si n 



1.09249254712821610000*cos 
-2 . 31789864232462420000*si n 



0.00981179087095048790*cos 
-0.04112489142801858300*sin 



-0 . 77748985612891275000*cos 
-0.31261244856835335000*sin 



-O.21395885460663036000*cos 
-0.49731561287915355000*sin 



-0.06326466246765757900*cos 
-0.01474582744236930400*sin 



-0 . 14793156423871748000*cos 
0 . 14046831528945039000*si n 



-0.24797168551359183000*cos 
0 . 10920604048146498000*si n 



-0 . 10529937270414663000*cos 
-0.09912804285685128400*sin 



0.17167131777537348000*cos 
0 . 04169883644591979000*si n 



-0 . 03604389302164717100*cos 
0 . 16030589009013940000*si n 



-0.05687946556362513700*cos 
-0 . 04157379940630168100*si n 



0 . 00746347100707194300*cos 
-0 . 07741448217908968000*si n 



0.04677965051303101900*cos 
-0.05168827145398219300*sin 



0.01881935686424155600*cos 
0.06393327808509162400*sin 



3 . 1415927*0*x/10 . 00000000) 
3. 1415927*0*x/10. 00000000) 



3 . 1415927*l*x/10 . 00000000) 
3 . 1415927*l*x/10 . 00000000) 



3 . 1415927*2*x/10 . 00000000) 
3. 1415927*2*x/10. 00000000) 



3 . 1415927*3*x/10 . 00000000) 
3 . 1415927*3*x/10 . 00000000) 



3 . 1415927*4*x/10 . 00000000) 
3 . 1415927*4*x/10 . 00000000) 



3 . 1415927*5*x/10 .00000000) 
3 . 1415927*5*x/10 .00000000) 



3 . 1415927*6*x/10 . 00000000) 
3 . 1415927*6*x/10 . 00000000) 



3 . 1415927*7*x/10 .00000000) 
3 . 1415927*7*x/10 .00000000) 



3. 1415927*8*x/10. 00000000) 
3 . 1415927*8*x/10 . 00000000) 



3 . 1415927*9*x/10 . 00000000) 
3 . 1415927*9*x/10 . 00000000) 



3 . 1415927*10*x/10 . 00000000) 
3 . 1415927*10*x/10 . 00000000) 



3 . 1415927*ll*x/10 . 00000000) 
3 . 1415927*ll*x/10 . 00000000) 



3 . 1415927*12*x/10 . 00000000) 
3 . 1415927*12*x/10 . 00000000) 



3 . 1415927*13*x/10 . 00000000) 
3 . 1415927*13*x/10 . 00000000) 



3. 1415927*14*x/10. 00000000) 
3. 1415927*14*x/10. 00000000) 
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I 



coil wave4. out 

N = 15 

A Term = -0.14926251694484238000*cos(3 .1415927*15*x/10. 00000000) 

B Term = 0.00430833019779906530*sin(3 .1415927*15*x/10. 00000000) 

N = 16 

A Term = 0.00400934689146752090*cos (3. 1415927*16*x/10. 00000000) 

B Term = -0.09107947879529254600*sin(3 .1415927*16*x/10. 00000000) 

N = 17 

A Term = O.04107199162903536100*cos(3 . 1415927*17*x/10. 00000000) 

B Term = 0. 044805 50109050864700*sin (3. 1415927*17*x/10. 00000000) 

N = 18 

A Term = -0. 02513483010903990200*cos (3 . 1415927*18*x/10 .00000000) 

B Term = O.O5328141397O335104OO*sin (3. 1415927*18*x/10 . 00000000) 

N = 19 

A Term = -O.02319750824937588300*cos(3 .1415927*19*x/10. 00000000) 

B Term = 0.01486304287795250400*sin(3 . 1415927*19*x/10. 00000000) 

N = 20 

A Term = -0. 02981677853568691700*cos (3 . 1415927*20*x/10. 00000000) 

B Term = -0.04877271975207420100*sin (3. 1415927*20*x/10. 00000000) 

N = 21 

A Term = 0. 08234647111292993400*cos (3 .1415927*21*x/10. 00000000) 

B Term = -0.02280237125775157500*si n(3 . 1415927*21*x/10. 00000000) 

N = 22 

A Term = -0.01755909709160662800*cos(3 .1415927*22*x/10. 00000000) 

B Term = 0.06396359476175053600*sin(3 .1415927*22*x/10. 00000000) 

N = 23 

A Term = -0.06736754162658392700*cos(3 . 1415927*23*x/10. 00000000) 

B Term = -0.02724784418945111100*si n(3 . 1415927*23*x/10. 00000000) 

N = 24 

A Term = 0.01342347891048435700*cos(3 . 1415927*24*x/10 . 00000000) 

B Term = -0.03836727371815250300*sin(3 . 1415927*24*x/10 . 00000000) 

N = 25 

A Term = 0.00649006668336018790*cos(3 . 1415927*25*x/10 .00000000) 

B Term = 0.01044439965961951000*sin(3 . 1415927*25*x/10 . 00000000) 

N = 26 

A Term = 0.01623500452876516400*cos(3 . 1415927*26*x/10 .00000000) 

B Term = 0. 03832471896222652800*sin (3 .1415927*26*x/10. 00000000) 

N = 27 

A Term = -0.05164227255180037200*cos(3 . 1415927*27*x/10. 00000000) 

B Term = 0.01833808584766311000*sin(3 .1415927*27*x/10. 00000000) 

N = 28 

A Term = 0.00453933015923544110*cos(3 . 1415927*28*x/10 . 00000000) 

B Term = -0.05186192618569127100*sin(3 .1415927*28*x/10. 00000000) 

N = 29 

A Term = 0.05528616183023992600*cos(3 .1415927*29*x/10. 00000000) 

B Term = 0.00540162460879087780*sin (3. 1415927*29*x/10. 00000000) 

N = 30 

A Term = -0. 02322276098659497400*cos (3 . 1415927*30*x/10. 00000000) 

B Term = 0. 03022148830030148600*sin (3 .1415927*30*x/10. 00000000) 
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coilwave4.out 

N = 31 

A Term = -0.02443069062556931700*cos(3.1415927*31*x/10. 00000000) 

B Term = -0 .01644699884497586800*sin(3 . 1415927*31*x/10 .00000000) 

N = 32 

A Term = -0 .01821613913237197400*cos(3 . 1415927*32*x/10 .00000000) 

B Term = -0.02863727439156048600*sin(3.1415927*32*x/10. 00000000) 

N = 33 

A Term = 0 .02587845850361290100*cos(3 . 1415927*33*x/10 .00000000) 

B Term = -0.00182897415157960690*sin(3.1415927*33*x/10. 00000000) 

N = 34 

A Term = -0.00380630386689837390*cos(3 . 1415927*34*x/10 .00000000) 

B Term = 0.04401568361827439000*sin(3.1415927*34*x/10. 00000000) 

N = 35 

A Term = -0.04296442464423103100*cos (3. 1415927*35*x/10 . 00000000) 

B Term = 0.00518168838519577830*sin(3.1415927*35*x/10. 00000000) 

N = 36 

A Term = 0.01947501514769441700*cos(3.1415927*36*x/10. 00000000) 

B Term = -0.0277016680279423S200*sin(3.1415927*36*x/10. 00000000) 

N = 37 

A Term = O.02715232505063041100*cos(3.1415927*37*x/10. 00000000) 

B Term = 0.00664869595287078880*sin(3.1415927*37*x/10. 00000000) 

N = 38 

A Term = 0.00772989692749737490*cos(3 . 1415927*38*x/10 .00000000) 

B Term = 0.02140801458595822600*sin(3.1415927*38*x/10. 00000000) 

N = 39 

A Term = -0.02910081959317530400*cos(3.1415927*39*x/10. 00000000) 

B Term = -0.01078718172161334900*sin(3 . 1415927*39*x/10 . 00000000) 

N = 40 

A Term = -0.00523878133819048740*cos(3 . 1415927*40*x/10 .00000000) 

B Term = -0.03573126874614390000*sin(3 .1415927*40*x/10. 00000000) 

N = 41 

A Term = 0.02360476292719618900*cos(3 . 1415927*41*x/10 .00000000) 

B Term = -0.00044265542481909019*sin(3 . 1415927*41*x/10 .00000000) 

N = 42 

A Term = -0.02080349760969470900*cos(3.1415927*42*x/10. 00000000) 

B Term = 0.02717379084896266000*sin(3.1415927*42*x/10. 00000000) 

N = 43 

A Term = -0.02483902897068750200*cos(3.1415927*43*x/10. 00000000) 

B Term = 0.00568638916591429100*sin(3.1415927*43*x/10. 00000000) 

N = 44 

A Term = -0.00414881756937270950*cos(3.1415927*44*x/10. 00000000) 

B Term = -0. 01870697070901052700*sin (3 . 1415927*44*x/10 .00000000) 

N = 45 

A Term = 0.03171448566130897300*cos(3.1415927*45*x/10. 00000000) 

B Term = 0.01053882640254597000*sin(3 . 1415927*45*x/10 .00000000) 

N = 46 

A Term = 0.00468919299202278940*cos (3. 1415927*46*x/10. 00000000) 
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6 Term = 



coil wave4. out 

0.02752936935293270500*sin(3.1415927*46*x/10. 00000000) 



N = 47 
A Term = 
B Term = 

N = 48 
A Term = 
B Term = 

N = 49 
A Term = 
B Term = 

N = 50 
A Term = 
B Term = 

N = 51 

A Term = 
B Term = 

N = 52 
A Term = 
B Term = 

N = 53 
A Term = 
B Term = 

N = 54 
A Term = 
B Term = 

N = 55 
A Term = 
B Term = 

N = 56 
A Term = 
B Term = 

N = 57 
A Term = 
B Term = 

N = 58 
A Term = 
B Term = 

N = 59 
A Term = 
B Term = 

N = 60 
A Term = 
B Term = 

N = 61 
A Term = 
B Term = 



-0.02272593657108745400*cos 
-O.O11122483244204255OO*sin 



0.01397743790271679200*cos 
-O.O24703504952082078O0*sin 



0.01021185128754429500*cos 
-0.00899956149044097400*sin 



-0.00294503367668482550*cos 
0 . 01949396818174421700*si n 



-0.03271457524882915400*cos 
-0 . 00089064520998476050*si n 



-0 . 00604495086443738140*cos 
-0 . 02244617615981126300*si n 



0 . 02708269626803621100*cos 
0.01774578110907269400*sin 



-0.01182899689958353900*cos 
0.01936693249441007200*sin 



-O.00515750086941284100*cos 
O.00213869906170930690*sin 



O.00300092187372739010*cos 
-0.02002084345616441500*sin 



0.02033055804109628800*cos 
-0.00781400247078987630*sin 



0.00233154693600249850*cos 
0 . 02200449020186051200*si n 



-0.03431615861781579800*cos 
-0.01476626517866104500*sin 



0.00793568465276016030*cos 
-0 . 01385648442370168500*si n 



0.00723145867476599620*cos 
0.00710424121692789390*sin 



3 . 1415927*47*x/10 . 00000000) 
3 . 1415927*47*x/10 . 00000000) 



3 . 1415927*48*x/10 . 00000000) 
3 . 1415927*48*x/10 . 00000000) 



3 . 1415927*49*x/10 . 00000000) 
3 . 1415927*49*x/10 . 00000000) 



3 . 1415927* 50*x/10 . 00000000) 
3 . 1415927*50*x/10 . 00000000) 



3 . 1415927*51*x/10 . 00000000) 
3 . 1415927*51*x/10 . 00000000) 



3 . 1415927* 52*x/10 . 00000000) 
3 . 1415927*52*x/10 . 00000000) 



3 . 1415927* 53*x/10 . 00000000) 
3 . 1415927*53*x/10 . 00000000) 



3 . 1415927* 54*x/10 . 00000000) 
3 . 1415927* 54*x/10 . 00000000) 



3. 1415927*55*x/10. 00000000) 
3 . 1415927*5 5*x/10 . 00000000) 



3 . 1415927*56*x/10 . 00000000) 
3 . 1415927*56*x/10 . 00000000) 



3 . 1415927*57*x/10 .00000000) 
3. 1415927*57*x/10. 00000000) 



3 . 1415927* 58*x/10 . 00000000) 
3 . 1415927* 58*x/10 . 00000000) 



3 . 1415927* 59*x/10 . 00000000) 
3 . 1415927* 59*x/10 . 00000000) 



3 . 1415927*60*x/10 . 00000000) 
3 . 1415927*60*x/10 . 00000000) 



3 . 1415927*61*x/10 . 00000000) 
3 . 1415927*61*x/10 . 00000000) 



N = 



62 
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coilwave4.out 



A 


Term = 


B 


Term = 


N 


= 63 


A 


Term = 


B 


Term = 


N 


= 64 


A 


Term = 


B 


Term = 


N 


= 65 


A 


Term = 


B 


Term = 


N 


= 66 


A 


Term = 


B 


Term = 


N 


67 


A 


Term = 


B 


Term = 


N 


68 


A 


Term = 


B 


Term = 


N 


69 


A 


Term = 


B 


Term = 


N 


70 


A 


Term = 


B 


Term = 


N 


71 


A 


Term = 


B 


Term = 


N 


72 


A 


Term = 


B 


Term = 


N 


73 


A 


Term = 


B 


Term = 


N 


74 


A 


Term = 


B 


Term = 


N 


= 75 


A 


Term = 


B 


Term = 


N 


76 


A 


Term = 


B 


Term = 


N 


77 


A 


Term = 


B 


Term = 



-0 . 00456024206520187680*cos 
0.01717558278867930800*sin 



-0.00974452387384814290*cos 
0 . 00864698139981097060*si n 



-O.00245335710225589140*cos 
-0 . 02417407162882102800*si n 



0.02970187298312643900*cos 
0.00696823948358570640*sin 



-O.00857732452333697870*cos 
0 . 01171639683764747300*si n 



-0 . 01736246711015310500*cos 
-0 . 01083087685738750400*si n 



0 . 00342093869984540530*cos 
-0 . 01129418504591160600*si n 



0.00420758510165385790*cos 
-0 . 00337744789368669840*si n 



O.00124039556629589020*cos 
0 . 02497602651608023200*si n 



-0.02052629372881104700*cos 
-0.00214241763108399170*sin 



0.00702708696474984870*cos 
-0.01390357489938043100*sin 



O.02090529260106091500*cos 
0.00819013187553470700*sin 



-0.00539747728947537020*cos 
0.00572335102581345190*sin 



-0 . 00974414609535226640*cos 
-0 . 00113711028677424060*si n 



-O.00231795833719622860*cos 
-0.02196368422837568000*sin 



0.01080788202668424300*cos 
0 . 00317909796239368610*si n 



3 . 1415927*62*x/10 . 00000000) 
3 . 1415927*62*x/10 . 00000000) 



3 . 1415927*63*x/10 . 00000000) 
3 . 1415927*63*x/10 . 00000000) 



3 . 1415927*64*x/10 . 00000000) 
3 . 1415927*64*x/10 . 00000000) 



3 . 1415927*65*x/10 . 00000000) 
3 . 1415927*65*x/10 . 00000000) 



3 . 1415927*66*x/10 . 00000000) 
3 . 1415927*66*x/10 . 00000000) 



3 . 1415927*67*x/10 . 00000000) 
3 . 1415927*67*x/10 . 00000000) 



3 . 1415927*68*x/10 . 00000000) 
3 . 1415927*68*x/10 . 00000000) 



3 . 1415927*69*x/10 . 00000000) 
3 . 1415927*69*x/10 . 00000000) 



3 . 1415927*70*x/10 . 00000000) 
3 . 1415927*70*x/10 . 00000000) 



3 . 1415927*71*x/10 . 00000000) 
3 . 1415927*71*x/10 . 00000000) 



3 . 1415927*72*x/10 . 00000000) 
3 . 1415927*72*x/10 . 00000000) 



3 . 1415927*73*x/10 . 00000000) 
3 . 1415927*73*x/10 . 00000000) 



3 . 1415927*74*x/10 . 00000000) 
3 . 1415927*74*x/10 . 00000000) 



3 . 1415927*75*x/10 . 00000000) 
3 . 1415927*75*x/10 . 00000000) 



3 . 1415927*76*x/10 . 00000000) 
3 . 1415927*76*x/10 . 00000000) 



3 . 1415927*77*x/10 . 00000000) 
3 . 1415927*77*x/10 . 00000000) 
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B Term = 



coil wave4. out 

0.02752936935293270500*sin(3.1415927*46*x/10. 00000000) 



N = 47 
A Term = 
B Term = 

N = 48 
A Term = 
B Term = 

N = 49 
A Term = 
B Term = 

N = 50 
A Term = 
B Term = 

N = 51 
A Term = 
B Term = 

N = 52 
A Term = 
B Term = 

N = 53 
A Term = 
B Term = 

N = 54 
A Term = 
B Term = 

N = 55 
A Term = 
B Term = 

N = 56 
A Term = 
B Term = 

N = 57 
A Term = 
B Term = 

N = 58 
A Term = 
B Term = 

N = 59 
A Term = 
B Term = 

N = 60 
A Term = 
B Term = 

N = 61 

A Term = 
B Term = 



-0.02272593657108745400*cos 
-0.01112248324420425500*si n 



0.01397743790271679200*cos 
-O.O24703504952082O78OO*sin 



0.01021185128754429500*cos 
-0 . 00899956149044097400*si n 



-0.00294503367668482550*cos 
0 . 01949396818174421700 *si n 



-0.03271457524882915400*cos 
-0.00089064520998476050*sin 



-0.00604495086443738140*cos 
-0 . 02244617615981126300*si n 



0.02708269626803621100*cos 
0.01774578110907269400-si n 



-0.01182899689958353900*cos 
0.01936693249441007200*sin 



-0.00515750086941284100*cos 
0.00213869906170930690*sin 



0.00300092187372739010*cos 
-0.02002084345616441500*sin 



O.02033055804109628800*cos 
-0 . 00781400247078987630*si n 



0.00233154693600249850*cos 
0.02200449020186051200*sin 



-0.03431615861781579800*cos 
-0.01476626517866104500*sin 



0.00793568465276016030*cos 
-0 . 01385648442370168500*si n 



0.00723145867476599620*cos 
0.00710424121692789390*sin 



N = 



62 



3. 1415927*47*x/10. 00000000) 
3 . 1415927*47*x/10 . 00000000) 



3 . 1415927*48*x/10 . 00000000) 
3 . 1415927*48*x/10 . 00000000) 



3 . 1415927*49*x/10 . 00000000) 
3. 1415927*49*x/10. 00000000) 



3. 1415927*50*x/10. 00000000) 
3 . 1415927*50*x/10 • 00000000) 



3 . 1415927*51*x/10 . 00000000) 
3 . 1415927* 51*x/10 . 00000000) 



3 . 1415927*52*x/10 . 00000000) 
3 . 1415927*52*x/10 . 00000000) 



3 . 1415927*53*x/10 . 00000000) 
3 . 1415927*53*x/10 . 00000000) 



3 . 1415927*54*x/10 . 00000000) 
3 . 1415927*54*x/10 . 00000000) 



3. 1415927*55*x/10. 00000000) 
3 . 1415927*55*x/10 . 00000000) 



3. 1415927*56*x/10. 00000000) 
3. 1415927*56*x/10. 00000000) 

3 . 1415927*57*x/10 . 00000000) 
3 . 1415927*57*x/10 . 00000000) 

3 . 1415927*58*x/10 . 00000000) 
3 . 1415927*58*x/10 .00000000) 



3. 1415927*59*x/10. 00000000) 
3 . 1415927*59*x/10 . 00000000) 



3 . 1415927*60*x/10 . 00000000) 
3 . 1415927*60*x/10 . 00000000) 



3 . 1415927*61*x/10 . 00000000) 
3 . 1415927*61*x/10 . 00000000) 
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coilwave4.out 

A Term = -0.00456024206520187680*cos(3.1415927*62*x/10. 00000000) 

B Term = 0.01717558278867930800*sin(3.1415927*62*x/10. 00000000) 

N = 63 

A Term = -0.00974452387384814290*cos(3.1415927*63*x/10. 00000000) 

B Term = 0.00864698139981097060*sin(3.1415927*63*x/10. 00000000) 

N = 64 

A Term = -0.00245335710225589140*cos(3.1415927*64*x/10. 00000000) 

B Term = -0.02417407162882102800*sin(3.1415927*64*x/10. 00000000) 

N = 65 

A Term = 0.02970187298312643900*cos(3.1415927*65*x/10. 00000000) 

B Term = 0.00696823948358570640*sin(3.1415927*65*x/10. 00000000) 

N = 66 

A Term = -0.00857732452333697870*cos(3.1415927*66*x/10. 00000000) 

B Term = 0.01171639683764747300*sin(3.1415927*66*x/10. 00000000) 

N = 67 

A Term = -0.01736246711015310500*cos(3.1415927*67*x/10. 00000000) 

B Term = -0.01083087685738750400*sin(3.1415927*67*x/10. 00000000) 

N = 68 

A Term = 0.00342093869984540530*cos(3.1415927*68*x/10. 00000000) 

B Term = -0.01129418504591160600*sin(3.1415927*68*x/10. 00000000) 

N = 69 

A Term = 0 . 00420758510165385790*cos(3 . 1415927*69*x/10 . 00000000) 

B Term = -0.00337744789368669840*sin(3.1415927*69*x/10. 00000000) 

N = 70 

A Term = 0.00124039556629589020*cos(3.1415927*70*x/10. 00000000) 

B Term = 0.02497602651608023200*sin(3.1415927*70*x/10. 00000000) 

N = 71 

A Term = -0.02052629372881104700*cos(3.1415927*71*x/10. 00000000) 

B Term = -0.00214241763108399170*sin(3.1415927*71*x/10. 00000000) 

N = 72 

A Term = 0.00702708696474984870*cos(3.1415927*72*x/10. 00000000) 

B Term = -0.01390357489938043100*sin(3.1415927*72*x/10. 00000000) 

N = 73 

A Term = 0.02090529260106091500*cos(3.1415927*73*x/10. 00000000) 

B Term = 0.00819013187553470700*sin(3.1415927*73*x/10. 00000000) 

N = 74 

A Term = -0.00539747728947537020*cos(3.1415927*74*x/10. 00000000) 

B Term = 0.00572335102581345190*sin(3.1415927*74*x/10. 00000000) 

N = 75 

A Term = -0.00974414609535226640*cos(3.1415927*75*x/10. 00000000) 

B Term = -0.00113711028677424060*sin(3.1415927*75*x/10. 00000000) 

N = 76 

A Term = -0.00231795833719622860*cos(3.1415927*76*x/10. 00000000) 

B Term = -0.02196368422837568000*sin(3.1415927*76*x/10. 00000000) 

N = 77 

A Term = 0.01080788202668424300*cos(3.1415927*77*x/10. 00000000) 

B Term = 0.00317909796239368610*sin(3.1415927*77*x/10. 00000000) 
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coil wave4. out 

N = 78 

A Term = -0.00665140754206137910*cos(3.1415927*78*x/10. 00000000) 

B Term = 0.01791275789818169400*sin(3.1415927*78*x/10. 00000000) 

N = 79 

A Term = -0.01793034566028110400*cos(3.1415927*79*x/10. 00000000) 

B Term = -0.00461813302022995O5O*sin(3.1415927*79*x/lO. 00000000) 

N = 80 

A Term = 0.00564425421353095640*cos(3.1415927*80*x/10. 00000000) 

B Term = -0.00367366629853589510*sin(3 .1415927*80*x/10. 00000000) 

N = 81 

A Term = 0.01541832536973611300*cos(3 .1415927*81*x/10. 00000000) 

B Term = 0.00062323749650626179*sin(3.1415927*81*x/10. 00000000) 

N = 82 

A Term = 0.00063109435374041913*cos(3 .1415927*82*x/10. 00000000) 

B Term = 0.01618455773494796600*sin(3 .1415927*82*x/10. 00000000) 

N = 83 

A Term = -0.01058925718458241500*cos(3.1415927*83*x/10. 00000000) 

B Term = -0.00628860531959835220*sin(3.1415927*83*x/10. 00000000) 

N = 84 

A Term = 0. 00414014139447912640*cos (3 .1415927*84*x/10. 00000000) 

B Term = -0.02011236363676289000*sin(3.1415927*84*x/10. 00000000) 

N = 85 

A Term = 0.00917247192479116010*cos(3 .1415927*85*x/10. 00000000) 

B Term = 0.00518839745243783110*sin(3.1415927*85*x/10. 00000000) 

N = 86 

A Term = -0.00648514514502907330*cos(3 . 1415927*86*x/10 . 00000000) 

B Term = 0.00543494720765342730*sin(3.1415927*86*x/10. 00000000) 

N = 87 

A Term = -0.01718658620480913200*cos(3.1415927*87*x/10. 00000000) 

B Term = 0.00280033832777892710*sin(3.1415927*87*x/10. 00000000) 

N = 88 

A Term = 0.00034537575385590845*cos(3.1415927*88*x/10. 00000000) 

B Term = -0.0108982588071767O500*sin(3.1415927*88*x/lO. 00000000) 

N = 89 

A Term = 0.01487869741345353900*cos(3.1415927*89*x/10. 00000000) 

B Term = 0.00608046246394765980*sin(3.1415927*89*x/10. 00000000) 

N = 90 

A Term = -0.00397363565827067460*cos(3.1415927*90*x/10. 00000000) 

B Term = 0.01888951788189512100*sin(3.1415927*90*x/10. 00000000) 

N = 91 

A Term = -0.00533687574762478110*cos(3 .1415927*91*x/10. 00000000) 

B Term = -0.00939670010485410550*sin(3.1415927*91*x/10. 00000000) 

N = 92 

A Term = 0.00463043960242224850*cos(3.1415927*92*x/10. 00000000) 

B Term = -0.00846707235118616610*sin(3.1415927*92*x/10. 00000000) 

N = 93 

A Term = 0.01110792542036384900*cos(3 .1415927*93*x/10. 00000000) 

B Term = -0.00380977408383048000*sin(3.1415927*93*x/10. 00000000) 
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N = 98 
A Term = 
B Term 
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